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Theoretical  and  experimental  investigations  of  a  novel  gas  dyn¬ 
amics  device  having  no  moving  parts  yet  performing  the  function  of  a 
compressor,  are  described.  This  device,  called  the  "Aerothermopressor", 
exploits  the  possibility  of  raising  the  total  pressure  of  a  high  speed 
gas  stream  through  cooling  of  the  gas.  VIhen  placed  at  the  exhaust  of 
a  gas  turbine,  the  Aerothermopressor  will  reduce  the  exhaust  pressure, 
thereby  improving  both  fuel  economy  and  power  capacity  per  unit  of  air 
flow,  Basic  elements  of  the  apparatus  comprise  a  nozzle  which  acceler¬ 
ates  hot  gas  into  an  evaporation  section;  a  water  injection  system  which 
delivers  finely-atomized  water  into  the  high-speed  stream;  an  evapora¬ 
tion  section  in  which  the  gas  is  cooled  and  most  of  the  water  evaporated; 
and,  finally-jS  diffuser  in  which  the  gas  stream  is  decelerated  and  the 
static  pressure  increased. 

Although  the  Aerothermopressor  is  simple  in  structural  arrangement, 
the  physical  processes  occurring  within  it  are  exceedingly  complex  in  their 
details.  The  simultaneous  effects  on  the  gas  stream  of  droplet  drag,  evap¬ 
orative  cooling,  area  variation  and  wall  friction  lead  to  many  regimes  of 
operation,  including  the  hitherto  unknown  passage  from  subsonic  to  super¬ 
sonic  speeds  in  a  constant-area  duct.  Theoretical  calculations  of  a  one- 
dimensional  nature,  involving  for  the  gas  stream  the  equations  of  continuity, 
momentum,  and  energy,  and  for  the  liquid  droplet  cloud  the  equations  of 
motion,  heat  transfer,  and  mass  transfer,  have  been  carried  out  on  a  high¬ 
speed,  electronic  digital  computer.  The  theory  reproduces  all  the  behavior 


patterns  of  experimontal  units  and  is  in  generally  good  quantitative 
agi'eement  with  the  experimental  data. 

The  results  of  exporjjnents  on  a  small-scale,  constant-area  unit  of 
2,13  inches  diameter  are  presented  and  compared  with  theoretical  calcula¬ 
tions.  The  experiments  and  theory  both  show  that  a  net  stagnation  rise 
is  possible  only  with  gas  flows  greater  than  about  2  Ib/sec;  below  this 
value  the  detrimental  effects  of  wall  friction  completely  absorb  the  gains 
due  to  cooling.  In  the  range  of  25  Ib/sec.  a  net  stagnation  pressuie 
rise  of  about  lOjS  seems  assured,  while  20%  seems  possible, 

Eai’ly  tests  on  a  recently-completed  medium-scale  unit  of  25  Ib/sec 
capacity  have  already  demonstrated  a  net  over -all  rise  in  stagnation  pres¬ 
sure. 

2.  KOMEMCL/tTURE 
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see  Eq.  8b 

specific  enthalpy  of  air 

specific  enthalpy  ''of  liquid  water  in  droplet 

specific  enthalpy  of  water  vapor 

film  coefficient  of  mass  transfer  (Bq.  0-2) 

film  coefficient  of  heat  transfer  (Bq.  C-5s) 

stagnation  enthalpy  of  gas  phase 

stagnation  enthalpy  of  air 

stagnation  enthal^r  of  water  vapor 

ratio  of  specific  heats  for  gas  phase 

length  of  evaporation  section 

latent  heat  of  water  vapor  at  temperature  7^ 

Mach  Number  of  gas  phase,  V/c 

Nusselt  Number  for  mass  transfer,  y,/8 

Nusselt  Number  for  heat  transfer, 

static  press\ire  of  gas  phase 

stagnation  pressure  of  gas  phase 

stagnation  pressure  of  mixture  (Bq.  A-15) 


Prandtl  Number^ 

heat  added  per  unit  mass  of  gas 
gas  constant  of  gas  phase 
gas  constant  of  air 
gas  constant  of  water  vapor 
universal  gas  constant 
relative  Reynolds  Number  of  droplet, 
Schmidt  Number^ 
time 
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absolute  temperature  of  gas  phase 

temperature  of  droplet 

stagnation  temperatwe  of  gas  phase 

gas  velocity 

droplet  velocity 

mass  rate  of  gas  flow 

mass  rate  of  air  flow 

mass  rate  of  water  flow 

mass  rate  of  water  injected 

molecular  weight  of  gas  phase 

molecular  weight  of  air 


molecular  weight  of  water  vapor 


fraction  evaporated. 


u/a^ 


see  Eq.  13 

longitudinal  distance  from  inlet  plane 

see  Eq.  D-8 

e/r  (see  Eq.  A-13) 

See  Table  III 

diffuser  loss  coefficient  (see  Eq.  D-15) 
compressor  efficiency 
turbine  efficiency 

mixture  stagnation  temperature  (see  Eq.  A-13) 


thermal  conductivity  of  gas  phase 
viscosity  of  gas  phase 
mar  t  density  of  gas  phase 
mass  density  of  water  in  droplet 

spatial  mass  density  of  saturated  water  vapor  at  droplet  surface 
spatial  mass  density  of  water  .vapcr  far  from  droplet  surface 
shear  stress  at  pipe  wall 

specific  humidity  of  gas  phase,  lb.  of  water  per  lb.  of  air. 
Initial  water-air  ratio, 

Subscripts 

/  at  inlet  of  evaporation  section 

Z  at  exit  of  evaporation  section 

3  at  exit  of  diffuser 
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£  water 

r  water  vapor 


3.  INTRODUCTION 


3.1  OMect 

•i't 

The  purpose  of  this  paper  is  to  present  both  an  introduction  and 
a  progress  report  on  a  novel  aerothermodynamic  device  which  performs  the 
function  of  a  compressor  but  which  requires  only  an  extremely  simple  mech¬ 
anical  structure  having  no  moving  parts.  Basically,  the  Aerothermopressor 
is  a  duct  within  which  atomized  water  evaporates  into  a  high-speed  stream 
of  high-temperature  gas,  thereby  inducing  a  rise  in  isentropic  stagnation 
pressure  of  the  gas  stream.  One  of  the  most  attractive  applications  of 
the  Aerothermopressor  Is  as  an  auxiliary  for  Improving  the  performance  of 
a  gas  turbine  plant. 

The  plan  of  the  paper  is  (1)  to  develop  the  basic  ideas  and  the  theory 
of  the  Asrotherroopressor,  (2)  to  illustrate  its  application  to  the  gas  tur¬ 
bine  plant,  (3)  to  discuss  the  factors  entering  into  the  design  of  an  Aero¬ 
thermopressor  and  the  considerations  for  handling  them,  (4)  to  present  the 
theoretical  and  experimental  results  thus  far  obtained,  and  (5)  to  evaluate 
the  future  prospects  of  the  Aerothermopressor  and  name  the  crucial  problems 
which  it  faces. 

3.2  Basic  Concent  and  Preliminary  Theory  of  Aerothermopressor 

A  glance  at  the  history  of  typical  inventions  shows  that  almost  with¬ 
out  exception  the  theory  of  a  device  follows  the  conception  of  its  mechan¬ 
ical  arrangement.  The  Aerothermopressor  seems  to  be  one  of  those  rare 
instances  wherein  a  theoretical  analysis  revealed  the  practical  possibil¬ 
ities;  in  fact,  the  mode  of  operation  of  the  Aerothermopressor  is  so  far 
removed  from  that  intuitive  sense  which  usually  underlies  invention,  that 


one  can  hardly  Imagine  lie  having  been  rationally  conceived  prior  to 
the  theoretical  advances  in  gas  dynamics  of  the  past  decade. 

]>l«A]]1n^Tl«rv  Analysis.  The  formula  which  motivated  the  develops 
ments  reported  herei  presented  by  ^apiro  and  Hawthorne  (ref.  1)  in 
1947, indicates  how  the  isenteopic  stagnation  pressure  of  a  gas  stream 
(Fig.  1)  is  affected  by  the  various  external  influences  which  may  alter 
the  state  of  the  stream! 


(Eq.l) 


This  equation,  which  is  derived  from  the  continuity,  momentum,  and 
energy  equations  for  a  perfect  gas,  is  the  result  of  a  one-dimensional 
analysis  of  a  simple  model  in  which  the  external  influences  considered 
Include  changes  in  the  stagnation  temperature,  pipe-wall  friction,  and- 
the  stepwise  injection  and  evaporation  of  liquid  with  a  concomitant  change 
in  molecular  weight  of  the  gas  phase.  The  quantities  appearing  in  this 
one-dimensional  analysis  may  be  thought  of  as  representing  certain  average 
properties  of  the  actual  duct  flow.  The  local  stagnation  pressure  ^ 
is  the  pressure  which  the  stream  would  reach  if  it  were  isentropically 
decelerated  to  rest  in  steady  flowj  similarly,  the  local  stagnation  temper¬ 


ature  T.  is  the  absolute  temperature  which  the  stream  would  reach  if 
it  were  adiabatically  decelerated  to  rest  in  steady  flow.  The  symbol  k 
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denotes  the  ratio  of  specific  heats;  M,  the  local  Kach  Number;  -f^the  pipe- 
wall  skin-ft-iction  coefficient;  D,  the  local  hydraulic  diameter  of  the 
duct;  longitudinal  distance  along  the  duct;  the  local  mass  rate 
of  flow  of  the  gas  phase;  K  the  local  molecular  weight  of  the  gas  phase; 
and  ratio  of  forward  component  of  liquid  injection  velcjcity  to 

gas  velocity.,  In  each  infinitesimal  interval  di  ,  the  mass  dlO'  of 
liquid  is  assumed  to  be  injected  into  the  gas  per  unit  time,  and  to  be 
completely  evaporated  and  mixed  at  the  end  of  the  interval. 

What  Eq.  1  shows  clearly  is  that  the  stagnation  pressvu-e  tends  to 
be  reduced  by  wall  friction,  by  the  aerodynamic  drag  associated  with  the 
injection  of  liquid  (assuming  ^  ^  ),  and  by  a  reduc¬ 

tion  in  molecular  weight.  Changes  in  cross-sectional  area  do  not  in 
themselves  alter  the  local  rate  of  change  of  stagnation  pressure.  Most 
significant  for  our  purpose,  however,  is  the  observation  that  a  reduction 
in  stagnation  temperature  {l.e.,  cooling  of  the  gas  flow)  tends  to  Increase 
the  stagnation  pressure.  Moreover,  Eq.  1  indicates  that  the  effects  noted 
are  the  result  of  a  dynamic  process,  as  evidenced  by  the  proportionality 
of  the  stagnation  pressure  change  to  the  square  of  the  I^ach  Number,  and 
that  they  may  not  be  produced  by  thermodynamic  changes  in  a  static  system. 
To  obtain  significant  changes  in  ^  ,  therefore,  it  is  evident,  even 

from  this  preliminai'y  analysis,  that  operation  at  high  subsonic  speeds  or 
supersonic  speeds  may  be  necessary. 

The  only  possibility  for  increasing  the  stagnation  pressure  is  through 
a  process  in  which  the  decrease  in  stagnation  temperature  can  be  accom¬ 
plished  in  such  a  way  that  the  adverse  effects  seen  in  Eq.  1  are  outweighed. 


Two  practical  poasibilities  come  to  mind!  the  first  involves  the 
use  of  ordinary  heat  exchangers,  and  the  second  involves  the  evapora¬ 
tion  of  a  volatile  liquid  into  the  air  stream. 

Futility  of  Feat  Exchangers.  For  an  ordinary  heat  exchanger  in 
which  the  air  flows  through  a  tube  having  cooled  walls,  Eq.  1  is  sim-  . 
plified  to  the  form 

n  2  L  X 


+  4i 


(Eq.  2) 


Detailed  analysis  of  this  equation  (ref.  1),  enqploying  Reynolds'  anal- 
ogy  between  friction  and  heat  transfer,  shows  that  the  stagnation  pres¬ 
sure  loss  produced  by  wall  friction  always  exceeds  the  stagnation  pres¬ 
sure  rise  attainable  by  cooling.  Consequently  heat  exchangers  can  not 

accomplish  the  desired  resxilt. 

CnoTlng.  In  order  to  determine  whether  evaporative 
cooling  offers  any  prospects  whatsoever,  the  analysis  will  for  this  pur¬ 
pose  be  simplified  in  a  favorable  way  by  assuming  that  the  frictional 
term  of  Eq.  1  is  negligible  compared  with  the  other  terms,  and,  conserva¬ 
tively,  that  the  liquid  is  injected  with  negligible  forward  velocity.  Then 
Eq,  1  may  be  written 


^  4. 

ft  2  Li 


(Eq.  3) 


The  energ7  equation  may  be  written  approximately,  ignoring  the 
kinetic  energy  increase  of  the  injected  fluid,  by  equating  the  de¬ 
crease  of  stagnation  enthalpy  of  the  gas  stream  to  the  enthalpy  rise 
of  the  evaporated  liquid: 


where  Cp  is  the  specific  heat  at  constant  pressure  of  the  gas  phase, 
is  the  enthalpy  of  the  injected  liquid  after  it  has  been  evapor¬ 
ated  and  brought  to  thermal  equilibrium  uith  the  gas,  and  ^  is  the 


enthalpy  of  the  Injected  IdqUid. 

fbr  the  change  irl  molecular  weight,  we  find,  from  the  definition  of 
molecular  weight  for  a  gas  mixture,  the  formula 


W  V 


(Eq.  5) 


where  is  the  molecular  weight  of  the  evaporated  liquid. 

Substituting  Eqs.  4  5  ''®  obtain 


l-W/Wy 

77^-^ 


lolw" 


(Eq.  6) 


FroD  this  relation  it  may  be  seen  that  the  criterion  for  a  rise  in 
stagnation  pressure  is  that 


WMr  -  i 


>  2 


(Bq.  7) 


Now,  other  things  being  the  sane,  liquids  with  large  latent  heats 
are  most  likely  to  meet  the  criterion  of  Bq.  7,  inasmuch  as 
is  of  the  order  of  magnitude  of  the  latent  heat.  Hence  water  seems  to 
be  the  desirable  choice  for  the  Aerothemopressor  because,  apart  from 
obvious  economic  reasons,  it  has  a  larger  latent  heat  than  almost  any 
other  fluid.  To  gain  some  concept  of  orders  of  magnitude.  Table  I  sets 
out  approxiimte  values  of  the  left-hand  side  of  Eq.  7  for  evaporation 
of  liquid  water  into  air. 


TABLE  I 


Tjj(deg.R) 

M=0 

H=2 

m 

500 

7.4 

7.4 

7.1 

6.8 

6.8 

1000 

4*4 

4.3 

4.0 

3.7 

3.6 

2000 

2.9 

2.7 

2.4 

2.1 

2.1 

Comparing  the  figures  in  TaM.e  I  with  the  criterion  of  Eq.  7, 
it  may  be  concluded  that  there  is  indeed  a  possibility  for  increas¬ 
ing  the  stagnation  pressure  by  evaporating  water  into  a  high-speed 
stream  of  air,  and  that  this  possibility  exists  over  a  wide  range 
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of  stagnation  temperature  and  Mach  Number. 

lnt«Tretlng  «ls  I  it  is  veU  to  recll  (a)  that  water  can 

te  evaporated  into  the  air  only  ao  lone  aa  the  latter  ia  not  saturated 
with  water  vapor,  and  (h)  that  the  amount  of  evaporation  which  can  he 
effected  per  unit  length  of  duct  depends  on  the  driving  temperature  dif- 
fereroe  essenti-a  for  heat  transfer  between  the  gas  and  the  liquid.  These 
considerations,  when  combined  with  the  values  of  Table  I  and  the  form  of 
EC  1  (including  the  wall-friction  tern),  denonstrate  that  the  ma:dmum 
rise  in  sugnation  pressure  will  be  obtained  when  the  stagnation  temper¬ 
ature  is  neither  e:coesslvely  high  nor  excessively  low.  and  with  Ifech 
tabers  which  are  neither  excessively  large  nor  excessively  smell. 
i-i  tnullcation  of  AerothprmaOTssffiJz^-SSS-TiaMnO-Bai 

Of  the  various  applications  of  the  Aerothermopressor  which  have  been 
considered,  its  uee  in  improving  both  the  perfcrnanoe  and  the  chsrscterie- 
tics  of  a  gas  turbine  plant  seems  to  be  one  of  the  most  premising  as  weU 

as  one  of  reasonably  immediate  application. 

The  «t  oh, ions  place  for  the  Aerothermopressor 

is  at  the  e.dx.ust  of  the  power  turbine  in  the  simple  gas  turbins  cycle, 
(Fig.  2a),  where  the  Aerothermopressor  is  supiplled  (at  no  expense, 
were)  with  a  stream  of  hot  air  at  sufflclenUy  high  temperature  to  mehe 
■  Significant  increases  in  stagnation  pressure  possible.  If  the  Aerothermo¬ 
pressor  produces  s  rise  in  stagnation  press, u-e  and  also  discharges  to  the 
atmosphere,  it  follows  that  the  stagnation  pressure  at  the  turbine  exit 
will  be  reduced  below  atmospheric,  thus  msHng  the  Aerothermopressor  a„^- 
Ogous  to  the  condenser  in  a  steam  power  plant. 


A  second  cycle  arrangement  is  sho'^/n  in  Fig.  2b •  Here  more  fuel 
is  burned  per  unit  of  air  flow  than  in  the  cycle  of  Fig»  2a,  the  com¬ 
pensating  advantages  being  that  the  pressure  at  the  tia’bine  inlet  is 
increased  while  at  the  same  time  the  mass  rate  of  flow  through  the 
turbine  is  augmented.  Nothing  stands  in  the  way  of  installing  a  second 
Aoi’othermopressor  at  the  turbine  exhaust,  as  suggested  by  the  dashed 
lines  of  Fig.  2b,  thus  compounding  the  improvements  in  performance. 

Effects  on  Performance.  In  Figs.  3b  and  3c  are  shown  tho  results  of 
cycle  caloilations  for  the  simple  cycle  of  Figs.  2a  and  3a,  illustrating 
the  advanttges  mentioned  above. 

It  may  be  seen  that  the  reduction  of  back  pressure  betters  the  per¬ 
formance  and  the  chax’acteristlcs  of  tho  gas  turbine  plant  in  a  number 
of  ways;  (a)  the  specific  fuel  consumption  is  reduced;  (b)  the  specific 
air  consumption  is  decreased,  i.e.,  tho  power  capacity  of  a  given  sisio 
of  machine  is  increased;  (c)  the  optimum  compressor  ratio  is  reducod  in 
inverse  proportion  to  the  square  root  of  the  pressure  ratio  across  the 
Acrothermopressor ;  and  (d)  tho  net  work  ratio  is  increased,  thus 
the  plant  perfoi’mance  less  sensitive  to  such  design  variables  as  turb¬ 
ine  inlet  temperature  and  compressor  and  turbine  cfficicncisn,  and  con¬ 
sequently  permitting  acceptable  performance  with  less  stringent  opeciti- 
cations  on  tliese  variables  than  would  otlierwise  be  possible. 


With  a  high  tiu’bine  inleu  teiriperatuvc  and  with  high  Vtilues  of  tu; 
and  compressor  efficiencies  (Fig.  3b),  the  percentage  improveraein^s  in 
economy  and  in  power  capacity  aro  approxircatGly  equal  to  the  per  sort  ■' 
rise  in  stagnation  pressure  across  tho  Aoroiliennoprosscr .  For  can'.,;: 


With  a  co.«pr6S3or  pressure  ratio  of  4,  .  stagnation  pressure  ratio 
across  the  AerothersK^pressor  of  1.2  produces  e  reduction  in  specific 
fuel  consineption  and  specific  air  consumption  of  20^. 

With  »  low  turbine  inlet  temperature  and  lou  component  efficiencies 
(Kg.  30),  the  improvements  in  performance  are  even  more  marked,  end 
mm/well  change  the  situation  from  one  of  marginal  net  power  output  to 

one  of  reu:?onably  acceptable  performance. 

rc,nenfsnn  nf  AerothMloECirPh  Rogoo^rafie^  In  the  embodiment 

of  Mgs.  2a  and  3a,  the  Aerothormopressor  must  compete  with  the  regenera¬ 
tor,  the  latter  being  the  conventional  scheme  for  utilising  the  energy 
in  the  high-temperature  exhaust  gas.  At  this  writing  it  does  not  seem 
that  the  Aerothermopressor  con  match  the  improvement  in  fuel  economy 
accessible  to  the  regenerator.  On  the  other  hand,  the  Aerothermopressor 
brings  about  a  reduction  in  specific  air  consumption  (l.e..  a  rise  in  power 
capacity)  of  magnitude  equal  to  the  reduction  in  spaoifio  fuel  consump¬ 
tion,  whereas  the  regenerator,  by  virtue  of  its  frictional  pressure  drop, 
produces  a  slight  loss  in  power  capacity.  If., at  may  he  even  more  import¬ 
ant  in  ..any  .ractlcal  sppllcations  is  that  an  effeotive  regenerator  rots 
the  gas  turbine  plant  of  the  co,.pactness  and  lightness  which  often  is  the 
outstanding  claim  of  the  gas  turbine  over  rival  prime  movers.  ,is  will  be 
seen,  the  Aerothermopressor  is  sufficiently  ...mall  and  light  as  not  eo  bo 
at  a  disadvantage  In  this  respect,  except  perhaps  for  automotive  applica¬ 
tions.  A  final  compensating  advantage  of  the  Aerothermopressor  is  that 
it  is  structurally  simple  and  not  subject  to  the  serious  problems  of 


fouling  and  corrosion  which  blight  the  status  of  the  regenerator. 

?■/.  meet  of  Size  on  Aerothermopressor.  FerfomilSg 

So  important  is  the  scale  effect  in  the  Aerothermopressor  that  it 
seems  desirable  to  conclude  this  introduction  with  an  explanation  of 
this  phenomenon  before  proceeding  to  more  detailed  considerations. 

For  given  gas  properties  at  the  beginning  of  the  evaporation  section 
and  a  given  initial  droplet  cloud,  the  necessary  length  of  the  evapora¬ 
tion  section  is  roughly  established  by  the'  residence  time  required  for 
the  droplets  substantially  to  be  evaporated.  If  the  gas  and  water  flows 
are  quadrupled,  for  example,  the  oross-sectional  area  of  the  duct  will. also 
be  quadrupled,  but  such  gas  properties  as  velocity,  pressure,  temperature, 
etc.,  will  at  corresponding  sections  be  nearly  the  same  as  before.  Accord¬ 
ingly,  the  required  residence  time  and  the  required  length  will  also  be 
unchanged,  but  the  length-diameter  ratio  of  the  duct  will  be  only  half 
as  large  as  in  the  first  Instance.  Since  the  stagnation  pressure  loss 
associated  with  pipe  ftiction  is  roughly  proportional  to  the  length-diameter 
ratio  (Eq.  1),  the  net  rise  in  stagnation  pressure  will  be  greater  for  the 
case  with  the  larger  air  flow.  This  Influence  of  air  flow  on  the  relative 
prominence  of  frictional  effects  is  furthermore  magnified  because  the  net 
rise  in  stagnation  pressure  of  Eq.  1  turns  out  to  be  a  comparatively 
small  difference  between  comparatively  large  quantities. 

What  all  this  amounts  to  is  that  below  a  certain  scale  of  air  flow 
the  Aerothermopressor  is  doomed  to  failure;  there  is  a  critical  size  for 
which  the  Aerothermopressor  will  be  on  the  verge  of  producing  a  rise  in 
stagnation  pressure;  and  above  this  critical  size  (which,  according  to  our 
present  knowledge,  is  in  the  diameter  range  of  2  to  3  inches )the  attainable 
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stagnation  presanre  rise  >nounts  rapidly  as  the  scale  of  air  flow  Is 
increased s 

.  rpHp.nnv  ns-  AmOTHERMOPREggOR 

T  TViftrmodynanilc  Llruit  of  Performance 

One  of  the  first  .prestlons  that  cones  to  mind  In  evaluating  a 
new  device  like  the  Aerothermopressor  Is,  what  la  the  naxMom  perl 
ance  of  whloh  It  Is  theoretically  oapable7  This  nay  he  answered  on 
pnnely  thernodynanlc  grounds  hy  considering  a  steady-flow  process  In 
which  the  streams  of  Inlet  air  and  Inlet  water  are  brought  together 
in  a  hypothetical  apparatus  having  no  themodynanlc  Irrevsrsihllltles, 
and  producing  an  exit  strsan  of  air  saturated  with  water  vapor.  As  In  the 
aotual  Aerothemopreasor,  It  Is  asauncd  that  the  reversible  apparatus 

does  not  exchange  heat  or  work  with  the  surroundings. 

The  analysis,  for  given  Initial  properties  of  air  and  water  requires 
that  (1)  the  stagnation  enthalpy  increase  for  the  water  ecpral  the  corres¬ 
ponding  decrease  for  the  air.  (11)  the  entropy  rise  of  the  water  equal  the 
corresponding  decrease  for  the  air,  and  (ill)  the  exit  atream 
ated.  Pertinent  results  of  the  calculations,  are  shown  in  Table  II 

for  a  variety  of  inlet  conditionse 

iThTISSySiV  and  ealculationrS^JTSde  hy  Mr.  Alve  Erlekson, 

Instructor  in  Mechanical  Englneoring.  Masoachusotts  Institute  of 


Technology. 


TABLE  II 


THKRMDDYNAMIC  LIMITS  OF  PERFORMANCE 


Inlet  Conditions 


Air 

Water 

X 

deg.F 

t 

psia 

T 

w 

deg.F. 

psia 

740 

10.0 

60 

U.7 

n 

U.7 

It 

It 

1040 

10.0 

ti 

II 

II 

12.0 

n 

It 

II 

14.7 

II 

n 

1340 

10.0 

it 

tf 

H 

14.7 

It 

It 

1040 

U.7 

60 

U.7 

It 

n 

100 

n 

1040 

U.7 

60 

U.7 

n 

It 

11 

3000 

Water/Alr 
Ratio  ■ 


Outlet  Conditions 
(Air  Saturated) 


T 

deg.F 


fi 

psla 


0.128 

0.125 

0.185 

0.183 

0.181 

0.2a 

0.234 

0.181 

0.187 

0.181 

0.182 


166 

180 

2U 

222 

230 

267 

284 

230 


230  89.8 


230 

231 


32.0 

45.0 

67.0 

78.6 

92.3 

142.5 

190.1 

92.3 


92.3 

93.6 


3.2 

3.1 

6.7 

6.6 

6.3 

U.3 

12.9 

6.3 

6.1 

6.3 

6.4 
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The  overall  stagnalion-presswe  ratio  of  the  gas  stream  (shovm  in  the 
last  column)  is  seen  to  depend  primarily  on  the  inlet  stagnation  temper¬ 
ature  of  the  air.  Variations  in  the  inlet  air  pressure,  inlet  water  pres¬ 
sure,  and  inlet  water  tempera tui'e  exert  relatively  small  influence  on 
the  overall  stagnation-pressm'e  ratio. 

So  attractive  are  the  figures  in  the  last  column  of  Table  II  that  a 
word  of  caution  is  necessary  lest  they  be  given  too  much  weight.  More  de¬ 
tailed  investigation  discloses  that  it  does  not  seem  possible  in  a  real 
Aerothermopressor  even  to  approach  the  reversible  mixing  of  air  and  water 
upon  which  Table  II  is  predicted.  The  main  sources  of  irreversibility  in 
the  Aerothermopressor  are  (i)  wall  friction,  (ii)  aerodynamic  drag  on  the 
accelerating  droplets,  (iii)  heat  transfer  from  the  aii’  to  the  droplets 
across  a  finite  temperature  difference,  and  (iv)  mass  transfer  of  water 
vapor  from  the  droplets  to  the  air  stream  across  a  finite  difference  in 
partial  pressure.  As  later  discussions  of  the  detailed  mechaniE’-ms  of  the 
Aerothermopressor  will  show,  there  seems  to  be  no  means  by  which  a  device 
having  the  general  character  of  the  Aerothermopressor  can  sjjaultanoously 
reduce  all  those  irreversibilities  to  small  proportions. 

A  corollary  of  these  remarks  is  that  the  last  coluni!!  oi  Table  XI, 
while  giving  absolute  thermodynamic  limits,  does  not  give  a  realistic 
view  of  what  the  actual  celling  to  Aerothermopressor  performance  might 
be.  Thus  we  are  presently  in  the  irritating  status  of  not  really  know¬ 
ing  what  is  the  best  performance  which  might  be  reasonably  expected  of 
the  Aerothermopressor,  and  it  appears  that  detailed  calculation^  like 


those  of  Section  4»7  may  be  necessary  for  Its  prediction. 

Still  a  further  corollary  Is  that  the  effects  of  inlet  air  and 
water  temperature  and  of  inlet  pressure  shown  in  Table  II  are  probably 
not  representative  of  their  effects  in  an  actual  Aerothermopressor.  In¬ 
asmuch  as  irreversibilities  play  such  a  strong  role  in  the  Aerotherrao- 
pressor,  the  influence  of  the  aforementioned  properties  might  be  realized 
most  forcefully  through  their  effects  on  the  Irreversibilities,  and  this 
of  course  is  an  aspect  totally  absent  from  the  analysis  underlying  Table  II. 
4.2  One-Dimensional  Analysis  of  Aerothermopressor 

Because  wall  friction  dictates  a  small  length-diameter  ratio  for  the 
Aerothermopressor,  the  continuously  distributed  injection  and  evaporation 
suggested  by  Fig.  1  does  not  represent  a  good  design.  In  order  to  approx¬ 
imate  TM-Himiin  droplet  residence  times,  it  seems  desirable  for  most  of  the 
water  to  be  injected  near  the  beginning  of  the  evaporation  section,  as 
suggested  in  Fig.  2a.  For  a  real  Aerothermopressor,  therefore,  Eqs.  6 
and  7  do  not  adequately  represent  what  occurs,  even  though  they  are  surely 
indicative  of  the  true  state  of  affairs  inasmuch  as  the  analysis  on  which 
they  are  based  does  in  fact  take  account  (albeit  according  to  a  somewhat 
distorted  model)  of  the  principle  physical  phenomena. 

peflnltlon  of  One-Dimensional  Model.  In  order  to  obtain  a  better  in¬ 
sight  into  the  details  of  the  Aerothermopressor  process,  however,  and  also 
to  provide  a  more  precise  theory  for  numerical  calculations,  it  is  necessary 
to  develop  an  analysis  which  is  in  better  accord  with  the  true  events.  In 
the  model  we  now  adopt,  the  flow  is  treated  as  one-dimensional  (i.e,,  no 
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radial  variations  in  any  properties  of  either  the  gas  stream  or  the  drop¬ 
let  cloud),  and  the  droplets  are  assumed  to  be  of  uniform  size.  Since 
atomization  occurs  very  quickly,  it  is  further  assumed  that  at  the  inlet 
plane  of  the  evaporation  section,  spherical  droplets  of  equal  size  are 
uniformly  distributed  over  the  cross-section  and  have  the  velocity  and 
temperature  with  which  they  left  the  injection  nozzles. 

At  any  section  of  cross-sectional  area  A,  the  gas  flows  at  the  mass 
rate  W  and  has  a  certain  velocity  V  ,  pressure  p  ,  temperatxire  T  , 
and  composition  cO  .  Simultaneously,  the  droplet  cloud  has  a  certain 
droplet  diameter  <L  and  its  own  mass  rate  of  flow  ,  velocity  , 
and  teraperatOTe  ^  distance  further  downstream,  the  afore¬ 

mentioned  properties  of  both  the  gas  stream  and  droplet  cloud  have  been 
altered  in  consequence  of  the  combined  effects  of  wall  friction,  change 
in  cross-sectional  area,  droplet  drag,  heat  transfer,  and  evaporation: 
the  latter  three  phenomena  result  respectively  from  the  differences  between 
the  two  phases  in  velocity,  temperature,  and  partial  pressure  of  the  water 
vapor. 

Influence  Coefficients,  For  the  infinitesimal  interval  dz  (Fig,4a), 
the  analysis  is  based  upon  a  simultaneous  solution  of  the  following  equa¬ 
tions  written  in  differential  form;  (i)  the  governing  equations  repres¬ 
enting  conservation  of  mass,  the  momentum  theorem,  and  the  first  law  of 
thermodynamics,  (li)  the  equation  of  state  of  a  perfect  gas  mixture  hav¬ 
ing  variable  specific  heat  and  obeying  the  Gibbs-Dalton  law  of  partial 
pressures  and  enthalpies,  and  (iii)  the  definitions  of  Mach  Number,  gas-phase 
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stagnation  tenperature  and  pressure^  and  nixtura  stagnation  terip- 
erature  and  pressure*  Details  of  the  analysis  are  given  in  Appendix  Af 
and  the  results  are  summarized  by  the  system  of  Influence  coefficients 
presented  in  Table  III.  The  influence  coefficients  are  sin^^ly  the  co¬ 
efficients  in  the  algebraic  relations  connecting  the  variables  of  the 
left-hand  column  vith  the  variables  of  the  top  row.  For  example, 

dM*_  dA  I  dQ.-^d\  . 

"  1  -  M*  A  1  -  M*  Cp\ 

and  so  on. 

Behavior  of  Stream  Properties  Iteder  the  Influence  of  Area  Change^ 

Heat  'I’ransfer.  Evaporation.  Wall  Frictiopf  and  Droplet  Acceleration «.  Apart 
from  their  utility  for  numerical  calculations,  the  influence  coefficients 
of  Table  III  reveal  clearly  how  the  stream  behaves  under  the  action  of 
area  change,  ftriction,  etc.  A  discussion  of  these  influence  coefficients 
at  this  point  will  help  later  to  explain  the  complex  and  varied  behavior 
patterns  ejdiiblted  by  the  Aerothermopressor  as  well  as  point  the  way  to 
rules  for  design. 

Detailed  investigation  of  orders  of  magnitude  shows  that  of  the 
variables  at  the  top  of  Table  III,  those  which  predominate  in  the  actual 
Aerothermopressor  process  are  the  terms  involving  area  change  (  dA/A  ), 
evaporation  {(hjJ'fw'  ),  wall  friction  (  4idi/D  ),  and  droplet  accel¬ 
eration  (  cl\^  y  ) .  How  these  effects  tend  to  change  the  stream 
properties  is  summarized  In  Table  IV.  Although  the  table  is  self-explanatory, 
some  additional  remarks  will  emphasire  the  moot  important  results. 


TABLE  m  -  INFLUENCE  COEPFICIENTS 


TABLE  IV  -  EEHAVIOR  OF  STREAM  PROPERTIilS 
UNDER  INFLUENCE  OF  AREA  CHANGE^  EVAPOR- 
ATTOW,  NALL  FRICTION.  AND  DROPLET  DRAG 


Area 

Evaporation 

Wall 

Liquid 

Increase 

produces 

,  (b) 

produces 

friction,  . 
produces'®' 

acceleration 

produces 

Mach  ^^u^lber,M 

subsonic 

decrease 

decrease 

increase 

(e) 

increase '  ' 

decrease 

supersonic 

increase 

\ 

,  (h) 

increase 

decrease 

Gas  Velocity, V 

subsonic 

decrease 

decrease 

.  (h) 

Incraase 

increase 

(o) 

increase 

decrease 

supersonic 

increase 

decrease 

Pressure,  p 

subsonic 

supersonic 

increase 

decrease 

]  (h) 

increase' 

decrease 

decrease 

Increase 

decrease 

increase^®) 

Temperatiure ,  T 

subsonic 

increase 

(h) 

decrease' 

increase 

decrease 

decrease 

.  (e) 

increase 

supersonic 

decrease 

Increase 

Gas  Stagnation 

subsonic 

nil 

decrease 

nil 

decrease'  ' 

(f) 

decrease 

Temperature,!^ 

supersonic 

nil 

decrease 

nil 

Mixture  Stagna- 

subsonic 

nil 

decrease 

nil 

nil 

tlon  Temperature, 
00 

supersonic 

nil 

decrease 

nil 

nil^e) 

Gas  Stagnation 

subsonic 

nil 

.  (h) 

increase' 

decrease 

decrease 

Pressure, 

supersonic 

nil 

.  (h) 

increase 

decrease 

decrease'  ' 

Mixture  Stagna¬ 
tion  Pressure, Pq 

subsonic 

nil 

,  (h) 

increase 

decrease 

(g) 

decrease  ® 

supersonic 

nil 

increase 

decrease 

(g) 

decrease  ® 

Notes:  (s’!  Oonosite  effects  for  area  deereasft. 

(b)  Opposite  effects  for  condensation.. 

(c)  Opposite  effects  are  impossible. 

(d)  When  1 ,  (l\^  >Oj when  y>i^d.Vg'^0 

(e)  Dependent  upon  magnitude  of  y  for  liquid  d^peleratlas.,. 

(f)  Opposite  effect  for  liquid  deceleration. 

(g)  Same  effect  for  liquid  deceleration. 

(h)  Based  on  ^  only,  and  generally  correct  for  in  excess 

two;  otherwise  effects  are  indeterminate. 


of 


Area  Change.  This  has  no  direct  effect  on  any  of  the  local  rates 
of  change  of  the  stagnation  properties.  Its  greatest  importance  is 
that  it  gives  the  designer  some  control  over  the  Mach  Number  and  the 
temperature  of  the  stream. 

HeaL Transfer  and  Chances  in  Liquid  Temperature.  The  term  in  dQ 
represents  external  heat  exchange  between  the  two-phase  sti'eam  and  the 
duct  wall}  while  the  term  in  olXj  represents  the  excess  of  a32  the  heat 
received  by  the  droplet  cloud  over  the  enthalpy  rise  of  the  evaporated 
liquid.  Net  energy  extraction  from  the  gas  stream  either  through  a 
negative  value  of  fiQ  or  through  a  positive  value  of  dfi/^  tends  to 
produce  a  gain  in  stagnation  pressure,  a  reduction  in  stagnation  temper- 
atwe,  and  a  change  in  Mach  Number  away  from  unity  (M  decraaoes  at  sub¬ 
sonic  speeds,  but  Increases  at  supersonic  speeds). 

The  terms  Involving  JtaT  represent  several  phenom¬ 
ena;  (1)  the  energy  exchange  associated  with  the  enthalpy  increase  of 
the  evaporated  liquid,  (ii)  the  energy  exchange  aseoclated  with  changes 
in  kinetic  energy  of  the  evaporated  liquid  as  its  velocity  changes  from 
that  of  the  liquid  to  that  of  the  gas,  and  (iii)  the  momentum  exchange 
associated  with  the  acceleration  of  the  evaporated  liquid.  For  most  inter¬ 
vals  of  ohe  Aerotherraoprossor  process  the  first  two  effects  (embodied  in 
the  term  )  are  by  far  the  most  important,  and  therefore  Table  IV 
shows  the  effects  due  to  the  term  alone.  With  this  assumption,  the 

criterion  for  a  rise  in  ^  is  that  /S  >  Z(\-‘^)  ,  which  in¬ 

dicates  that  the  stream  Is  more  apt  to  rise  in  stagnation  pressure  after 
the  droplet  cloud  has  been  accelerated  than  in  the  period  immediately 
following  injection. 
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Wall,  flrlctlon.  Of  all  the  Independent  variables  in  Table  III, 
the  wall  friction  term  is  the  only  one  which  cannot 

assume  negative  values — that  Is,  the  term  is  by  definition 

always  positive,  and  the  friction  coefficient,  according  to  the  second 
law  of  thermodynamic 8,  may  never  be  negative.  Consequeniily,  both  stag¬ 
nation  pressures  are  reduced  with  any  Increase  in  length,  thus  pointing 
to  the  prime  requirement  of  accomplishing  the  evaporation  in  the  minimum 
possible  distance. 

Chanj;93  in  Liquid  Velocity.  The  term  in  must  be  Inter¬ 

preted  with  the  necessary  condition  that  the  liquid  Is  accelerated  only 
while  it  is  traveling  more  slowly  than  the  gas  (  )>  ®nd  is  decel¬ 
erated  only  while  traveling  more  rapidly  (  ^  >  1  ) ,  It  represents 

two  physical  phenomena;  (i)  The  drag  (or  thrust,  as  the  case  may  be)  of 
the  liquid  droplets  on  the  gas  stream,  and  (ii)  the  work  done  by  the  gas 
oti  the  liquid  when  the  latter  is  accelerated.  By  virtue  of  its  manner 
of  definition,  the  mixture  stagnation  temperattire  is  unaffected  dir¬ 

ectly  by  changes  In  liquid  velocity,  whereas  the  gas  stagnation  temperature 


reflects  the  work  effect,  liquid  acceleration  reducing  X  and  docsloratlon 
incroasing  T]!  .  Liquid  drag  a Iways  decreases  the  gas  stagnation  pres¬ 

sure  g  •  liquid  thrust  (negative  dV^  )  usually,  but  not  always,  increas¬ 
es  g  ,  (the  e;Lception  being  when  ^  is  quite  large  compared  with  unity). 
Of  particular  interest  js  the  fact  that  the  mixture  stagnation  pressure  ^ 
is  altered  by  changes  in  liquid  velocity  only  when  y  is  different  frojn 
unity.  Strangely  enough,  both  liquid  drag  and  liquid  thrust  decrease  ^ 

This  is  so  because  changes  in  ^  are  the  result  of  the  irrovorsibilitios 


in  the  acceleration  or  deceleration  process,  and  the  irreversibility  reduces 
whether  the  aerodynamic  force  is  occasioned  by  the  droplets  moving  faster 
than  the  gas  or  vice  versa.  When  the  acceleration  or  deceleration  is  brought 
about  so  slowly  that  the  droplets  follow  closely  the  speed  of  the  gas 
stream  (y  tending  toward  unity),  the  acceleration  or  deceleration  tends 
to  occur  reversibly,  and  the  loss  in  ^  approaches  zero. 

Molecular  Weight,  As  Indicated  hy  ‘5^  tw«  d\N/w  is  controlled 
by  the  evaporation  term  Ahr/vr  and  by  the  composition  of  the  gas  phase  . 
together  with  the  molecular  weights  of  air  and  water  vapor.  The  evaporation 
of  water  (molecular  weight  18)  into  air  (molecular  weight  29)  reduces  the 
molecular  weight  of  the  gas  phase.  With  dW/W  negative.  Table  III  shows 
that  the  change  in  molecular  weight  induced  by  evaporation  tends  to  decrease 
the  stagnation  temperatures  and  pressures  and  to  drive  the  Mach  Number  to¬ 
ward  unity. 

SEgSj-fis;  Hgat*  since  k  for  water  vapor  is  about  1.3  as  compared 
with  about  1.4  for  air,  evaporation  yields  negative  values  of  dk/k  . 

This  in  turn  tends  to  increase  the  Mach  Number,  decrease  the  stagnation 
temperatures,  and  Increase  the  stagnation  pressures. 

Subsonic  yg.  Supersonic  Speeds.  Examination  of  Tables  III  and  IV  shows 
that,  almost  without  exception,  the  effects  of  heat  transfer,  wall  friction, 
etc.  on  the  Mach  Number,  gas  velocity,  pressure,  and  temperature  are  of 
opposite  signs  at  subsonic  and  supersonic  speeds.  This  is  a  well-known 
phenomenon  of  gas  dynamics,  but  it  is  of  particular  interest  and  import¬ 
ance  here  because  the  multiplicity  of  phenomena  occurring  simultaneously 
in  the  Aerothermopressor  may,  by  any  of  several  combinations,  bring  about 
either  choking  or  continuous  transition  through  the  speed  of  sound  (ref.3), 


More  will  be  eald  about  this  later. 

Regarding  the  four  stagnation  iffoperties,  on  the  other  hand,  the 
effects  observed  at  supersonic  speeds  are  the  same  as  those  at  subsonic 

speeds. 

A.3  De8criIrf■^on  of  AerothennoPTeggOT 

With  the  help  of  the  concepts  developed  in  the  preceding  section, 
it  is  now  possible  to  explain  the  many  types  of  flow  patterns  ejdiibitei 
by  the  Aerothermopressor. 

For  slfflplioity,  we  begin  with  an  Aerothermopressor  having  a  long 
constant-area  evaporation  section  and  with  a  gubafifilfi 
plane  of  water  injection  (Fig.  5).  In  the  entry  plane  the  gas  flows  at 
high  velocity  and  temperature,  and  the  liquid  is  Injected  at  low  velocity 
and  temperature.  The  most  laqportant  phenomena  are  those  of  droplet  drag, 
evaporation,  and  waU  friction.  Although  these  three  effects  of  course 
occur  simultaneously,  each  is  predominant  over  some  portion  of  the  process. 
An  appreciation  of  this  fact  is  essential  for  a  qualitative  understanding 

of  the  AerothermopresBor  operation, 

R«»gimB  1  -  Dronlet  Drag  Predominant.  During  the  droplet  acceleration 

period  which  accompanies  and  follows  atomization,  wall  friction  is  coa^jar- 
atively  negligible  because  the  large  velocity  difference  produces  huge 
accelerations  and  excellent  conditions  for  heat  transfer  to  the  droplet 
cloud.  Although  the  effects  of  heat  transfer  and  evaporation  are  enormous, 
they  are  nevertheless  outweighed  by  the  effects  of  droplet  drag.  Vie  may 
therefore  speak  of  Regime  I  as  the  zone  in  which  drag  is  so  predominant 
that  it  controls  the  qualitative  behavior  of  the  stream.  According  to 


Table  IV,  therefore,  the  Mach  Number  increases,  the  static  pressure 
decreases,  the  stagnation  pressures  decrease,  the  gas  velocity  increases, 
and  the  static  and  stagnation  temperatures  decrease. 

At  first  the  velocity  and  temperature  differences  between  the  gas 
and  liquid  are  so  great  that  they  produce  enormous  rates  of  evaporation, 
droplet  velocity  increase,  gas  temperature  decrease,  and  liquid  temper¬ 
ature  increase.  These  high  rates,  however,  diminish  the  driving  velocity 
and  temperatTore  differentials,  with  the  consequence  that  the  aforementioned 
rates  decrease  as  the  process  proceeds. 

In  the  case  of  the  liquid  temperature,  an  Increase  is  reflected  in 
a  greatly  magnified  evaporation  rate,  which  in  turn  tends  to  cool  the 
droplet  (the  vapor  pressure  varies  approximately  as  the  fifteenth  power 
of  the  absolute  temperature  in  the  applicable  range),  Tlierefore  the 
droplets  rapidly  tend  to  approach  the  wet  bulb  temperature  of  the  gas 
stream,  in  which  state  the  heat  received  from  the  gas  is  equal  to  the 
enthalpy  rise  of  the  evaporated  material.  The  local  wet  bulb  temperature 
does  not  vary  much,  and  so  the  droplets  remain  at  nearly  constant  temp¬ 
erature  (in  the  neighborhood  of  150® F)  for  virtually  the  entire  process. 
Regime  II-Dronlet  Evaporation  Predominant.  As  the  velocity  dif¬ 
ference  between  gas  and  liquid  decreases  (y  ten  ng  toward  unity),  the 
effects  of  droplet  drag  decay  much  more  rapidly  than  the  effects  of 
evaporation,  the  latter  being  substantial  even  for  zero  velocity  differ¬ 
ence.  A  point  is  reached,  therefore,  where  the  evaporative  phenomena 
counterbalance  the  drag  phenomena,  and  there  ensues  a  zone  in  which 
evaporation  predominates  and  controls  the  behavior  of  the  stream. 
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I,  .how. 

BtaUc  «A  stagnation  proa«l«8  lnore.se,  tM  ga.  velocity  decreases, 
snd  the  stagnation  end  etatlc  t«.peraturee  decrease. 

dn  nnexpected  resnlt  1.  Indeed  hy  the  gas  deceleration  In  Rega*  1 1. 
at  point  in  the  vicinity  of  the  tr«.eltlon  fro.  Reglne  I  to  Regl»  . 

the  gas  velocity  1.  dlnlniehed  to  the  prevlonsly-lnoreased  droplet  veloc  y. 

.11,  the  ga.  f»««v  deoelerate.  tecanse  of  .«por.tlon,  the  ll.ld  tends 
to  retdln  It.  ««nt«.,  and  thus  the  U-pId  -eon  find.  ItseU  trave  ng 

.or,  rapidly  than  the  ga.  «ad  conaegnently  being  d«elerated.  for  ^e 

OP  .,!«  n  the  llonld  veloolty  1.  ella.tly  greats  th«.  the  gas  v^f. 
nlth  the  veloolty  spread  being  oonst^atly  n^.  Oarrespondln^f,  t^ 
.o-nt,  ««rgy  ratnm.^  to  th,  ga.  .tr.«.  Indnce  ch«.gee  In  gas  s^» 

properties  opposite  to  those  erperlenoed  In  Regis*  I.  for  esas®l..  e 

nsgatlv,  drag  In  Reglr*  n  assist.  In  the  recovery  of  the  stagnation  pres.nr 
lost  in  Regl-  1-  ^t  should  he  noted  that  these  effeots  ere  proportions 
to  the  Of  11^  tRVOlved,  which  steadily  decreases  as  evaporation 
psoceeds.  ««.»  1...  le  returned  to  the  strea.  In  Regis*  11  than  tahen  fro. 

it  in  Rfigine  I* 

In  the  course  of  Regime  II,  t,ne 
Psap^mA  TTlv.Wa31  Fr1lT^1^r  In 

esap^^atlon  rate  heco*..  ever  snaller  as  the  s^face  area  of  the  droplets 
decreases  and  as  the  velocity  *.1  teaperature  differences  pronotlng  h» 
transfer  ere  eaten  away,  mii^tely  a  stage  Is  reached  where  pipe  frlc 
tlon.  Which  hitherto  has  Inrhed  In  the  heoR^ound  as  a  weeh  contender, 

takes  over  by  default. 
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Regime  III,  where  wall  friction  is  controlling,  Table  IV  indioatei 
that  the  Itach  Number  and  gas  velocity  once  again  rise,  and  the  static 
and  stagnation  pressures  fall# 

Little  need  be  said  about  Regime  III  except  that  it  has  no  place 
in  a  well-designed  Aerothermopressor.  Nothing  is  to  bo  gained,  and  much 
can  be  lost,  by  allowing  the  evaporation  section  to  be  so  long  that  evap¬ 
orative  cooling  is  outweighed  by  pipe  friction.  The  stream  should  be 
allowed  to  enter  a  decelerating  diffuser  some  distance  before  the  trans¬ 
ition  from  Regime  II  to  Regime  III. 

Transition  Through  Sneed  of  SounAi.  As  observed  in  Pig.  5  and  Table  IV, 
the  predominance  of  drag  in  Regime  I  drives  the  Mach  Number  toward  unity, 
while  the  predominance  of  evaporative  cooling  in  Regime  II  drives  the  Mach 
Number  away  from  unity.  This  permits  a  continuous  transition  through 
Mach  Number  unity  at  the  cross-over  point  where  droplet  drag  and  evapora¬ 
tive  cooling  nullify  each  other  (ref.3);  what  happens  is  analogous  to  the 
passage  through  the  speed  of  sound  in  an  isentropic  Laval  nozzle,  wherein 
the  tendencies  of  the  contracting  section  to  force  the  Mach  Number  toward 
unity  and  of  the  diverging  section  to  force  it  away  from  unity  are  exactly 

nullified  at  the  geometric  throat. 

The  longitudinal  Mach  Number  distributions  are  shown  in  Fig.  6a  for 
a  subsonic  entry.  For  a  low  initial  Mach  Number,  M  rises  from  a  to  b 
(drag  controlling),  then  decreases  (evaporation  controlling).  An  in¬ 
crease  of  initial  Mach  Number  from  a  to  c  increases  from  b  to  d  the  Mach 
Number  at  which  evaporation  begins  to  control  over  drag.  There  is  a 
unique  initial  hhch  Number  (depending  of  course  on  such  parameters  as 
initial  stagnation  temperature,  water-air  ratio,  etc.),  called  the  critical 


l>^oh  Humber,  for  ^loh  th#  cross-over  point  occurs  exactly  et  Mach 
Mumber  unity.  This  condition  yields  curve  e-f,  downstream  of  which  the 
Mach  Number  nay  either  ascend  into  the  supersonic  region  or  descend  into 
the  subsonic  region,  for  either  event  satisfies  the  requirement  of  Table  IV 
that,  in  the  region  where  evaporation  controls,  the  Mach  Number  proceed 
away  from  uriity.  Whether  the  flow  follows  the  subsonic  curve  f-g  or 
the  supersonic  curve  f-h  depends  on  the  back  pressure  of  the  system. 

Just  as  in  the  case  of  the  Imval  nozsle.  Thus,  contrary  to  widely-held 
belief,  it  is  poesiblo  for  a  flow  to  pass  from  subsonic  velocity  to  super¬ 
sonic  velocity  in  a  duot  of  constant  area. 

Ifathematically,  the  transition  through  point  f  may  be  interpreted 
in  terms  of  the  algebraic  relation  for  dM  /M*  represented  by  the 
first  line  of  influence  coefficients  of  Table  III.  For  our  present  purposs 
we  nay  write  this  as 


dM^  Qfe) 

a*  “  i-M* 


(Bq.  8a) 


where 


ir.  which,  for  a  given  case,  G  is  a  function  of  €  alone  inasmuch  as 

all  the  stream  properties,  etc.,  and  their  derivatives,  dM/(ii,dT/dz:, 

etc.... are  functions  of  2  alone. 

From  what  has  been  said  before,  the  sign  of  Q  depends  on  whether 
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drag  or  evaporation  is  controlling.  The  sign  of  cIMYm^  , 
however,  depends  not  only  on  the  sign  of  G,  but  also  upon  whether  the 
flow  is  subsonic  or  supersonic,  as  summarized  in  Table  V. 


IabIo...Y  -  31gn_oiLdM^/dz 


M  <  1 

M  =  1 

M  >  1 

G  >  0  (drag  controls) 

+ 

o« 

- 

G  »  0 

0 

indeterminate 

0 

G  <  0  (evaporation  con¬ 
trols) 

- 

+ 

For  subsonic  flow,  the  table  clearly  shows  a  maximum  in  the  curves 
at  the  point  where  (i.e.  evaporation  and  drag  evenly  balanced),  a 

rising  Mach  Number  in  Regime  I,  and  a  falling  Mach  Number  in  Regime  II. 
What  is  of  additional  Interest  is  that  jg  of  the  indetermin¬ 

ate  form  OfO  when  Q  passes  through  zero  as  M  becomes  unityj  it 
is  this  mathematical  singularity  which  permits  passage  through  Mach  Number 


xinity.  The  slope  of  the  Mach  Number  curve  at  point  f  may  be  found  by 
applying  L'Hospital's  Rule  to  Eq.  8a,  from  which  it  is  found  (ref.  3) 
that 


(Eq.9) 


Thus,  when  M“i  and  ^  “0  ,  the  flow  may  be  continued  only  if 

is  decreasing,  and  in  that  case  there  are  two  solutions  having  slopes 
of  equal  magnitude  but  opposite  sign,  one  proceeding  to  supersonic  speed 
and  the  other  to  subsonic. 


of 


The  erltioel  curve  e-f-g  and  e-f-h  is  a  elngalar  curve  in  the  faml)- 
eolutlou.  of  B,.  8a,  and  the  point  f  le  knovm  ae  a  eaddle  point  singular- 


ity. 

Ileferring  again  to  Fig.  6a,  an  increase  of  initial  Mach  Bi»- 
her  fro.  e  to  j  causes  the  Mach  HuD*er  to  reach  unity  while  drag  is  stiU 
controlllng(l.a-.  with  0  still  positive).  Onder  these  conditions  it  is 
iapossihle  for  the  solution  to  ho  continued  (ref.  3),  and  we  say  that  the 
flow  is  chohed.  In  practice,  the  flow  represented  hy  the  curve  J-k  is 
poeslW.  only  if  the  duct  .r«i  is  lnoreas«i  at  k  or  earlier.  If  the  duct 
contirme.  at  constant  area,  a  steady  flow  with  initial  Mach  taber  3  may 
«t  he  realised,  and  the  highest  initial  Mach  Hurnher  which  is  possible  is 

that  at  e,  hence  the  tern  critical  Mach  Number. 

^oek  Waves .  When  the  flow  becomes  supersonic  on  the  branch  f-h  of 
ng.  6a,  atm  another  feature  in  the  form  of  normal  .shocks  may  add  to  the 
variety  of  poseihle  flow  patterns.  These  take  the  flow  from  supersonic 
ep.«i  to  subsonic  speed,  after  which  the  Mach  Number  decreases  according 
to  Table  V.  Thus,  depending  on  the  back  pressure,  such  Mach  Number  dis¬ 
tributions  as  o-f-n-l-m  or  o-f-n-p-q  nay  bo  observed. 

Entry.  All  the  dlsousslon  relating  to  Figs.  5  and  6«-efor 

to  a  nbsonlc  initial  Mach  Number.  For  a  supersonic  Mach  Number  at  the 
plane  of  injection  the  altuation  is  analogous:  there  wiU  be  the  same  three 
predominating  Begimes  Already  described,  but  the  directions  of  change  will 
in  many  instances  bo  different.  A  doecriptlve  picture  of  the  flow  may, 
howevw.  be  constructed  from  Tables  IV  and  V,  and  Fig.  6b  shows  some  of 
the  Mach  Number  dietributiona  which  may  arise. 
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Tendency  of  Congtantrarea  AerotfegErMsoH 

Perceptive  readers  may  have  already  noticed  in  ^ig.  6a  a  nem'otic 
inclination  of  the  constant-area  Aerothermopressor  ..  ;ard  either  self- 
emasculation  or  suicide. 

Referring  to  Eq.  1,  or  to  the  last  two  lines  of  Influence  coefficients 
of  Table  III,  it  is  seen  that  the  fractional  gain  in  stagnation  pressure 
is  mathematically  represented  by  the  product  of  M  and  an  expression 
representing  the  net  difference  between  the  adverse  effects  of  dreg  and 
wall  friction  and  the  favorable  effects  evaporative  cooling,  an  ex¬ 
pression  which  for  brevity  we  now  call  E  }.  The 

term  E  is  positive  when  evaporative  cooling  is  predominant,  i.e.,  when 
the  stagnation  pressure  is  rising.  Now  E  is  influenced  rather  strong¬ 
ly  by  the  local  toch  Number  because  the  rate  of  evaporation  is  controlled 
by  heat  transfer,  and  the  latter  in  turn  depends  on  the  difference  between 
the  static  temperatures  of  the  gas  and  the  liquid.  Thus,  other  things 
being  equal,  E  is  a  maximum  at  Mach  Number  zero,  and  decreases  with  in¬ 
creasing  Mach  Number  until  it  becomes  negative,  at  which  condition  the 
stagnation  pressure  decreases. 

The  product  M^E  ,  which  is  a  meastire  of  the  fractional  change 
in  stagnation  pressure,  therefore  starts  at  zero  for  zero  Mach  Kumbor,  then 
increases  as  M  increases,  and,  after  passing  through  a  maximum,  decreases 
until  it  finally  becomes  negative  at  a  sufficiently  high  Mach  Number  (occ 

Pig.  7a). 


ft1hsnn^e  Impotence.  In  all  the  auhsonlc  curves  of  Fig.  6a,  the 
predominanoe  of  evaporative  cooling  caueee  the  Mach  Kuniber  to  decrease 
in  Regina  II.  By  tending  to  reetraln  the  gas  temperature  from  decreas¬ 
ing,  this  Mach  Number  reduction  in  turn  is  favorable  for  maintaining 
high  ratoe  of  evaporation,  and  thue  the  tendency  toward  decrease  in 
Ffach  Number  is  further  enhanced.  The  net  result  is  that  of  evaporative 
cooling  is  strongly  ascendant  over  dreg  and  wall  friction.  Though  this 
may  seem  desirable,  it  is  of  little  avail,  since  the  potency  of  the 
process  Is  vitiated  by  the  low  Mach  Number  which  the  high  value  of  E 
produces,  and  the  gains  realisable  from  evaporative  cooling  are  to  a 
degree  forfeited. 

In  brief,  the  process  is  self-defeating  because  at  the  very  moment 
when  evaporative  cooling  stands  at  the  center  of  the  stage  ready  to  produce 
large  gains  in  stagnation  pressure,  it  finds  itself  with  equal  fervor  reduc¬ 
ing  the  Ifach  Number  to  the  point  where  what  should  be  a  shout  is  only  a 

whisper J 

g^ipersonlc  Suicide.  Suppose  now  that  the  flow  has  achieved  the 
supersonic  curve  e-f-h-n  of  Fig.  6a.  Having  arrived  in  Regime  II, 
evaporative  cooling  pushes  the  flow  to  higher  and  higher  Mach  Numbers, 
both  effects  inducing  a  rapid  drop  in  the  gas  temperature.  As  a  con¬ 
sequence,  the  evaporation  rate  is  seriously  reduced.  Ultimately,  the  rate 
of  evaporation  is  reduced  to  such  a  negligible  amount  that  the  stream 
would  scarcely  change  state  any  further  except  for  the  effects  of  wall 
friction  which  are  now  predominant.  To  be  sure,  the  decelerating  in¬ 
fluence  of  wall  friction  permits  a  slight  revival  in  the  rate  of  evap¬ 
oration^  but  this  only  produces  a  battle  of  attrition  in  which  evaporation 


occurs  only  Insofsr  os  the  chauseo  hrousht  shout  by  wsll  friction  per- 

mit  It  to,  and  no  net  gain  ie  possible. 

The  suicldsl  tendency  is  nou  reveslsd:  At  eupereonlc  speeds,  the 

predoBinsnce  of  evsporstl,.  cooling  Inoressss  the  Msoh  hunber,  psvlng 
the  usy  for  wsll  friction  to  bsoow  prsdonlnsnt  end  thus  negstlng  the 

benefits  of  further  evaporation. 

. . .  nf 

in  the  role  of  psyohlstrlst,  do  to  .««r.  s  «r.  heslthy  behsvlor  of 
the  A,rother«,«.sor.  with  th.  sin  of  extrscting  the  n^clnu.  poselUe 
rise  in  stsgnstlon  pressnrel  He  »,  In  foot  do  s  greet  desl,  by  con¬ 
trolling  the  mch  Hunber  distribution  either  through  the  medium  of 

in  oroee-sectlonsl  ores  of  the  evsporstion  section  or  through 
the  Introduction  of  sddltlonsl  wster  st  downstresm  ststlons.  The  former 


possibility  is  discussed  below. 

in  subsonic  flow,  s  decresse  In  ores  1.  st  first  required  to  prevent 
the  Hsoh  Humber  from  felling  too  rspldly.  At  eupereonlc  speeds  the  ares 
must  slso  t.  diminished  In  order  to  prevent  the  Msoh  Humber  from  becom¬ 
ing  too  high,  in  either  osse,  however,  ss  the  process  proceeds,  s  point 
is  reached  at  which  it  is  desirable  to  hsve  a  falling  Meoh  Humber  In 


order  to  mrdntain  the  evaporation  rate,  and  this  ultimately  dictates  an 
increase  in  area  at  some  subsonic  speed.  Thus  the  evaporation  section 
merges  continuously  with  the  diffuser,  and  there  Is  no  real  line  of  do- 
marcation  between  the  two. 

Kbom  this  discussion  emerges  the  significance  of  area  control  as  a 


Mora  Is  said  on  this  Important 


means  for  achieving  ms-ximum  performance. 
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subject  Ic  Sections  4*5  end  4*6. 

Msoh  Nunber  control  by  additional  vator  injeotion  depends  on  the 
fact  that  innedlately  after  injection  the  drag  of  the  additional  vater 
tends  to  predoninate,  and  this  provides  a  means  whereby  the  Mach  Number 
may  be  pushed  towards  xuiity* 


4.5  Appro>rf»nfttfl  Kfathod  for  Selecting  Best  Area  Variation  in. .Evaporation 
Section 

The  discussion  in  the  preceding  section  me^  be  given  an  approximate 
analytical  form  which  provides  a  design  method  for  selecting  a  variation 
in  oross-sectlonal  area  which  is  at  least  in  the  neighborhood  of  that  for 
optimum  performance. 

Detailed  calculations  of  the  type  shown  in  Section  4*7  indicate 
that  the  total  distance  occupied  by  Regime  I  is  only  a  few  inches,  and 
that  only  a  small  fraction  of  the  evaporation  occurs  during  this  period. 

In  Regime  II,  where  most  of  the  evaporation  occurs,  the  droplets  have  a 
speed  very  near  that  of  the  gas  stream  (i.e.,  ^  )»  liquid 

temperature  is  nearly  constant  at  a  value  differing  little  from  the  local  wet 
bulb  temperature,  dy  assuming  that  ^  liquid  temper¬ 

ature  is  constant  and  equal  to  an  average  value  of  the  wet-bulb  tempera¬ 
ture,  Eq.  B-9  expressing  the  fractional  riso  in  stagnation  pressure  per 
unit  of  humidity  increase,  may  be  derived,  as 


shown  in  Appendix  B. 

Design  Procedure  for  Regime  II.  As  explained  in  Appendix  B,  the 
choice  of  the  parameters  and  ifv  permits  the  construc¬ 

tion  of  the  curves  such  as  shown  in  i'ig,  7a,  For  each  value  of  OJ  , 
the  curves  show  clearly  the  insignificant  gains  achievable  at  low  Mach 
Numbers,  and  the  rapidly  deteriorating  performance  at  excessively  high 


Maoh  Numbers. 


To  appreciate  Fig.  7a,  vs  first  ask^  vhat  is  the  criterion  for  best 

performance?  The  answer  is  simply  that  the  net  stagnation  pressure  rise 

of  the  entire  process  be  a  maximum.  Now  this  does  not  necessarily  re- 

rise 

quire  that  the  most  stagnation  prossur^be  realized  for  each  mcremont 
evaporated,  for  a  sacrifice  at  an  early  stage  may  augment  the  potential¬ 
ities  at  a  later  stage.  Consequently,  the  maximization  of  the  overall 
stagnation  pressure  rise  is  no  simple  matter,  transcending  by  far  in  dlf- 
ficulty  the  conventional  variational  problem  of  the  calculus.  Dealing  with 

this  in  its  full  complexity  necessarily  hinges  on  the  numerical  Integra- 
in 

tions  discussedy^^Section  4.7.  With  the  present  simplified  design  procedure, 
there  is  no  better  choice  than  to  make  the  plausible  assumption  that  a 
design  which  maximizes  the  stagnation  pressure  rise  for  each  interval 
of  evaporation  will  in  fact  also  approximate  the  design  which  gives  max¬ 
imum  over-all  nise  in  stagnation  pressure. 

The  foregoing  statement  is  equivalent  to  saying  that  the  dashed  line, 
labelled  "design  cur/e"  in  Fig.  7a,  passing  through  the  maximum  points 
of  the  curves  of  constant  Cd  ,  yields  at  least  a  first  approximation 
to  the  best  Aerothermopressor  design  in  Regime  II.  IVansf erring  each 
point  of  the  "design  curve"  to  Fig.  7b,  we  obtain  the  solid  line  representing 
the  optimum  variation  of  M  with  CJ  .  The  solid  c\irve  of  Fig,  7c  is 
determined  by  plotting  the  value  of  each  point  of  the 

"design  curve"  against  the  corresponding  value  of  OJ  j  this  curve  represents 
the  raaxlrauin  fractional  rise  in  stagnation  pressure  per  unit  of  evaporation. 
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By  Integrating  the  area  undw  the  curve  of  fig*  7o,  the  change  in  stag¬ 
nation  preeeure  aaaoclated  with  a  given  aaount  of  evaporation  nay  then 
be  computed,  ftom  the  equation  of  continuity  it  is  then  possible  to 
oaloulate  the  croee-aectlonal  area  for  eaoh  valve  of  (a)  .  Then,  using 
Eq.  B-B,  the  longitudinal  oocrdlnat©  of  each  point  may  be  detei-niried  by 
numarioal  integration,  than  eatabllehlng  the  shape  of  this  portion  of  the 
Aerothermopreesor  • 

The  method  described  is  found  to  lead  to  the  result  that  the  Mach 
Number  must  be  decreased  as  the  evaporation  proceeds.  This  generally 
dictates  a  reduction  in  cross-sectional  area  in  the  early  part  of  Regime  II, 

followed  fay  an  increase  in  the  later  part. 

leffeats  of  Initial  Paraaetura  on  Bast  Design.^  Figs.  7b  and  7c  illus¬ 
trate  with  dashed  lines  the  comparative  positions  of  the  '•design  curves" 
which  result  flrom  fig.  7a  for  different  values  of  the  initial  parameters 
mere  favorable  to  effective  operation.  Such  more  effective  operation 
OOMS  ft*om  higher  values  of  \  and  t  and  ftrom  lower  values  of  d, 

and  f /D  . 

Under  these  olrcuBUitanees  the  curves  of  constant  <*)  in  Fig.  7a 
are  dlaplaoed  upwards  and  to  the  right.  This  in  turn  moves  the  corres¬ 
ponding  "design  curves"  of  Figs.  7b  and  7o  upwards*  Thus,  with  more  favor- 
abla  Initial  parameters,  the  level  of  operating  Mach  Numbers  is  higher, 
and  greater  inoreases  in  stagnation  pressure  per  xmit  of  evaporation 
may  be  aohieved. 

Bvaluatlon  of  aimnlified  Design  PTMCdur^^  iha  principal  defect  of 
the  method  outlined  above  is  that  it  le  incapable  of  revealing  the  correct 
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Nosselt  Humber  to  be  used  at  each  stage.  To  determine  Na  requires 
an  evaluation  of  the  velocity  differenoe  between  the  phases,  and  this 
can  be  found  only  from  a  oonsideration  of  the  entire  past  history  of 
the  procesa#  The  anaOLyais  gives  recognition  only  to  the 

present  state  of  affairs,  and  does  not  look  to  the  past  or  to  the  future. 
During  the  course  of  the  process,  the  value  of  Na  change,  oontlnuously, 
and  a  series  of  charts  like  that  of  fig.  7a,  each  with  a  different  value 
of  Na  ,  would  be  really  necessary  for  establishing  the  beet  design 
(assuming  that  1*1  were  known  at  each  point).  This  limitation  is 
more  Important  than  it  might  at  first  seam  to  be,  because  the  comparatively 
small  difference  between  the  total  beneficial  effects  of  evaporative  cool¬ 
ing  and  the  total  detrimental  actions  of  drag  makes  a  change  in  Na  from, 
say  2  to  4,  most  significant  in  reorientating  the  optimum  variation  of 


M  with  . 


What  all  this  means  is  that  the 


ol?/?dw 


method  cannot  give 


the  Aerothermopressor  the  advantages  which  might  accrue  from  control 
of  the  Susselt  Humber,  ito  emample,  by  choosing  what  might  appear  to  bo 


a  comparatively  unfavorable  Mach  Humber  distribution  at  one  stage,  the 
ensuing  acceleration  or  deceleration  of  the  stream  might  be  used  to  estab- 
llsh  large  velocity  differences  and  correspondingly  large  Husselt  Numbers, 
with  the  oondequenoe  that  Increased  gains  in  stagnation  pressure  realised 


at  a  later  stage  might  more  than  compensate  for  the  mediocre  performunoe 


accepted  at  the  earlier  stage. 

Thus  the  analysis  advises  the  Aerothermopressor  to 
»ll,e  for  the  present  only,  with  no  thought  for  the  future",  and  conse 
quently  cannot  by  itself  give  the  final  answer  to  the  question  of  how 


the  Aerothennopreesor  should  be  designed.  On  the  other  hand,  it  cer¬ 
tainly  represents  a  starting  point  for  a  more  adequate  approach  such 
as  that  outlined  in  Section  4.7,  and  it  is  hoped  that  ultimately  a 
study  of  specific  results  generated  by  the  rather  involved  procedure 
of  Section  4.7  will  in  turn  suggest  some  simple  rules  of  thumb  with 

which  the  results  of  the  comparatively  simple  ^5 

my  be  modified  to  be  snore  closely  in  accord  with  the  best  design. 

A.6  Some  Design  Coneld«-JtloDfl 

Msarmingly  slaple  as  it  is  in  structural  arrangement,  the  Aero- 
thermopresaor  is  host  to  a  bewildering  variety  of  physical  events. 

These  phenomena,  constituting  fundamentally  mass  conservation,  momentum 
effects,  and  energy  effects  for  both  the  droplet  cloud  and  the  gas  stream 
are  strongly  inter-related.  Consequently  the  designer,  in  attempting 
to  give  due  weight  to  each  factor  influencing  performance,  is  faced 
in  the  Aerothermopressor  with  a  large  number  of  conflicting  tendencies. 
The  best  design  is,  as  usual,  a  oompromlse  among  these.  What  is  baf¬ 
fling  is  that  at  present  there  seems  to  be  no  s imple  way  for  evaluating 
the  relative  ln?)ortance  of  each  effect  and  thus  of  arriving  at  clear- 
cut  deiign  rules.  No  doubt  a  combination  of  sufflolent  experimental 
data  with  worked-out  numerical  cases  such  as  those  of  Section  4.7  will 
someday  help  in  this  regard. 

Having  Bads  thsae  ratber  discouraging  remarks,  us  may  Immedlatsly 
sdd  that  mich  is  indasd  known  In  a  qualitative  way  about  tbs  Influsboss 
of  tbs  Bost  Important  design  variables.  This  knowledge  Is  Important 
not  only  beoauae  It  reveale  pitfalls,  but  also  becauae  It  uaually 


suggestK  how  experimental  information  can  be  used  in  improving  a  given 
design.  Therefore  we  shall  discuss  below  the  major  elements  of  the 
Aerothermopressor,  state  the  decisions  facing  the  designer  in  each, 
and  present  the  major  considerations  by  which  the  design  choices  may 
be  made. 

Design.  Specif icatloni?.  In  a  specific  situation,  we  are  confronted 
by  design  values  of  air  flow  (  ),  initial  stagnation  temperature  (  ), 

and  initial  stagnation  pressure  (  ),  although  the  magnitudes  of  the 

two  latter  quantities  are  in  some  degree  dependent  on  the  performance  of 
the  Aerothermopressor  itself.  In  addition,  there  may  bo  a  range  of  these 
three  parameters  specified  in  the  area  of  off-design  operation. 

Control  of  Aerothermopressor.  Once  a  particular  unit  has  been  built, 
at  least  three  procedures  are  possible  for  adjusting  to  off-doslgn  con¬ 
ditions  or  for  compensating  for  imperfections  in  designs 

(l)  By  far  the  easiest  scheme  is  to  control  the  rate  of  water  in¬ 
jection  at  inlet,  thereby  governing  the  water-air  ratio.  Presumably  it 
would  be  routine  procedure  for  the  operator  to  adjust  the  water  flow  to 
the  point  of  maximum  stagnation-pressure  rise.  At  best,  however,  this 
method  can  take  care  of  off-design  conditions  or  lapses  in  design  only 
approximately, 

(ii)  "TrLmming"  of  the  lengthwise  fech  Number  variation  may  be 
accomplished  by  means  of  variable-area  devices  such  as  a  central  plug 
or  a  series  of  concentric  sleeves  located  along  the  axis  of  the  Aerotherm- 
opressor.  This  would  considerably  detract  from  the  mechanical  simplicity 
which  the  Aerothermopressor  otherwise  possesses. 
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(lu)  Trimins-  of  l.n8tl»U.  M«h  tot®  vTloticn  My  bo 

th.  injootlon  of  .orilluy  wator  at  oovorol  station. 

acng  th.  ^.poratlon  sootlon.  Even  th^gh  the  offoctlvoneos  of  tho 
,„*aiu7  «!«•  BiV  b.  sUghtUr  l-palred  throngh  a  loss  In  rosiaenoo 
tl»..  th.  .ff.oti«n...  of  th.  Bln  l.)«t.a  vator  Bight  b,  angmonted 


Buff identic  to  effect  a  net  gain. 

ftlirntilinlinn - fo"n°n«'t»a  IW  porpo...  of  isolating  ana  a.al- 

iBg  ,.p«.at.l,  with  th.  prlnolpl.  phono.™.,  »  1.  oonvonlont  to  think 
of  th.  i™otl»r«.pr»..r  a.  c^prlslng  the  foUovlng  oloBonts:  (1)  th. 

Uir  »™i.  nhloh  «.ort«.t..  th.  hot  air  ffo.  th.  lov.-.P«a  upstr«. 
™,tl<».  into  th.  hlgh-.p..a  sartlon,  (11)  th.  «at.r  Injection  systoBi 
(ill)  th.  -aocelwatlon"  ...tlon,  la  vhloh  occur  .toBl«itlon  and  th. 

■ajor  portion  of  th.  ar.pl.t  .cc.l«.tlon;  (Iv)  th.  scaporatlon  ssctlon. 
la  uhlch  «.pcr.tlon  1.  g™««Uy  oontroUlng  and  vher.  nnst  of  th. 
„.poratlon  «.cur.,  «.d  (v)  th.  dlffusar,  vhlch  ultlBately  d.c.l.rat.s 

th.  etr...  to  low  .p«d.  «.d  tomlnat..  th.  procs. 

■  Bach  of  th...  ccBpononts  will  now  be  consld.r.d  soparatsly. 

.  4irjadSa  TWO  luostlon.  arise  her.:  (1)  how  high  should  b.  th. 
mch  »UBh«.  M,  at  th.  nossl.  .nt  (which  1.  also  th.  plan,  of  wate 
injection),  ted  (11)  vhat  shape  should  th.  nossl.  have? 

[hr  th.  apwlfled  values  oC  W  ,  ,  «”d  Xi  '  diamster 

at  th.  nossl.  .sit  d.tw.ln.6  th.  Maoh  llUBb«  1^.  If  I*!  io  oxo.sslv.ly 

large,  the  stagnation-pressure  loss  associated  with  th.  droplet  acceler¬ 
ation  Bay  b.  prohlbltl,.,  ho.««.  If  It  1.  too  small,  th.  droplet  dla- 
«.t«.  will  b.  cOBp.ratlv.ly  le-ge,  c«d  the  required  length  of  the  evap¬ 
oration  section  wiu  be  escesslv.,  leading  to  Intclerabl.  friction^ 


losses.  Ths  best  value  of  M|  in  each  case  Dusty  therefore»  be  deter¬ 
mined  either  by  experiment  or  by  caleulations  of  the  type  illustrated  in 
Section  4<7* 

The  influence  of  the  Mach  Number  level  on  the  length  required  for 
evaporation  may  be  seen  fr^  the  follovinfj  approximate  considerations  i 
Fbr  air  atomisation,  the  droplet  diameter  <i^varies  inversely  with  the 
air  velocity  V  (ref.  4).  Since  the  distance  travellod  ly  each  droplet 
during  its  acceleration  period  is  a  very  small  ft'sction  of  the  total  dis¬ 
tance  L,  required  for  evaporation,  L  B«y  b®  approximated  by  suppeslng 
that  the  droplet  moves  with  the  same  velocity  as  the  air.  Plrora  consid¬ 
erations  of  heat  transfer, 


For  very  low  Reynolds  Numbers,  mA  •  Substituting  this  value 

of  ,  and  integrating  ths  equation,  we  find  that  ths  time  for  oosi- 

plete  Bvspotratioo  follows  ths  law  d  .  In  addition,  ths  length  L 

is  proportional  to  the  product  of  V  and  t  •  Therefore, 


Lw-  Vt  d,  l/v 


and  the  required  length. is  seen  to  be  Inversely  proportional  to  the 
velocity.  For  a  given  length,  a  greater  fraction  of  total  evaporation 
occurs  at  high  air  speeds  than  at  low,  even  though  high  air  speed  means 
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a  ooaparatlvely  small  residence  time. 

Iteg^dlng  th.  .tap.  of  tta  .ir  noirt.,  th.  tal.  con.ldor.tiota 

fcor.  «•«  that  (1)  it  .houU  ta  tanglttaltaUr  «•  rt»rt  ..  potalU. 

.0  ttat  -Mlv.  p«t  of  tta  tator  lojootion  .jrrtta  «)r  ta.lly  to 
pl.o«l  in  .  lov-Toloolty  tagloo,  «d  W  «  taoold  tavoroo  prosour, 

gr«U«.t.  00  th.  taU.  ItagO  onoo^  to  Induo.  bo»otary-lv«  .oporotioo. 

STOtan.  Tta  r«pil«»nt.  of  th.  »tar  loj.otlon 
.^tao  ta.  ttat  (1)  th.  droplrt  oloud  produota  rtiooM  tao.  th.  taoHoot 


potalhl.  vol—tarfta.  tata  drop  .1..,  (U)  th.  drop  .1..  .p.otr-  ho 
„ch  ttat  th.  tau  fh«otlon  of  th.  drop.  Iflog  .tor.  th.  Tola»-tarfao. 
MU  .11.  1.  tat  too  ltag«,  (Ul)  th.  droplrt  cloud  .hould  ta  unlforoly 
dlitathutad  ovta  th.  ta...  taoUon.  tad  (iv)  th.  pta..ltl.  .tagootloo- 
priamr.  lo«.  ..MoUtod  vlth  th.  urodjoulo  drag  of  tta  injtatloo  ay.- 

tem  should  be  a  minimum. 

Without  pr.oi«dlng  th.  po..ltUlt3r  of  other  tarugueot.  vhich  m, 
proT.  ur.  f.Yor.bl.,  urt  of  th.  .ffort  thu.  fta  ha.  taao  glvu  to  a 
"porouplo."  typ.  of  Injactor,  ooiaprl.ii«  ta.ll  tuta.  lo  parallel,  with 
tta  water  InJ.otta  at  low  .port  into  th.  high-.p.*i  air  appr«doat.ly 
at  the  .nt  of  the  air  no..!..  UtU.  ure  ou  ta  aaid  now  about  the 
injoutlon  .yatta  .uept  that  th.  laok  of  pr.ol«  Infortatlon  oonoemlng 
th.  drop  .!»  produeod  1.  ptahap.  th.  wuktat  llnh  oonneotlng  theory  and 


experiment . 

Th.  beat  watar-alr  ratio  roaolt.  Iron  a  wingwoniM  between  two  oppoa- 
ing  affect..  That  la,  a.  tta  water  flow  1.  lnor««d,  th.  aoceltaation 
arag  ia  lncr....d  (unfavorabl.) ,  but  tta  Mount  of  ta.peratlon  ooourrlng 


0 


in  a  given  distance  is  also  increased  (favorable).  Fortunately  this 
is  a  matter  which  in  each  case  can  be  very  easily  settled  on  an  oper¬ 
ating  Aerotherraoprsssor. 

oration  Section.  Within  the  very  few  inches  of  length  con¬ 
stituting  the  acceleration  section,  the  droplet  cloud  is  formed  and 
accelerated  approximately  to  the  gas  speed,  and  of  the  order  of  1/4 
the  total  evaporation  occurs.  Because  of  the  exceedingly  short  length 
(ref.  5)  and  the  extra  complexity  of  the  processes  in  this  region,  it 
may  prove  unfeasible  to  attempt  any  Mach  Number  control  by  means  of 
area  variations  in  this  zone. 

A  considerable  loss  of  stagnation  pressure  occurs  in  the  accel¬ 
eration  section.  The  loss  of  stagnation  pressure  associated  only  with 
the  acceleration  of  the  droplet  masses  (disregarding  heat  transfer  and 
evaporation)  is  illustrated  by  the  following  values  (ref.  5)  for  a 
water-air  ratio  of  0.2  (the  loss  is  approximately  linear  with  water- 

air  ratio): 


M,  Poz  IPo, 

0.3  0.978  0.988 
0,4  0.957  0.977 
0.5  0.927  0.961 
0.615  (choked)  O.846  0.933 


These  losses  are  of  the  same  order  as  the  net  gain  in  the  Aerothermo- 
pressor,  thus  emphasizing  strongly  the  necessity  of  careful  design  at 
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every  point  In  order  to  achieve  acceptable  performance. 

For  the  limiting  case  of  vary  low  Mach  Number,  a  formula  for  the 
irrecoverable  loss  of  stagnation  pressure  due  to  droplet  acceleration 
may  be  derived.  From  momentum  considerations,  the  pressure  change 
'Pb-P*  )  due  to  the  acceleration  of  the  drops  from  zero  speed  to 
the  gas  speed  V  without  area  change  is  given  by 


■p  -  p  ss  M  yoK*  (Eq.  9) 


and  is  equal  to  the  loss  in  gas  stagnation  pressure  during  the  accel¬ 
eration,  If  the  mixture  is  then  reversibly  decelerated  to  zero  speed, 
the  pressure  attained  is  calculated  by  considering  the  two-phase 

mixture  to  be  homogeneous  with  the  appropriate  density; 


(Eq.  10) 


Accordingly  the  net  lose  In  stagnation  pressui'e  is  given  by 
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tox  2  ~  K  2. 


(Eq.  11) 


That  is,  the  minimum  conceivable  net  loss  in  stagnation  pressure  is  one 
dynainit!  hear!  nultiplled  by  the  wator-aih'  I'atio,  and  3  s  exactly  ono-hslr 
tlie  loss  in  mi)  stagnation  pressiu'e  during  ilio  accelsi’ation,. 


The  best  water  injection  velocity  is  governed  by  the  opposing 
considerations  that  (i)  a  low  watei-  injection  velocity  results  in 
small  drops  but  large  acceleration  loss,  whereas  (ii)  a  high  water  in¬ 
jection  velocity  gives  large  drops  but  small  acceleration  loss. 

The  required  number  of  injection  points  depends  on  the  angle  of 
spread  of  the  droplet  spray,  JVora  the  data  obtained  thus  far,  it  appears 
that  the  best  number  probably  lies  within  the  range  of  1/2  to  2  injection 
points  per  square  inch  of  cross-seotion. 

Rvanoration  Section.  The  primary  considerations  for  the  design  of 
the  evaporation  section  have  already  been  discussed  in  Sections  4»4  and  4*5. 

Diffuser .  No  precise  line  of  demarcation  in  terms  of  the  Aerothermo- 
pressor  process  can  be  made  between  evaporation  section  and  diffuser,  in— 
asnruch  as  significant  evaporation  and  droplet  deceleration  effects  occur 
within  the  diffuser,  A  substantial  portion  of  jhe  evaporation  may  occur 

B 

in  that  portion  of  the  evaporation  section  which  has /geometrically  diverging 
cross-section.  Moreover,  the  gas  decelerates  and  the  pressure  rises  in 
virtually  the  entire  length  of  the  evaporation  section. 

Granting  these  facts,  a  stage  is  ultimately  reached  in  the  Aerothermo- 
pressor  process  where  the  main  objective  is  to  decelerate  the  stream  and 
recover  pressure  efficiently.  This  part  of  the  apparatus  is  arbitrarily 
called  the  diffuser.  The  best  conventional  diffuser  is  a  cone  of  approx¬ 
imately  6“  included  angle.  A  consideration  of  the  concurrent  phenom¬ 
ena  of  evaporation  and  droplet  deceleration  occurring  in  the  Aerothermo- 
pressor  diffuser,  however,  can  be  interpreted  to  suggest  a  somewhat  longer 
diffuser  in  order  better  to  exploit  these  additional  phenomena.  The  dif¬ 
fuser  of  smaller  angle  makes  the  revived  evaporation  in  the  diffuser  more 
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effective  because  it  causes  the  evaporation  to  occur  at  relatively 
higher  Mach  Numbers;  and  it  adds  to  the  efficiency  of  droplet  decel¬ 
eration  by  causing  the  droplets  to  remain  more  nearly  at  the  gas 
speed  (ref.  6). 

Ontumim  Mach  Number  Distribution  for  ibtire  Aerpinermopressop. 

Guided  by  m  combination  of  the  qualitative  considerations  presented  in 
this  Section,  together  with  approximate  design  method  for  Regime  II 
given  in  Section  4.5  and  the  results  obtained  by  the  more  complete  theor¬ 
etical  method  discussed  later  in  Section  4-7 i  the  nature  of  the  optimum 
Jtach  Number  distribution  for  the  entire  Aerothermopressor  may  be  ration¬ 
ally  anticipated.  The  dlscuseion  is  illustratei  schematically  in  Fig.  7d. 

Jbon  a  oonsideratlon  of  the  importance  of  maintaining  a  large  ten^j- 
erature  difference  for  heat  transfer,  we  may  confidently  expect  that  the 
position  of  Mach  Number  (point  b)  should  occur  near  the  begin¬ 

ning  of  the  process  where  the  humidity  is  low  and  the  stagnation  temp¬ 
erature  high. 

Now,  referring  to  the  approximate  method  for  Regime  II,  large  val¬ 
ues  of  the  relative  Reynolds  Number  can  be  established  only  in  the  part 
of  the  process  where  the  llach  Number  is  high  and  vdiere  large  accelera¬ 
tions  and  decelerations  are  possible.  Thus,  in  the  neighborhood  immed¬ 
iately  downstream  of  point  b,  the  curve  of  M  vs.  U)  will  be  that  of 
Fig.  7b  Tor  a  Nuscelt  Number  of  about  4.  Near  the  end  of  the  evapora¬ 
tion  region  the  Mach  Number  is  necessarily  amaller,  the  speed  of  sound 
is  lower,  and  large  accelerations  or  decelerations  are  no  longer  possible# 
Together  with  the  fact  that  the  droplet  diameter  is  also  smaller,  all 
this  implies  that  the  relative  Reynolds  Number  is  quite  small.  Therefore, 


near  the  end  of  the  evaporation  section  (point  c),  the  curve  of  M  vs.  O) 
should  be  like  that  of  Fig.  7b  for  a  Nusselt  Number  of  about  2, 

Assuming  that  point  c  represents  the  location  vhere  the  rate  of 
sttignation  pressure  rise  has  virtually  vanished,  the  curve  c-d  will 
represent  the  diffuser.  The  latter,  as  e^qplained  previously,  will  most 
probably  be  a  cone  of  half-angle  somewhat  less  than  3®. 

The  best  inlet  Mach  Number  (point  a)  cannot  be  foreseen  easily,  as 
it  depends  on  the  best  couqjromise  between  small  droplet  size  and  small 
acceleration  drag.  In  some  cases  the  inlet  Mach  Number  might  be  the 
highest  of  the  entire  process,  and  then  M  would  decrease  continuously. 

The  value  of  U)  at  point  b  is  determined  by  the  advantage  of 
making  the  Mach  Number  rise  from  a  to  b  occur  as  quickly  as  possible, 
with  the  restraint,  however,  that  too  rapid  an  acceleration  will  maintain 
such  large  differences  in  velocity  that  large  losses  in  will  be 

produced.  It  seems  extremely  important  to  choose  the  location  and  ^fach 
Number  of  point  b  wisely,  for  the  nature  of  the  entire  process  calls  for 
the  maximum  exploitation  of  the  benefits  of  cooling^the  very  region  where 
evaporation  is  most  predominant,  where  the  temperature  differential  is 
still  large,  and  where  there  exists  the  possibility  of  operating  locally 
at  large  Mach  Numbers.  Immediately  downstream  of  point  b  is  where  the 
maximuni  gains  are  to  be  won,  and  if  they  are  not  won  there,  they  are  for¬ 
ever  forfeited. 

In  summary  then,  the  best  Mach  Number  distribution  for  subsonic  entry 
appears  to  call  for  a  soaring  rise  of  Mach  Number  in  the  acceleration 
zone,  followed  by  a  gradual  decrease  in  Mach  Number  in  Regime  II  (involv¬ 
ing  a  gradual  transition  from  the  curve  MU^4-  to  that  for  l*lu  *  2.  ), 


and  concluded  by  e  fairly  rapid  decrease  of  t-tech  Mumbor  in  the  dlffuaer 


Because  of  its  complexity,  its  many  diverse  design  considerntions, 


and  the  fact  that  only  by  great  care  can  a  poor  design  be  avoided,  a 


successful  future  for  the  Aerothermopressor  may  depend  as  irtuch  on  the 
develOFoent  of  a  rational  theory  as  on  experijnental  development.  The 
latter  goes  slowly,  and  nay  be  of  greatest  significance  in  helping  to 
font  a  theoretical  design  method  which  adequately  represents  Aerotherno- 
preBS(^  perforiKmos.  Then  we  may  look  to  the  theory  for  such  matters 
as  determining  the  optimum  cross-sectional  area  variation,  extrapolat¬ 
ing  performance  to  unknown  ranges  of  supply  pressures  and  teiDperaUu-es, 


and  predicting  off-design  performance. 


The  analysis  considered  is  purely 


ona-dlnenslonal*  in  aiiy  cross-sectional  piano  (i)  the  gas  velocity, 


temperature,  and  pressure  are  uniform,  and  (ii)  the  droplot  velocity, 
temperature,  slwj^and  spatial  density  are  uniform.  In  addition  it  is 
assumed  that  the  droplets  are  atomized  immediately  upon  injection,  a:.d 


that  no  subsequent  atomization  or  agglomeration  occurs. 

nf  the  Analysis.  At  the  entry  plane  of  the  Aerotherwo- 

pressor  (section  1),  the  initial  conditions  may  be'choracterized  by 
seven  quantities,  of  which  the  first  six  are  intensive  and  the  last  ex- 


tensive i 


ft. 

(ii)  X. 

(iii)  M, 

(iv)  12^ 


,  inlet  gas  stagnation  pressure 
,  inlet  gas  stagnation  temperature 
,  inlet  Mach  Kumbar 
,  inlet  water-air  ratio 
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(v) 

,.  inlet  liquid  temperature 

(vi) 

,  inlet  liquid  velocity 

(vil) 

,  inlet  cross-sectional  area 

Note  that  P„  >  X.  .  M,  ,  and  Aj  together  determine  the  air 
flow  rate  'K  ,  no  that  could  be  used  in  place  of  Aj  ..  Further¬ 

more,  for  a  given  style  of  water  injector,  we  may  anticipate  that  the 
listed  variables  v/ill  determine  the  initial  droplet  diameter  . 

The  duct  itself  may  have  a  shape  which  might  be  characterized  by 
relating  the  area  ratio  to  the  distance  2  from  the  inlet  plane;  For 
example,  we  might  write 


(Eq.  12) 


where  OC^  ,  OC^  ,  etc.,  are  dlmonsionleas  numbers  defining  the  duct 
shape. 

At  any  section  2  dounstrenm  of  the  inlet  j)lane,  having  a  local 
cross-sect: onal  area  A  ,  the  state  of  the  stream  may  be  characterized 
b"  six  intensive  quantities; 


(i) 

f 

,  gas  pressure 

(ii) 

T 

,  gas  toMperaturo 

(iii) 

M 

,  Mach  Number 

(iv) 

Ct) 

,  spoeiflc  huMidity 

(v) 

Tt 

,  liquid  temporature 

(vi) 

,  licpjid  velocity 

Qo¥Tnlng  Bnatlmn.  Tfatts^  for  opooiflod  valms  of  tho  ••von 
inltlid  varlabloo  and  a  8|)aolflad  duot  ahapa^  alz  pfaTaioallf  Indapanlant 
af}oatlmt  ara  naoaaaarjr  for  aatabliriilng  tha  six  rattulrad  propertlaa 
at  aaoh  atatlon  Z  •  niaea  six  aquatlono  ara  fVindaMntally  (1)  tha 
aqoationa  of  oontlnuityi  iMBBntiaif  and  ettargy  ter  tha  gas  atraaaf  and 
(11}  tha  aquatlona.of  naaa  tranafar,  haat  tranafarf  and  drag  for  the 
droplet  olottdi 

Tha  ^upaa  goaarnlBg  a<|uatlona  for  the  gaa  atraaa  ara  darlvad  In 
Appaadlz  Ai'  and  thalr  alaailtanaoua  aolutlon  baa  already  bean  atnnarlead 
hgr  aaana  of  tha  Inflnanoa  Ooafflolant  aquatlooa  of  Table  111.  In  theae 
aquatlona  the  tama  and  are  dapandant  on  anas  tranafar  (evap¬ 
oration)  lyon  tha  droplatai  tha  tern  la  dapandant  vq>on  ohangaa 

in  taaparatnra  and  hnnldlty  of  tha  gaa  atraan)  tha  tarn  In  gIA|  la 
dapandant  on  tha  Maa  tranafar  and  haat  tranafar  for  the  droplatai  and 
the  tam.in  ci\^  la  dapandant  on  tha  drag  foroaa  batvaan  tha  droplet 
and  gaa.  Tharaforay  a  oomplata  analytloal  fomulatlon  of  tha  problaai 

.  9 

raqulraa  that  the  aqnatlona  of  Table  III  be  augnantad  by  the  three  equa- 
tlona  of  nasa  tranafar,  haat  tranafar,  and  drag  for  tha  droplata.  Tha 
derivation  of  the  latter  equations  ia  praaentad  in  i^ndix  G. 

Wimerlcal  flolntioa  of  I  The  syatan  of  ordinary,  alniltan- 

aphs,  non-llnaar,  differential  aquatlona  referred  to  above  aaana  to  be 
vail  beyond  tha  poaaibllity  of  analytloal  integration,  and  nuaarloal 
aolutlon  la  therefore  naoaaaary.  After  a  brief  exparlanoa  with  working 
aolntlona  by  hand  oaiqmtatlon,  it  baeana  evident  that  real  prograsa 
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could  be  -id.  only  by  .  -oh  -re  r.pld  cloulatlon  procod.«.  Iho 
problsm-was  accordingly  prcgrai«.cd  for  the  Wilrlvind  I  computer,  a 
hlgh-apeed,  electronic,  digital  machine  at  the  Haa.achua.tte  Inatltute  . 

of  Technolo^a 

Although  extended  dlacuaalcn  of  the  ««..  of  aclutlcn  la  luappr^rlate 
here,  It  may  be  mentioned  In  paaalng  that  a»e  of  th.  difficulties  en¬ 
countered  d«-lng  the  development  of  thl.  numerical  treatment  involved 
(0  adeciuately  large  hlgh-ep.«i  memory,  (11)  tolerably  small  truncation 
error  without  exooaslvely  large  computational  time,  (ill)  tendency  of 
the  nuamrlcal  solution  toward  Serl-a  InaUbUlty,  a».  (Iv)  .stabllah- 
ment  Of  the  alngular  .elution  paaalng  oontl-ou.ly  through  Kaoh  hu.^ 

unity.  More  details  are  given  by  OavrU  (ref.  8). 

At  th.  ire»,nt  time  the  f<«,rth-ord.r  Runge-Kutta  method  of  numerl 

leal  integration  1.  uaad.  l»t  the  general  scheme  of  calculation  may  be 
more  easily  Illustrated  by  considering  the  flrat-order  (Euler)  method  of 
numerical  Integration.  In  the  latter  th.  basic  Idea  Is  to  calculate  the 
increment  at  a  point  using  only  the  first  to-  of  the  Taylor  Series  expansion 
about  the  point,  for  .:»mple,  AM**  {dU/dijlX.  ^PPoee  that  the 
initial  values  of  the  problem  are  given,  and  that  the  fraction  evaporated, 
ym  (o/jI2j  ,  which  is  taken  as  running  variable,  la  divided  Into 

many  small  intervals  .  The  finite  difference  equation,  are  then 

solved  for  «ich  interval,  and  the  propertlea  computed  for  the  «-d  of  the 
interval  ere  then  available  for  use  in  the  next  interval. 

TO  Illustrate.  Imagine  that  the  caloulatlcna  have  been  carried  to  a 

certain  value  of  X  .  where  the  values  of  W  ,  <^,4-8  ’f’ 

r  1/  M  z  ,  &n6  4  are  all  knovm.  Than,  for  an  increment 

in  humidity  AX  !  the  values  of  Aoi  and  AUJ  are  calculated  from  Fq.C-1; 


/i1|  l8  oonqnited  from  Eq,  C-5b  vith  Bqs.  C-3»  C-3a»  C-4»  C-6,  and  C-7| 

M  la  found  from  Eq.  C-5a  with  Eqa.  C-6  and  C-7;  le  obtained 

from  Eq.  C-8  with  Eqa.  C-9;  and,  finally,  tba  distance  t  in  which 
the  foregoing  changes  in  droplet  properties  occur,  is  determined  by  Eq.6-10 
Assuming  that  the  duct  shape  is  Imown  may  be  found  for  the  corres¬ 
ponding  ah  ,  All  the  increments  In  the  top  row  of  Table  III  sre  now 
known  for  the  step  under  consideration,  and  the  corresponding  changes 
in  gas  [UTopertiaa,  ,  Ap  ,  AV ,  etc.,  may  be  reckoned  from  the 
equationsof  Table  III.  Having  carried  out  those  calculations,  the  drop¬ 
let  and  gas  properties  are  con^letely  determined  at  the  end  of  the  hum¬ 
idity  Interval,  and  all  is  in  readiness  for  repeating  the  procedure  within 
the  next  humidity  interval. 

It  should  be  noted  here  that  the  solution  is  propagated  by  the  fact 
that  at  any  section  t  ,  where  all  the  stream  properties  are  known,  it 
is  posslblo  to  calculate  the  first  derivative  of  each  property  by  means 
of  the  oquatioitfof  Appendix  C  and  Table  III,  The  Increment  in  each  prop¬ 
erty  can  therefore  be  determined  by  any  of  the  conventional  numerical 
procedures.  One  of  the  crucial  aspects  of  such  a  calculation  is  that  each 
sucoessive  value  depends  on  previous  vedues. 

Present  Status  of  Numerical  Calculations.  At  the  present  writing 
a  considerable  number  of  integrations  has  been  carried  out  for  the  con¬ 
stant-area  Aerothermopressor  and  a  smaller  number  for  the  variable-area 

case.  Currently,  work  is  proceeding  on  a  promising  soheme  for  eystem- 
the 

atloally  optimising^parameters  and  the  area  variation  controlling  per¬ 
formance,  but  this  has  not  yet  progressed  to  the  point  where  results 
may  be  reported.  Many  of  the  results  thus  far  obtained  have  been  presented 


by  Gavril  (ref.  8).  A  few  are  presented  here  for  illustration,  and 
in  Section  5  we  shall  present  additional  computed  results  for  comparison 
with  the  experimental  results* 

Constant-Area.  Subsonic  Aerothermporessor  Process* .  fig.  8  shows 
the  results  of  numerical  calculations  for  the  case  of  a  constant-area 
unit  with  an  initial  Mach  Number  of  0.7,  a  vater-air  ratio  cf  0.20,  initial 
stagnation  conditions  of  1600  R  and  14.7  psia,  a  friction  parameter  f/ 
of  0.004  ft  (corresponding  to  0.004-  with  a  duct  of  12*  dia¬ 
meter),  and  an  Initial .  droplet  diameter  of  13.8  microns  (as  conpited 
from  the  Nukiyama  and  Tanasawa  formula,  ref,  4),  All  the  varied  phenom¬ 
ena  previously  described  as  occun*lng  In  Regimes  I,  II  and  III  are  seen 
here  quantitatively.  Noteworthy  points  Include:  (i)  Regime  I  is  of  the 
order  of  only  0.1  ft.  in  length,  although  the  actual  delay  in  atomiza¬ 
tion  will  increase  this  somewhat  j  (il)  the  liquid  overtalces  the  gas  in 
about  0.4  ft}  (lii)  about  10^  of  the  water  is  evaporated  in  Regime  I  and 
about  30%  at  the  end  of  the  droplet  acceleration  period}  (iv)  the  1%  loss 
in  stagnation  pressure  incurred  in  Regime  I  is  not  recovered  until  .bout 
50^  of  the  water  Injected  is  evaporated,  thus  placing  the  burden  of  win¬ 
ning  a  net  gain  on  the  remaining  part  of  the  process}  (v)  the  net  rise 
in  stagnation  pressure  is  about  (vi)  the  required  length  of  evapora¬ 
tion  section  is  of  the  order  of  4  feotj  and  the  droplet  residence  time 
is  of  the  order  of  0.01  seconds, 

Gonf.-tnnt-Arna.  Supersonic  Aerothormonressor  Process.  Fig.  9  is  for 
the?  sane  case  as  Fig.  8  except  that  the  initial  Mach  Number  is  1.5  and 
tho  corresponding  initial  droplet  diameter  is  7.34  microns  (based  on  the 


tha  risky  sxtrapolation  of  the  Nukijuks  end  Tansssva  drop-slse  forauls  to 
aiq>ersonlo  speeds)*  Here  agsln  the  dlstinot  phenonena  of  Regies  If  II 
and  in  are  olearly  evldentf  althou^  they  are  generally  of  opposite  ehar- 
aoter  to  the  subsonle  case*  The  lengths  of  Regies  I  and  of  the  evapora¬ 
tion  tone  are  roughly  the  saee  as  those  in  Fig*  8.  The  liquid  teiqperitnref 
hovoreTf  is  about  70**P  as  ooepared  with  about  for  the  subsonic  oasSf 
due  prlearlly  to  the  lower  pressure  level.  Free  a  consideration  of  the 
Miob  Nueber  and  teeperature  variations f  It  appears  that  the  net  stagna¬ 
tion  pressure  rise  of  soee  could  be  considerably  augnented  by  decreas¬ 
ing  the  cross-sectional  area  of  the  duct  at  about  1  foot  from  the  Inlet  In 
such  a  way  that  the  Mach  thniber  vovld  be  reduced. 


atlon*  With  regard  to  oertaln  Initial  paraaeterSf  nanely,  ^  cl^ , 

t  It  can  be  said  with  assurance  that  a  given  change  will  be 
either  undesirable  or  desirable.  S^clflcallyf  the  performanoe  will  In- 
prove  as  m,  and  d/Q  decrease  and  as  Tq^  increases  (except  at 
very  high  values  of  the  latter).  This  will  generally  be  true  irrespective 
of  the  values  of  such  other  parameters  as  ^2^  i  M|  *  “foi  * 

,  and  area  variation,  although  the  effect  on  performance  assoc¬ 
iated  with  alterations  In  any  of  the  variables  ,  and 


will,  of  course,  depend  greatly  on  the  values  of  all  the  para¬ 


meters. 


Considering  on  the  other  hand  the  p^arqneters  i  M.  »  Hi  r 

only  state  vaguely  that  whether  or  not  Increases  in 


these  will  improve  or  detiact  from  the  net  stagnation  pressure  rise 
depends  on  the  specific  values  of  these  variables  as  well  as  those 
of  {  /D|  ,  cLp  ,  X,  ,  and  of  the  area  variation.  The  reason  for  this 
indeterminancy*  as  explained  before,  is  that  an  increase  (or  decrease) 
of  any  of  the  variables  ^1  both  favor¬ 

able  and  unfavorable  aspects j  in  some  circumstsujces  the  favorable  aspects 
will  hold  the  upper  hand,  in  others  the  unfavorable  aspects  will  predom¬ 
inate  • 

These  remarks  show  that  the  effects  of  ,  M,  ,  ,  and  A( 

can  hardly  be  deduced  from  a  few  comparisons ,  and  that  hasty  conclusions 
drawn  from  a  few  cases  might  be  dangerous.  What  is  more,  a  really  mean¬ 
ingful  study  of  the  effects  of  ,  M,  $  Vq,  »  %  l^\ 

a  selection  of  the  optimum  area  variation  for  each  set  of  initial  variables 

an 

for  otherwise  a  favorable  case  might  tmfairly  be  compared  with/unfavorable 
one.  This  really  implies  that  determination  of  the  optimum  area  varia¬ 
tion  and  optimum  values  of  the  initial  parameters  must  be  made  simultan¬ 
eously.  The  studies  already  made  may  be  far  from  these  desiderata,  so 
no  attempt  will  be  made  here  to  evaluate  systematically  the  .influence 
of  ,  M|  ,  .and  /V,  . 

Regarding  ,  dj,  ,  and  T«,  ,  it  is  well  to  repeat  that  the 

magnitudes  of  the  effects  on  performance  are  in  fact  strongly  dependent 
on  all  the  other  initial  variables  and  on  the  area  variation,  although  the 

direction  of  the  effects  is  not  in  doubt.  Recognizing  in  advance  that  the 
* 

results  of  one  comparison  may  not  be  meaningful  in  other  circramstances, 
we  nevertheless  shall  present  certain  specific  results  for  the  sake  of 


llluatrating  the  influenee  of  tije  Tarlahlee  4/Dj  f  t  and 

RffiiBt  nf  Slae.  With  Yalwa  of  the  other  variables  remaining  the 
same,  Jig*  10  ■lw««  variations  in  rate  of  air  flow  on  sub¬ 

sonic,  constant-area  operation.  Tbe  parameter  “ffej  shown  In  Fig.  10 
varies  appro3d»tely  Inversely  with  the  square  root  of  the  air  flow. 

Fbr  a  value  of  In  the  nel^borhood  of  0.004  ,  the  three  values 

of  il^  correspond  to  duct  diameters  of  approximately  oo  ,  12  , 

and  2* 

The  curves  of  ter^ature,  velocity,  and  distance  vs.  fraction 
evaporated  are  nearly  Independent  ot  t  The  curves  of  p  , 

p  ,  and  M  ,  however,  show  a  strong  dependency.  As  "f/Dj  in¬ 
creases,  Regime  III  sets  In  earlier  (Regime  III  Is  absent  when  f /D|  *  0  ), 
and  the  ir«»*<ew  stagnation  pressure  Is  decreased.  A  net  rise  In 
seems  Impossible  If  exceeds  a  certain  value.  There  seems  to 

be  little  difference  between  the  cases  “f  and  f/Oj*  0.(X>4>, 

yet  the  former  would  be  very  superior  if  the  optimum  area  variation  were 

employed  fof  eaoh  oaae. 

pfAftt  of  IVonlat  Slse.  Assuming  that  the  initial  droplet  else 
may  be  varied  independently  of  ^  through  the  injector  deaign,  or  that 
the  droplet  else  produced  differs  from  that  of  the  Muklyama  and  Tanasawa 
formula  (ref.  4),  Fig.  11  ahowa  how  a  oaae  of  amall-acale,  subsonic,  con¬ 
stant-area  operation  would  be  altered  If  the  droplet  siie  were  either  one- 

I 

half  or  twice  aa  large  aa  the  Nuklyana  and  Tanaaawa  value  (  ). 

The  Influence  of  la  Been  to  be  enormous,  and  would  be  even  greater 
if  each  case  had  the  optimum  area  variation.  For  80^  h*action  evaporated. 
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crucial  to  effective  porfonnaacc  than  an  atoiai*a. 

(i)  nothing  1.  »r.  0  ^  ^ 

.  tho  octnoi  ..  -  opec- 

th<«ry  h*!  enporlnont  1. 

BffTiirt  'i»tffl’-<lr  BbUs 
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,ff.ot.t»t«-.lrr.tl.o«..P«i«“« 

„ater  ro^alrnd  to  Mtur.t.  the  olx 
.ting  thot  h.«  tho  ono-n  ,  ..n  „, 

.troM  1.  i«  «»  nolg»torh^  "  doflolonoy  and 

as  corro^-nd,  roopoPt  «  y.  i„„a„,s,  Sai!!  the  to*?  “« 

.  eonsld^ahl.  n  and  III  aaoh  set  In 

evaporation  phenomen  h  stance  ^  •  Evapora- 

,  i  .  but  at  about  the  same  distance 

et  earlier  values  »  0  ■  0.1  and  is  most  force- 

tion  urely  aohlevss  sow  prad™i  ^  ^  .trtaental 

foiforA,'®-^  •  .,t  penlomanoe  ese.s 

“  “  "T^risTconsldele  pnaotloal  ocnse,uenoe  that  the  net 
ationi  although  it  l  ID  are 

as  veil  as  the  position  of  ma=dn«  ^ 

Btagnatlon  pressure  rise 

relatively  insensitive  to  A, 


■nirough  Snaftd  Of  Sguml.  Fig.  .13a  illustrates  constant- 

» 

s^aa  operation  at  rarious  subsonlo  initial  Mach  Numbers,  and  shows  the 
extreme  eensltltity  of  the  flow  to  the  initial  Mach  Number  when  it  la 
in  the  neighborhood  of  the  critical.  To  pass  through  the  speed  of  sound, 

the  calculations  glrs  0.74Z89 . .  whereas  M,  -  0.74  gires 

purely  subsonic  operation  with  a  peak  Mach  Number  of  0.91  while  M,*0.75 
leads  to  choking  at  O.l  .  On  the  supersonic  branch  of  the  crit¬ 

ical  curve  the  rapid  rise  of  Mach  Number  soon  brings  on  Regime  III,  and  frio 

tlonal  choking  ensues  in  the  vicinity  of 

Fig.  .13b  lUustrates  the  corresponding  curves  of  gas  stagnation 
pressure.  The  corresponding  static  pressure  variations  are  shown  in 
Pig.  27,  where  they  are  ooinmured  with  experimental  data. 

phonks .  Onoe  supersonic  speeds  have  been  reached,  there  exists  the 
possibility  of  normal  discontinuities  (shocks)  in  the  gas  phase.  These 
may  occur  at  various  positions,  depending  on  the  back  pressure  level.  In 
the  calculations  it  was  assumed  that  the  gas  stream  passed  through  an 
adiabatic,  normal  shock  as  though  the  droplet  cloud  were  not  present,  and 
the  properties  of  the  droplet  cloud  on  the  downstream  side  of  the  shock 

were  assumed  the  same  as  those  on  the  upstream  side. 

Figs.  .13a  and  13b  show  operation  on  the  supersonic  critical  with 
a  shook  at  ,  the  performance  downstream  of  the  shock  being 

typically  subsonic.  Following  a  slight  drop  in  stagnation  pressure  across 
the  shook  itself,  there  is  a  rapid  rats  of  rise  brought  on  by  the  decel¬ 
eration  of  the  droplets  and  a  high  evaporation  rate  resulting  froa  com¬ 
paratively  largo  difforonoes  in  velocity  and  temperature. 
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rnntrol  of  Cross^Sectlonaltoga.  How  considerable  are  the  gains 
to  be  won  by  proper  control  of  the  cross-sectional  area  is  illustrated 
in  Fig.  U.  On  this  chart  the  base  run  for  purposes  of  comparison  is 
the  constant-area  supersonic  critical  curve  with  Xi 

.  i4.7  psla ,  -f /D,  =  0.004  and 

d  -  dl^-  *  13.4/“  ,  for  which  the  critical  Mach  Humber  is  M,  -  0.727. .... 
The  stricing  feature  of  this  constant-ai-ea  case  is  its  earlier  mentioned 
suicidal  character:  after  passing  through  Mach  Number  unity,  the  Mach 
Number  is  pushed  to  such  high  values  by  evaporative  cooling  that  the 
temperature  potential  for  heat  transfer  is  severely  depressed  and,  ailier 
about  48^6  faction  evaporated,  wall  ffiction  gains  control,  and  the  steg- 
nation  pressure  declines  rapidly  before  it  is  ever  able  to  ascend  even 

to  its  initial  value. 

By  inserting  a  normal  shock  at  ^ 

g.3  ten^.tur.  Is  sot  allowed  to  fall  too  rapidly,  sod  there  is  a  notable 
gain  in  perfomance,  the  net  stasnation  pressnre  rise  being  more  than  ll». 
With  operation  at  constant  area  on  the  subsonlo  orltloal  curve,  the 
rise  in  stagnation  pressure  is  seen  to  be  In  the  neighborhood  of 


7  or  8/t. 

In  the  light  of  the  discussions  of  Sections  4.4  and  4.5  and  of  other 
oonsiderations  too  lengthy  for  discussion  here,  the  Htch  Stunber  distribu¬ 
tion  shown  in  ng.  U  was  described,  .n  this  figure  the  solidlines  represent 
.operation  on  the  subsonlo  end  supersonic  eritloal  curves  for  a  unit  of 
about  12"-dlai.etsr.  The  variable  area  case  (dashed  line)  is  compounded 
of  (1)  constant-area  operation  on  the  supersonic  crltlonl  up  to  ]C’^0.y>  , 

(11)  variable  area  operation  from  t-'O.'SO  to  with  the 

Mach  Humber  distribution  Illustrated,  and  (ill)  operation  in  a  conical 
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dlffuBer  of  e^includod  angle  for  X'>0.54-  *  It  may  be  seen  that 
the  maximum  area  reduction  ie  25$  of  the  initial  area,  and  that  the 
diffuoer  begins  shortly  after  the  position  of  minimum  area,  Vfhat  ie 
of  greatest  practical  importance  is  that,  in  the  variable-area  case  ehown, 
a  stagnation  pressure  rise  of  more  than  18$  is  indicated. 

The  foregoing  example  is  only  an  isolated  instance  presented  mere¬ 
ly  to  Bhow  the  importance  of  Mach  Number  control,  and  is  illustrative  of 
the  improvements  to  be  gained  whether  this  control  is  brought  about  by 
area  control  or  by  lengthwise  water  injection.  In  this  example,  the 
Mach  Number  variation  was  chosen  arbitrarily  and  would  be  the  optimum 
only  through  coincidence.  As  mentioned  earlier,  more  systematic  ways 
of  determining  the  optimum  are  under  development.  Only  fl-aguientary 
results  of  this  more  systematic  study  are  available,  the  most  iiiqportant 
practical  conclusion  thus  far  being  that  with  >01  “  ^  and  an 

air  flow  of  about  20  Ib/sec,  a  stagnation  pressure  rise  of  more  than 
20%  seems  to  bo  attainable  provided  that  the  drop  size  is  that  of 
the  Nukiyaraa-Tanasawa  forimi]-a  (ref.  4)* 

4.8  Approximate  Discontinuity  AnalyBis 

For  real  Insight  into  the  details  of  the  Aerothermopi-essor  process, 
01*  for  knowledgeable  optimization  of  the  area  vai’iation,  wa  have  no 
better  procedure  at  the  moment  than  the  extended  numerical  integrations 
described  in  Section  4.7.  But  for  the  approximate  computation  of  the 
order  of  magnitude  of  attainable  performaiico,  a  compai-atively  simple 
algebraic  method  nay  be  constructed. 

Rather  than  dealing  accurately  with  the  history  of  the  process, 
the  approximate  analysis  proceeds  on  the  assumption  that  in  the  evapora¬ 
tion  section  the  curve  of  preseiu-e  vs,  area  may  be  represented  as  a 


-  64 


straight  lino.  The  core  of  the  analysis  is  then  the  Integral  equations 
of  continuity,  iDoinentua,  and  energy  oonneoting  the  inlet  and  exit  states 
of  the  evaporation  section.  At  low  Mach  RubAmts  an  especially  single 
form  of  the  analysis  becomes  usabls. 

Details  of  the  analysis  are  recorded  in  Appendix  D,  and  a  typical 
set  of  counted  results  is  shown  in  Fig.  15.  Mure  extended  analyses 
of  the  discontinuity  type  together  with  «  wealth  of  nunerloal  results 
are  presented  in  ref.  8. 

In  Fig.  15  the  solid  ouxves  represent  the  low  Maoh  Number  analysis. 
For  eomparison,  the  dashed  curve  shows  how,  for  the  ease  of  no  iViction 
in  the  evaporation  section  and  diffuser,  the  results  are  modified  idien 
the  more  accurate  iterative  oaloulation  procedure  of  Appendix  D  is  used* 
With  the  inlet  parameters  shown,  an  abscissa  of  1.0  of  course  corresponds 
to  constant  areai  1.24  corresponds  to  equal  Kach  Number  at  inlet  and 
exit;  and  1.51  corresponds  to  equal  pressures  at  inlet  and  exit. 

A  glance  at  the  resulta  of  Section  4>7  shows  that  the  linear  prea.. 
3ure->area  relatlonahlp  Is  a  orudo  approximation.  Nevertheless  it  is 
believed  that  aa  far  as  gross  orders  of  magnitude  are  concerned,  the 
curves  of  Fig.  15  are  signifioant  in  respect  to  the  attainable  stag¬ 
nation  pressure  rise  and  In  dlBplaylrjg  the  effects  of  artta  change, 
wall  friction,  and  diffuser  loss.  Among  the  oonolusions  that  my  bo 
drawn  are  the  following; 

(1)  With  a  diffuser  loss  coefficient  of  0,1^ a  dlamater  of  1  foot, 
a  leiiixth  of  6  feet,  and  f  *  0.004  ,  a  stagnation  presaui^a  rine 

approaching  Z0%  is  ppBsihla, 

(ii)  A  net  rise  in  stagnation  preeaure  is  not  posaiila  if  fL/D 
axcaeds  approximately  0,075  (for  i  ^  0.005 


,  thiB  linilte  L/p  to  15 ). 


(Ill)  Area  decrease  la  definitely  beneficial,  and  the  optimum 
amount  of  area  contraction  Increases  as  the  frictional  effects  are 
minimised. 


1.,  EXPERIMENTS  WITH  A  SMALUSCALEt 
gONSTAMT.AREA  AEROTHFSMQPRPAqnR 

— ShiSStiYt  of  awn-iSoale 

Tlia  experiments  described  below  were  carried  out  on  a  scale  of  air 
flow  so  small  that  there  was  no  expectation  of  a  net  rise  in  stagnation 
pressinre.  Instead  the  goals  were  as  follows: 

(I)  To  explore  quickly  and  economically  the  diverse  phenomena  and 
regimee  of  operation  in  the  Aerothermopreetior, 

(II)  To  obtain  e^qjerimental  data  for  testing  the  validity  of  the 
theoretical  caloulationa . 

(III)  To  estimate  by  extrapolation  what  would  be  the  prospects  of 
effective  operation  on  a  large  scale,  with  the  purpose  of  assessing  the 
advisability  of  an  expensive,  medlum-eoale  experiment, 

(Iv)  To  acquire  informatJon  ncceaaary  for  the  Intelligent  design 
of  such  a  medium-Boalo  unit. 

All  these  aims  were  achieved.  The  methode  and  results  are  described 

In  ful]  detail  by  Wadleigh  (ref.  9),  and  a  summary  of  this  work  Is  presented 
here. 

5.aZ  Deta^p  9^'  ApparatuB 

Powered  by  a  steain-jot  ejector  at  the  downstream  end, 
■atire  system  (Fig.  16)  operated  at  Bub~atmoHphorie  prenyuie,  Air 
ft-om  the  atmoephere  was  drewn  into  a  Eji,0-rirod  combuetion  chamber j  from 
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the  letter  the  heated  air  paaaed  through  tvo  uilring  elbowa  ead  a  etralght- 
enlng  length  to  the  upstream  stagnation  seotioni  thence  through  a  con¬ 
verging,  accelerating  oossle  Into  the  evaporation  duct,  and,  after  being 
cooled  by  vater  «,apor.tlon,  Into  a  conical  diffuser  of  S'  included 
angle,  finally,  the  air  floeed  through  a  largo  quench  tank  to  the  Jet 
pnsp  Vi.  a  control  valve  by  means  of  uhich  the  back  pressure  on  the  Aero- 

therinopressor  was  adjusted. 

T..t  l>.ct.  AS  the  objective  »as  to  obtain  basic  data  rather  than 
maxlmom  p«.formB«.,  the  svaporatlon  duct  «s  of  constant  area,  2.13 
in  *es  in  dlMsster  and  73  Indhss  long,  the  test  section  vas  fltt«i 
with  tventy  statl.  pr.«»r.  taps  on  the  sail,  and  sltb  eight  access  port. 

for  introduction  of  measuring  probes. 

uea.e  TCectlon  System.  Atomisation  of  the  so-called  "gas  atomisa- 

tlon"  typo  sas  achlnred  by  introducing  sater  at  low  speed  into  the  high- 
velocity  air  stream  at  the  beginning  of  the  ^13  -inch  evaporation  duct. 
Since  the  type  of  water  injection  system  which  would  beat  yield  the  desir¬ 
able  features  of  droplets  together  with  uniform  spatial  distribu¬ 

tion  could  hardly  be  predicted  in  adv«.ce,  prellmln«fy  emperlments  were 
osrrled  out  with  some  twenty  different  injection  schemes  (ref.  9),  in¬ 
volving  oombinatlons  of  axial-flow  water  Jeta,  radial-flow  water  Jets 
originating  in  the  pip.  wall,  and  water-sheet  formation  nossles.  The 
merit  of  each  injection  eyst»n  was  Judged  by  compsirlng  the  net  change 
in  stagnation  presaurs  fPcn  upstrea.  stagnation  section  to  downstrMm 
quench  tad,,  with  all  other  variable,  at  the  inlet  plane  of  the  evap¬ 
oration  duct  held  constant.  Very  mibstantlal  dlfferenoee  in  pa^formanoe 
were  found.  Ihe  beet  injection  system,  and  that  used  for  all  further 
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vork|  va«  an  axial  injaotlon  aystan  In  which  the  water  waa  Introduced 
at  the  inlet  plane  of  the  S.13  >lnoh  duct  neana  of  aeven  tubea 

of  dinanalona  0.125 *-0D,  0.110  -ID  ,  with  one  tube  on  the  axia 
and  the  othera  diapoaed  ayumetrioally  on  a  circle  of  1.25  -inohea  dianeter* 
Manaurementa.  Ihe  following  maaaureinenta  wera  nadet 

(i)  Stagnation  preaaure  and  atagnatlon  taaq^ature  in  the  up- 
atreaa  atagnatlon  aeetlon  ware  neaaurad  directly  by  an  iopaot  tuba  and 
a  multl-ahieldad  thermocouple. 

(ii)  Stagnation  inreaaure  at  the  diffuaer  exit  waa  naaaured  an 
iiqiact  tuba. 

(lii)  Air  flow  rate  waa  inferred  from  the  upstream  atagnatlon  prop¬ 
erties  and  the  preaaure  drop  aoroaa  the  inlet  bell  mouth,  aaaumlng  for  ^e 
latter  a  discharge  coefficient  of  unity. 

(iv)  Water  flow  waa  naaaured  by  a  rotametar  and  a  aharp-edged  or¬ 
ifice, 

(v)  Entranoe  air  Mach  Number  (at  the  plane  of  water  injection)  waa 
calculated  fToia  the  iaentroplo  law  connecting  it  with  the  upstream  atagna¬ 
tlon  pressure  and  the  pressure  drop  across  the  inlet  ball  mouth. 

(t1)  Longitudinal  static  pressure  distributions  along  the  length 
of  the  test  section  were  determined  by  connecting  the  wall  static  taps  to 
a  multiple  manomater  board. 

(vii)  The  amount  of  water  evaporated  at  the  sections  where  access 
porta  Were  available  was  determined  by  traversing  these  sections  with  a 
special  sampling  probe  designed  for  withdrawing  a  san^le  of  the  gas  phase 
only  and  by  measuring  the  humidity  of  the  air-vapor  aan^le  by  means  of  a 
Foxboro  "Dewcel'^  indicator.  Details  of  the  humidity-measuring  system  and 


of  the  novel  eanpllng  probe  are  given  in  ref*  9* 

Control  of  Varlablea.  In  the  apparatue  described,  the  slgnifl^t 
parameters  which  could  be  controlledi  were  (i)  the  upstream  stagnation 
pressure,  *1^  (il)  the  upstream  stagnation  temperatu^'s,  } 

(ill)  the  water-air  ratio,  12^  and  (Iv)  either  the  initial  Mach  Number, 
Mj  ,  or,  in  the  case  of  choked  flow  where  was  uniquely  deter¬ 
mined  by  ,  and  Q.^  ,  the  position  of  a  normal  shock  in 

the  test  sedtion* 

These  were  varied  in  a  systematic  manner  during  the  course  of  the 
tests . 

5t3 . 

We  present  here  selected  results  which  Illustrate  the  theoretical 
considerations  of  Section  4  and  which  also  show  the  Influences  of  the 
controllable  parameters.  Here  again  the  reader  is  cautioned  that  the 
effeota  of  these  parameters  in  a  small-scale,  constant-area  apparatus 
are  generally  not  the  same  as  in  a  large-scale  unit  having  an  optimum 
area  variation. 

Comparison  of  Flows  With  and  Without  Evaporation,  figure  1?  Ulus- 
trates  in  an  unmistakable  way  the  powerful  influence  of  that  evaporative 
cooling  which  is  the  heart  of  the  Aerothermopressor  process. 

The  throe  runs  with  1^,  "  I500*R  constitute  essentially  adiabatic 
flows.  With  no  water  injection  (  X2q  *  0  )  the  process  is  the 

one-dimensional  Fanno-line  process  (ref.  2^  chapter  6)  and  exhibits  a 
constantly  falling  static  pressure  typical  or  subsonic  flow,  with  the 


rate  of  pressure  decrease  increasing  because  the  Mach  Number  is 
iiicr-eased  by  friction.  As  water  is  introduced,  with  12^  first 
equal  to  0.1  and  then  O.Z.  ,  a  larger  drop  in  pressure  than  before 
occurs  near  the  inlet  because  of  droplet  dragj  after  the  first  few 
inches,  the  droplet  drag  is  small,  and  the  pressure  drop  is  primarily 
that  of  pipe  friction.  However,  the  rate  of  frictional  pressure  drop 
is  greater  at  the  larger  values  of  12^  ,  because  tho  droplet  drag  in¬ 
creases  the  Mach  Numbers  at  which  friction  occurs  throughout  the  duct. 

For  the  case  Tqi  *“530f2  and  2  ^  in  fact,  the  flow 

at  the  duct  exit  is  nesirly  choked,  i.e.  the  I^ch  Number  is  nearly  unity. 

The  case  of  ^  ■■  1500*12.  is  again  essentially  adia¬ 
batic  and  is  directly  comparable  with  the  case  of  at  1^|*530 

The  small  difference  in  the  rates  of  pressure  drop  for  these  two  Fanno- 
llne  runs  is  present  because  of  a  difference  in  initial  Mach  Number  and 
also  because  the  pipe  friction  factor  f  is  altered  by  the  difference  in 
length  Reynolds  Number  between  the  two  cases. 

Turning  to  the  curves  for  “  1500  K  with  water  injection, 

we  see  pronounced  differences  owing  to  the  evaporative  cooling  which  is 
now  present.  Whereas  for  1^,“’5S0‘‘R  (i.e.,  no  evaporation)  water  in¬ 
jection  Increases  the  pressure  drop,  with  Xi  **  ^  (1»®»»  strong 

evaporation)  water  injection  decreases  the  pressure  drop.  It  seams  safe  to  say 
also  that  the  Mach  Number  rise  here  when  evaporation  occurs  is  less  than 
when  evaporation  does  not  occur. 

Effect  of  Water-Air  Ratio,.  Certain  aspects  of  the  phenomena  Just 
discussed  are  brought  out  more  clearly  in  Figure  18,  for  which  the  initial 
Mach  Number  of  0.65  is  higher  than  that  for  the  runs  shown  in  Figure  16, 
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thus  magnifying  the  pressure  gradients.  The  fact  that  cooling  at  sub¬ 
sonic  speeds  does  indeed  tend  to  decrease  the  Mach  Number  was  elearljr 
evidenced  by  this  series  of  runs,  inasmuch  as  it  was  not  even  possible 
to  operate  at  «0.6?  with  no  water  injection,  the  initial  Mach 
Number  for  choking  idien  operating  dry  being  less  than  0.&5  , 

For  each  y^ve  of  ,  the  three  characteristic  Regimes  dis¬ 

cussed  in  Section  4*3  are  observed: 

(1)  Regime  I.  For  a  distance  of  about  two  inches  the  pressure  falls 
rapidly  because  of  droplet  drag.  The  magnitude  of  this  initial  pressure 
drop  is  roughly  proportional  to  the  rate  of  water  flow,  i.e.  to  . 

(il)  RwHem  tt.  a  rapidly  rising  pressure  occurs  in  the  next  20  Inches 
or  so,  where  evaporative  cooling  is  predominant.  The  maximum  rate  of 
pressure  rise  increases  as  increases,  oMefly  because  an  increased  12^ 

means  a  greater  amount  of  heat  transfer  surface  in  the  droplet  cloud.  For 
such  a  anell  £L^  as  0.10  ,  the  net  rise  is  rather  small, and  reaches 

its  peak  at  <  ^  10  inches}  in  this  case  there  is  scarcely  enough  eurface 
to  accomplish  much  cooling,  and  the  droplet  cloud  vanishes  comparatively 
rapidly  because  the  driving  tenperature  differential  is  not  decreased 
very  much  by  cooling  of  the  gas.  With  values  of  between  O.Z 

and  0.4  the  pressure  rise  exceeds  the  loss  in  Regime  I,  such  that  at 
f  S  25  inches  the  pressure  is  hl^er  than  in  Ihe  inlet  plame.  The 
nwytimiin  pressure  in  the  duct  and  also  at  the  dlffhser  exit  is  obtained 
with  ^  0.3  y  this  being  tiie  water-air  ratio  for  which  the  balance 
between  pressure  rise  due  to  cooling  and  pressure  drop  due  to  droplet  accel¬ 
eration  and  wall  friction  is  an  optimum.  For  low  valuea  of  ,  the 
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dlBtance  ^  to  tho  point  of  maxlmizn  preaeure  increases  as  £1^ 

increases  because  the  droplet  cloud  persists  for  greeter  distances.  For  _ 

high  values  of  ,  on  the  other  hand,  the  distance  decreases  as  ^2,^ 

increases  because  the  gas  stream  tends  towards  saturation  awre  quickly.  j 

(ill)  ttt.  After  the  rate  of  evaporation  has  been  reduced 

Buffioiently  by  the  falling  gas  tenqisrature  and  diminishing  droplet  i 

diameter,  wall  friction  gains  control.  The  pressure  begins  to  drop, 
and  the  rate  of  pressure  fall  grows  as  evaporation  becomes  weaker  and 

'i 

weaker . 

Passage  Throu>rh  Mach  Number  Ifeltv.  Figure  W  shows  the  pressure  ^ 

distributions  obtained  with  constant  values  of  Tm  )  f  and  ,  | 

*  I 

but  with  varying  initial  Mach  Numbers,  ^ 

The  curves  for  M,  “0.48  and  M,“0.65  are  for  values  of  back  | 

pressure  so  high  that  the  entire  flow  is  subsonic.  Qualitatively  | 

similar  to  those  of  Fig.  18,  they  demonstrate  as  wall  the  dynamic  i 

character  of  the  events  in  the  sense  that  the  pressure  gradients  grow  ^ 

as  the  Mach  Number  ascends*  ( 

I 

All  the  remaining  curves  of  Fig*  19  are  for  the  critical  initial 
kfach  Number,  M,  =0.78  ,  for  which  evaporation  becomes  predominant 

over  drag  .lust  at  the  point  where  the  local  Mach  Number  is  unity.  The 
lowermost  curve  in  Fig.  19  is  for  the  greatest  possible  length  of  super¬ 
sonic  now  in  the  test  section,  and,  except  for  a  shock  near  the  exit, 
corresponds  to  the  supersonic  critical  curve  shown  in  Fig.  13.  The 
branch  curves  rising  precipitously  from  this  lowermost  curve  represent 
normal  shocks  followed  by  subsonic  flow. 
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ShnnV  Configurations.  Considering  In  more  detail  the  lowermost 
curve  in  Pig.  19,  it  la  seen  that  the  pressure  falls  and  the  Mach  Bumber 
rises  to  unity  in  Regime  I,  up  to  Z  Inches;  the  pressure  continues 
to  fall  and  the  Mach  Numcor  rises  to  substantial  supersonic  values  in 
Regime  II,  up  to  ^  *  28  inches;  then,  because  of  wall  ftriction, 
the  pressure  begins  to  rise  and  the  Mach  Humber  to  fall  in  Regime  III, 
up  to  *  -  55  inches.  However,  the  duct  length  is  too  great  for 
supersonic  flow  to  continue  to  the  exit,  for  in  the  absence  of  shocks, 
the  Mach  Number  would  reach  the  choking  value  of  unity  before  the  exit. 
Therefore,  a  shook  must  stand  somewhere  in  the  duct.  For  very  low  back 
pressures  (e.g.,  condition  2  in  Fig.  19),  the  shock  is  at  its  further¬ 
most  possible  position  downstream,  that  is,  at  the  position  for  which 
the  subsonic  flow  following  the  shock  becomes  choker^  at  the  duct  exit. 

In  this  ease  the  shock  lies  in  the  range  of  «  between  55  and  62  Inches, 
and  M»i  at  the  duct  exit;  the  flow  in  the  diverging  diffuser  (now 
operating  as  a  nozzle)  then  accelerates  to  supersonic  speeds,  and  is 
again  decelerated  to  subsonic  speeds  by  means  of  a  shock  in  the  "dif¬ 
fuser",  the  shock  having  a  loss  in  stagnation  pressure  sufficiently 
high  to  match  the  prescribed  low  back  pressure.  The  lowermost  cuiiro 
therefore  exhibits  an  amazing  range  of  gas  dynamics  phenomena,  com¬ 
prising  two  continuous  accelerations  through  the  speed  of  sound  (albeit 
by  different  mwchanisms),  two  normal  compression  shocks,  three  distinct 
regions  of  subsonic  flow,  and  two  distinct  regions  of  supersonic  flow. 

As  the  back  pressure  is  increased,  say  from  condition  2  to  condi¬ 
tion  3,  the  shock  in  the  "diffuser"  moves  upstream  to  lower  supersonic 
Mach  Numbers,  but  the  flow  in  the  duct  is  completely  unaltered.  When 


the  back  preaaure  reaches  condition  however,  the  "diffuser"  shock 
disappears  In  the  duct  exit  where  M  »  1,  and  this  then  represents  the 
limit  of  operation  with  choked  flow  at  the  duct  exit,  and  with  the 
furthermost  downstream  position  of  the  shook  In  the  duct. 

A  further  rise  in  back  pressure  over  that  for  condition  4  there¬ 
fore  causes  the  shock  in  the  duct  to  move  upstream*  As  the  back  pres¬ 
sure  is  increased  the  shock  continuously  and  suooesslvely  moves  from 
"  position  1  to  positions,  5,  6,  7  and  8.  For  all  the  latter  modes  of 
operation  the  flow  following  the  shook  and  at  the  duct  exit  is  subson¬ 
ic.  Following  the  shook,  the  pressure  at  first  rises  because  of  subsonic 
cooling  and  deceleration  of  the  liquid  droplets,  and  subsequently  falls 
because  of  subsonic  friction. 

Because  of  the  now-famlllar  shock-boundary  layer  Interaction,  the 
shocks,  rather  than  being  disoontinuities,  occupy  a  distance  of  about 
8  inches,  corresponding  to  about  4  duct  diameters. 

Stable  vs.  Unstable  Operation  at  the  Critical  Mach  Number.  Attempts 
to  operate  on  the  subsonlo  branch  of  the  critical  curve  or  with  shocks 
further  upstream  than  the  position  corresponding  to  8  were  unsuccess¬ 

ful,  When  the  back  pressure  was  raised  slightly  above  that  of  condition  8 
(corresponding  to  curve  8),  the  entire  flow  became  highly  unsteady,  with 
large  surges  in  pressure,  flow  rate,  and  shock  position. 

This  strange  behavior  may  be  explained  by  reference  to  the  schematic 
diagram  of  Fig.  20,  which  refers  to  a  small-soale  Aerothermopressor  In 
which  there  is  at  all  times  a  net  loss  in  stagnation  pressure  from  inlet 

to  diffuser  exit.  Consider  an  experimental  arrangement  in  which  P  , 

'oi 

,  and  are  fixed  and  in  which,  at  least  for  the  present,  steady 


flov  Is  assiuaed  possible*  Theiii  ss  Increases  from  serOf.  ;tbe 

value  of  P  at  the  diffuser  exit  will  at  first  decrease  because  the 

«  Z 

net  loss  in  Pq  io  roughly  proportional  to  ,  This  loss  contin¬ 

ues  until  Mj  reaches  its  critical  value,  after  which  M,  nay  increase 
no  further,  and  wo  nay  use  the  position  of  the  shock  as  the  controllable 
variable  governing  perfomanoe*  Fkon  the  theoretical  consideration  of 
Section  4,  there  will  be  a  particular  position  of  the  shock  (point  a) 
in  Fig.  20  for  boot  perfomanoe.  In  Fig*  19,  therefore,  the  curve 
of  -p  versus  shock  position  will  first  rise  to  point  a,  and  will 
subsequently  fall  as  the  shock  moves  further  downstream. 

Now,  operation  at  a  point  to  the  rig^t  of  a,  ouch  as  at  b,  is  stable 
in  the  sense  that  a  slight  variation  in  the  jiyiposed  back  pressure  will 
induce  a  change  in  Aerothermopresoor  operation  of  such  nature  as  to  make 
the  pressure  adjust  to  the  iaqjosed  back  pressure.  For 

example,  suppose  that  steady  operation  prevails  at  b,  and  that  the  in^ 
posed  back  pressure  is  suddenly  increased  slightly*  This  will  gener¬ 
ate  a  ymell  compression  wave  which  will  propagate  upstream  until  it 
reaches  the  shock,  and  its  interactiosi  with  the  shook  will  force  the 
latter  slightly  upstream.  But  the  latter  event  Increases  the  delivered 
back  pressure  and  thus  the  Aerothermopressor  adjusts  stably  to  the  imposed 
dlstxirbanoe. 

Operation  at  a  point  such  as  c  on  the  portion  of  the  curve  between 
a  and  e,  on  the  other  hand,  is  unstable.  If,  after  an  assumed  steadf 
operation  at  c,  the  imposed  back  pressure  is  suddenly  raised  slightly, 
the  generated  compression  wave  will  on  reaching  the  shock  again  move 


the  latter  upstream.  But  this  reduces  the  delivered  exit  pressure, 
with  the  result  that  conditions  are  further  ft-om  equlllbriuin  than 
before,  and  the  point  of  operation  is  driven  toward  e.  There  ensues 
an  oscillatory  mode  of  lare®  amplitude,  the  character  of  which  depends 
upon  the  transient  characteristics  of  the  air  and  water  now  systems, 
that  is,  on  the  way  in  which  »  Xi  »  ^  react  to  the 

changes  in  air  How  pursuant  to  M,  dropping  below  its  critical  value. 

These  phenomena  are  reminiscent  of  the  behavior  of  conr^erging-diverging 
diffusere  for  supersonic  now.  It  is  well  known  that  the  operation  of 
such  a  diffuser  is  stable  when  the  shock  is  in  the  divergent  portion 
(in  which  region  a  downstream  displacement  of  the  shock  reduces  the 
delivered  stagnation  pressure),  but  unstable  when  the  shock  is  in  the 
convergent  section  (in  which  region  a  downstream  displacement  of  the 
shock  increases  the  delivered  stagnation  pressure). 

of  Performance  Coeffloiegt.  Inasmuch  as  all  the  phenom- 

ena  occurring  within  the  Aerothermopressor  are  strongly  dynamic  in  char¬ 
acter,  a  measure  of  overall  performance  is  required  from  which  most  of 
the  sensitivity  to  Ifcch  Number  has  been  removed.  Being  mindful  of  the 
Corm  of  th.  Influonoe  coofflolento  for  dp^  /p^  ohown  in  TaUo  HI. 
we  form  the  over-all  performance  coefficient 
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where  P  le  the  stagnation  pressure  at  the  diffuser  exit* 

Althou^  the  use  of  this  dimensionless  figure  of  performance 
does  in  fact  illuminate  the  results.  It  seems  in  place  once  again  to 
state  that  the  value  of  this  performance  coefficient  depends  not  only 

on  p  f  f  -^0  *  ^1  *  length  of  evaporation  section, 

*01 

but  also  on  the  actual  duct  diameter  and  on  the  area  variation,  becausei 
both  friction  and  area  change  act  to  alter  the  local  Mach  Number. 

Effect  of  Ratio  on  Qyer-All  Performance.  The  curves  of 

Fig.  21  are  typical  of  the  variation  in  o\i'er-all  performance  when  the 
water  injection  rate  is  varied  and  the  initial  stagnation  conditions 
and  Mach  Number  all  held  constant*  For  each  set  of  initial  parameters 
( K.  X.,  ?„  ),  there  is  an  optimum  value  of  .  Below  this 

optimum  ,  the  total  amount  of  cooling  is  weak  compared  with  wall 

friction;  above  this  value,  the  added  droplet  drag  cannot  be  offset  by 
evaporative  cooling  because  of  the  tendency  of  the  stream  to  reach  satur¬ 
ation  too  quickly.  The  performance  is  less  sensitive  to  variations  in 
above  the  optimum  value,  however,  than  to  variations  below  the  optimum, 
i,e.,  eMessive  water  is  not  as  detrimental  to  performance  as  a  deficiency. 
The  curves  are  fortunately  quite  flat  near  their  peaks,  and  so  the  pre¬ 
cise  adjustment  of  the  water  flow  is  not  critical. 

As  might  bo  expected  from  the  theoretical  considerations,  an  increase 
in  1^,  requires  an  increase  in  ^  •  For  example ,  with  "p  «  1 4-.  8 

psla,  the  table  below  shows  the  optimum  values  of  at  several  temper¬ 

atures  { 


^01' 


1200®R 


ISOO^R 


1800®R 


0.16 


0.24 


0.31 


Figure  21  ehows  also  that  the  optimum  water-air  ratio  is 

substantially  independent  of  0  . 

■01 

No  significant  trend  of  optimum  water-air  ratio  with  InitiSl 
Mach  Number  was  observed. 

Effect  of  Initial  Stagnation  Temperature  on  Over-all  Performanoe*. 
from  Fig.  21  it  may  be  seen  that,  with  a  given  initial  Mach  Number,  the 
best  performance  attainable  with  the  appropriate  value  of  12^  is 
improved  as  the  initial  stagnation  temperature  is  raised.  The  same 
result  is  obvious  also  from  Fig.  22,  which  shows  only  the  best  per¬ 
formance  at  each  Mach  Number  and  temperature.  An  Increase  in  the 
temperature  difference  available  for  cooling,  as  well  as  an  Increase 
in  the  saturation  humidity,  account  in  a  simple  manner  for  the  better 
performance  at  high  temperatures.  On  the  same  grounds  it  is  easy  to 
understand  also  why  the  optimum  water-air  ratio  rises  as  the  Initial 
stagnation  temperature  increases,  at  least  in  the  range  of  the  latter 
under  consideration. 

Effect  of  Initial  Mach  Number  on  Over-AJ,!  Feyfpymance,  Natiirally, 
the  net  loss  in  stagnation  pressure  in  the  small-scale  Aerothermopressor 
increases  as  the  Mach  Number  level  is  Increased.  However,  this  gives 
a  false  picture  of  how  the  performance  of  a  large-scale  unit  might  depend 
upon  initial  Mach  Number,  A  better  representation  is  given  in  Figure  22, 
where  the  over-all  performance  coefficient  is  plotted  against  initial 
Mach  Number.  The  right-hand  end  of  each  curve  represents  the  choking 
Mach  Number  for  the  specified  values  of  *  loi  *  * 

the  corresponding  value  of  the  ordinate  is  for  the  shock  position  yielding 


the  maxlioum  final  sta^atioh  pressure* 


What  Figure  22  shows  la  that  the  over-all  performance  Improves 
as  the  initial  Mach  Number  increases.  If  the  wall  fflction  were  some¬ 
what  less  than  in  the  case  here  considered,  curves  such  as  those  in 
Figure  22  would  eJLso  have  negative  values  of  the  ordinate  at  low  Mach 
Numbers j  however,  they  would  differ  in  that  they  would  have  positive 
values  at  high  Mach  Numbers.  Thus  there  is  in  general  a  certain  in¬ 
itial  Mach  Number  below  which  a  rise  in  V  cannot  be  achieved,  and 

'0 

above  which  the  net  rise  in  stagnation  pressure  mounts  rapidly  as  the 

initial  Mach  Number  is  increased.  The  explanation  of  this  behavior 

is  related  to  the  effect  of  M,  on  droplet  diameter;  as  M|  increases, 

a 

the  droplet  diameter  decreases,  and/greater  fraction  of  the  water  is 
evaporated  in  the  given  length  of  duct,  as  described  previously  in  Sec¬ 
tion  4*6  and  as  shown  later  in  Figure  26. 

Effect  of  Initial  Stagnation  Fraasure  on  Over-all  Performance,.  FigureJ 
21  and  22  also  show  a  considerable  worsening  in  perfoivsance  as  the  inlet 
stagnation  pressure  is  reduced.  For  fixed  values  of  M  ,  \  ,  and 

*  It 

12^  ,  the  influence  of  decreasing  may  be  felt  in  several  ways, 

namely;  (1)  through  an  increase  in  initial  droplet  diameter  produced  by 
decreased  atomization  forces;  (ii)  through  a  reduction  in  heat  trans¬ 
fer  coefficient  and  in  drag  effects  produced  by  lessened  relative  Rey¬ 
nolds  Numbers;  (Hi)  through  an  increase  in  pipe  friction  factor  pro¬ 
duced  by  the  decrease  In  pipe  Reynolds  Number;  (Iv)  through  the  lowered 
wet-bulb  temperature  and  the  Increased  vapor-carrying  capacity  of  the 


stream  before  satiaration  is  reached?  and  (v)  through  an  inorease  in  the 
mass  transfer  coefficient  associated  with  an  inorease  in  the  mass  diffus- 

ivlty. 

As  heat  transfer  rather  than  mass  transfer  is  controlling  in  the 
evaporation  process  effects  (Iv)and  (v)  are  probably  of  small  consequence. 
Furthermore,  it  is  hard  to  Imagine  effect  (iil)  as  accounting  for  the 
results  seen  in  Figs.  21  and  22. 

Some  additional  relevant  facts  appear  in  Big.  23,  which  shows  the 
longitudinal  pressure  distributions  at  two  pressure  levels  with  compar¬ 
able  (but  slightly  different)  values  of  M|  and  •  Taking  note 
of  the  comparative  values  of  M,  and  ,  what  is  most  striking  Is 
the  enormous  difference  in  the  rates  of  pressuwrise  at  the  beginning 
of  Regime  II.  In  the  theoretical  calculations  of  ref.  8,  the  effect 

of  P  under  conditions  of  unchanged  was  found  to  be  extremely 
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email,  i.e.,  the  net  influence  of  item  (li)  in  the  foregoing  list  was 
found  to  bo  unimportant.  Accordingly  the  evidence  points  to  item  (i)  aa 
being  controlling,  i.e.,  that  the  effect  of  initial  stagnation  pressure 

is  strongest  through  its  influence  on  drop  size. 

Evaporation  Profiles.  For  fixed  Initial  conditions,  Fig.  2/.  shows 
typical  profiles  of  specific  humidity  at  four  different  transveree  cross- 
sections  of  the  flow. 

Near  the  inlet,  at  *  =  4.5  ,  the  humidity  profile  reflects  the 

deficiencies  of  the  water  injection  system  in  achieving  uniform  distri¬ 
bution.  However,  even  in  the  short  distance  of  4.5  Inches,  an  astonish¬ 
ing  amount  of  evaporation  has  occurred,  thus  testifying  to  the  tremendous 


heat  trrnsfar  surface  of  a  droplet  cloud. 

Further  dovmatream,  a  rapid  straightening  of  the  hum5dity  profile  is 
observed,  until  at  £“€^5  the  profile  is  nearly  flat.  By  virtue  of 
the  large  droplet  inertia,  it  does  not  seem  likely  that  the  turbulence 
of  the  stream  does  much  to  redistribute  the  droplets  in  a  more  uniform 
way.  What  the  data  suggest,  rather,  is  that  the  turbulent  migrations 
of  the  gas  Itself  rapidly  mixes  those  portions  of  the  gas  in  contact  with 
many  droplets  with  t'lose  portions  in  contact  with  few.  Or  from  a  differ¬ 
ent  point  of  view,  each  piece  of  gas  courses  turhulently  through  the  dif¬ 
ferent  parts  of  the  droplet  cloud,  and  thus  experiences  the  average  wat¬ 
er  distribution.  If  this  picture  is  true,  uniformity  of  the  spray  dis¬ 
tribution  on  a  flne-grairod  scale  may  not  be  critically  Important,  so 
long  as  uniformity  on  a  larga-grained  scale  is  achieved.  Vfeight  is  added 
to  this  picture  of  highly  turbulent  gas  motion  by  the  observation  in 
Pig.  24  that  the  gas  humidity  at  71 =  0.5  inches  docreases  as  the  gas  goes 
from  z=«13.5  in,  to  z  =  40.5  inches;  this  can  only  be  accounted  for  by 
the  mixing  of  the  very  humid  gas  at  71 «  0.5  inches  with  the  much  less 
humid  gas  near  the  walls. 

Ultimately , the  humidity  is  higher  netir  the  walls  than  at  the  center. 
This  Booms  to  be  the  result  of  the  velocity  profile  in  the  duct,  which 
grants  to  droplets  near  the  walls  a  longer  residence  time  than  to  those 
at  the  center. 

Totig-itudinal  Rate  of  Evauoration.  By  integrating  humidity  profiles 
such  as  those  of  Fig.  24,  VJO  obtain  the  open-point  cxirves  of  Fig.  2$, 
shewing  row  average  humidity  increases  with  longitudinal  distance.  There 


IB  OTident  here  the  very  rapid  Initial  rate  of  evaporation  followed 
by  a.  decreasing  rate  which  continues  to  fall  as  the  stream  tends  toward 
saturation.  At  a  given  duct  position,  the  specific  hui-idlty  increases 
as  the  water  injection  rate  is  Increased.  The  fraction  evaporated,  how¬ 
ever,  decreases  as  SI,  increases  because  the  increased  water  concen¬ 
tration  causes  the  stream  to  approach  saturation  more  rapidly. 

For  the  case  0,-0  ,5056  of  the  water  is  evaporated  in 

the  first  10  inches,  and  9556  is  evaporated  in  a  distance  of  40  inches. 

For  the  case  which  is  slightly  greater  than  the  saturation 

value  of  about  0.21,  approximately  50j6  of  the  saturation  amount  is  evap¬ 
orated  in  the  first  10  inches,  and  about  85^  of  the  saturation  amount 
in  40  inches. 

In  addition  to  the  direct  humidity  measurements,  it  was  possible  to 
infer  the  local  ftraction  evaporated  from  a  knowledge  of  initial  conditions 
from  the  measured  local  static  pressure,  and  from  a  one-dimensional  anal¬ 
ysis  involving  the  equations  of  continuity,  momentum,  and  energy  between 
the  inlet  and  the  section  in  question  (refs.  9  and  8).  Three  additional 
assumptions  were  necessary,  to  wit 5  (i)  the  droplet  moves  at  the  gas  speed 
(ii)  the  value  for  the  peeudo  friction  factor  used  in  the  analysis  may  be 
taken  at  an  average  value  for  f  of  0.004,  and  (ill)  the  droplet  temp¬ 
erature  is  130®F.  With  these  assumptions  the  curves  marked  by  the  full 
points  in  Figure  25  were  obtained.  Because  of  assumiptlon  (i)  above,  the 
discontinuity  analysis  can  be  fairly  accurate  only  at  considerable  dis¬ 
tance  from  the  inlet.  Where  the  aesumption  Ve/V^*^  is  a  fairly 
good  one,  the  validity  of  the  calculated  humidity  will  depend  on  the 
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choice  of  the  value  used  for  f,  V/ith  conventional  values  of  the  friction 
factor,  this  discontinuity  analysis  leads  to  values  of  humidity  which  are 
too  high.  A  comparison  of  the  direct  humidity  data  with  those  inferred 
from  the  pressure  measurements  indicates  similar  trends  and  orders  of  mag¬ 
nitude,  with  the  quantitative  agreement  improving  as  S  increases. 

Total  Hater  Evaporated.  Figure  26  shov;B  how  the  specific  humidity 
noar  the  end  of  the  duct  (  «£  =  675"  .  )  varies  with  v/ator-air  ratio 

(  )  and  initial  Mach  Kumter  (  ). 

For  a  fixed  value  of  M,  ,  the  exit  specific  humidity  increases 
with  A  because  the  total  amount  of  heat  transfer  surface  in  the  drop¬ 
let  cloud  is  proportional  to  .  As  12^  proceeds  above  0.3,  how¬ 

ever,  the  exit  specific  humidity  tends  to  level  off  owing  to  the  stream 
having  become  nearly  saturated  at  a  specific  humidity  of  about  0.21,  such 
that  the  exit  humidity  can  scarcely  increase,  no  matter  how  much  water  is 
fed  in. 

How  ijq)ortant  a  role  is  played  by  drop  size  is  shown  strikingly  in 
Fig.  26,  where  it  may  be  seen  that  for  fixed  values  of 
-Clj,  >-  the  exit  humidity  increases  as  Mj  increases.  This  despite 
the  fact  tliat  the  residence  time  of  the  droplets  is  leas  at  the  higher 
Mach  Numbers.  As  shown  in  Section  4*6,  the  length  required  for  coagplete 
evaporation  varies  approximately  in  inverse  proportion  to  the  initial  Mach 
Number;  in  a  given  length  the  total  amount  evaporated  increases  as  Mj 
increases. 

Wall  fyiction  Factor,  , In  order  to  interpret  the  small-scale  exper¬ 
imental  data  and  to  predict  the  performance  of  large-scale  units,  knowledge 
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of  tho  effectivo  friction  factor  f  is  desirable .  The  value  of  f  nay  well  be 
dlffarent  fton  that  for  ordinary  fully-developed  pipe,  flows  for  several 
reasons;  (i)  the  water  cloud  probably  tends  to  make  the  velocity  profile 
more  uniform;  (11)  evaporation  alters  the  density  and  viscosity  profiles; 
(ill)  near  the  pipe  inlet  the  boundary  layer  surely  does  not  fill  the  pipe; 
(Iv)  in  much  of  tho  evaporation  section  the  static  pressure  is  rising, 
which  is  opposite  to  the  falling  pressure  characteristic  of  simple  pipe 
flow;  (v)  lSQ}lngement  cf  water  on  the  walls,  crashing  of  droplets  into 
the  boundary  layer,  and  tearing  of  water  off  the  walls,  all  Dost  have  some 
Influence  on  the  skin-friction  coefficient*  'Hiese  effects  are  separate 
and  complex,  and  for  the  present  it  is  necessary  to  simply  lump  them  all 
into  an  "effective"  ftrlctlon  factor. 

Values  of  the  effective  friction  factor  under  various  clroumstanees 
were  evaluated  from  three  types  of  experiments; 

(i)  Dry  runs  at  low  and  high  temperatures. 

(ii)  Vet  runs  at  low  temperature  and  therefore  with  no  evaporation* 
(ill)  Vet  runs  at  high  temperature,  with  considerable  evaporation* 
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For  type  (i),  the  friction  factor  could  be  determined  from  the  Fhnnoollne 
analysis  (ref,  2,  C3iapt.  6).  In  addition,  a  friction  factor  ^ 
types  could  be  estimated  from  a  "dlBcontlnultyiL 
conservation  laws  of  mass,  momentum,  an^ 
exit  of  the  constant-area  dnc*- 
pressure  and  h\^ 
example^  th| 

Deta'^ 
ar'i 


rT 


<•  S4 


deg.R 

1500 

1500 

520 

520 

1500 

1500 

4 

Pq^,  psia 

U.8 

10.9 

14.7 

14.7 

U.6 

14.6 

^0 

0 

0 

0 

0.30 

0.15 

0.25 

Me 

0.597 

0.597 

0.595 

0.500 

0.500 

0.500 

i”*  (discontinuity) 
analysis 

0.0034 

0.0034 

0.0032 

0.0035 

0.0040 

0.0046 

£  (Fanno-line) 

0.0038 

0.0038 

0.0035 

mm  m  ^ 

175,000 

128^000 

605,000 

502,000 

154,000 

154,000 

0.0041 

0.0044 

0.0032 

0.0033 

0.0043 

0.0043 

The  last  Wo  rows  rspreient,  rsspeotlysly,  the  pipe  KeynoUs  Donlier  asd 
the  usual  friction  factor  tor  fully-developed  turbulent  pipe  flow  at  that 

Reynolds  Number e 

In  the  dry  rune  the  dlBContlnulty  analysis  (vhloh  inrolves  cerUlm 
assumptlona  leading  to  aasll  error8_ln  thfl  nO?1ITMUtf>on  of  ^ 
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•i,L  Conmariaon  of  Theory  and  Experiment 

We  turn  not/  to  a  compai-ison  of  the  experimental  results  of  Section  5.3 
with  those  computed  using  the  theoretical  method  of  !3ection  4.7. 

Figure  27  shows,  in  addition  to  the  observed  curves  of  Figure  25 
for  ’  ,  the  corresponding  theoretical  curves  computed  with 

f  =»  0.005  and  the  Huklyama  and  Tanasawa  drop  sise  of  19  miorons  as  well 
as  with  slightly  smaller  drop  si’ses.  The  measured  and  theoretical  curves 
are  of  the  same  general  character,  and  the  agreement  is  not  disappointing 
when  one  considers  the  broad  sweep  of  physical  phenomena  which  the  theory 
attempts  to  model* 

In  Jig.  28  are  shown  the  theoretical  curves  corresponding  to  the  ob¬ 
served  critical  curves  of  ilg.  19.  The  theoretical  critical  Mach  Number 
of  0.74  con^ares  favorably  with  the  measured  value  of  0.78,  and  the  theory 
gives  a  gratifyingly  accurate  picture  of  the  different  regimes  of  oper¬ 
ation,  including  the  interposition  of  shocks  with  subsequent  subsonic  opei> 
ation. 

These  results  speak  well  for  the  theory  as  a  means  of  estimating  perform- 


e-scale  unit,  and  as  a  basis  for  determining  the  optlaum  area 
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Wity  dl.tril»tl»n  1.  q»lt.  InwnsltlT.  to  cl»m«es  In 
comparatively  large  wall  ftrictlon  In  the  snall-acale  constant 

t,„d,  to  fore,  tHo  ^  ">•  “* 

autn^tio.  in  .  ur*0-«al.  nnlt  »,«lr..  o,*.*..  to 

„oo^l.h  tto  ««  r.»at.  Ihi.  »nn.  that  th,  c-^l-  to  h.  n^a 

^10.  1.  hotn.™  th.  oon.t«>t^«  nnlt  a  • 

in  >a.loh  th.  a-  1.  con-tantOy  n»hxo.a  In  »ch  a  nay  a.  to  «!«  th. 
„„  ,^h  »-h.r  at  «ch  vain,  of  *  a.  In  th.  a«U-a.al.  nnlt. 

^.a  «.  th.  for-goln,  n^aonln,.  v.  nay  prnc«a  a.  foUona.  Inaa- 
.^h  a.  all  th.  ph.no»n..  .xo.pt  that  of  nail  fiction,  ar.  allk.  In  ^ 
0.n,.  oon.ia.»a,  n.  nay  (naln,  Tahl.  Ill)  th,  p«f»r«n»  .«f- 

flcient  as 
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(Bq.  13) 


^  y  inolnaoa  th,  eff.cta  of  aU  tom,  nthav  than  that  of  naU  friction, 
la  th.  for  th,  ca,,,  oonald,r«,.  »on,  as.»lng  f  „nal  to 
1.0010,  I  any  h,  conpntoa  for  th,  anall-al,  nnlt,  and  th,  foUonlng 
aeanlta,  corr„po.^ind  to  th,  pointa  of  ohUlnod. 
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TABLE  VI  -  EXTRAPOUTION  Og 
smaLL-SCMI^FQRMM^ 

Tq^,  deg.R 

1200 

1500 

1800 

0.777 

0.780 

0.783 

Poi» 

14.7 

14.7 

U.7 

w.ilb/aeo 

B 

0.755 

0.677 

0.612 

f 

0.0040 

0.0040 

0.0040 

k 

1.35 

1.35 

1.35 

Dfinohes 

2.13 

2.13 

2,13 

Vd 

34.5 

34.5 

34.5 

^2. 

0.17 

0.24 

0.31 

[  <P03- 

“Poi^/poi*^  ]  nut; 

-0.23 

-0.16 

-0.13 

y 

"0.15 

0.22 

0.25 

Vd  "  20 

0.95 

1.00 

1.02 

I/D  »  15 

0.99 

1.03 

1.05 

l/D  »  10 

1.02 

1.0? 

1.08 

ft/ 

1/D  -  5 

1.05 

1.10 

1.12 
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The  results  under  the  dashed  line  in  the  table  represent  extrapolated 
performance  for  the  large-scale  vinlt. 

To  gain  some  estimate  of  slzsi  wa  obsez*ve  iVom  the  pressure  and  • 
humidity  data  that  a  length  of  about  five  feet  is  required.  A  dia¬ 
meter  of  one  foot  would  correspond  to  an  air  flow  of  about  25  Ib/sec, 
or  that  of  a  gas  turbine  of  about  1500  net  horsepower.  With  5» 
a  stagnation  pressure  reduction  at  the  turbine  exit  of  about  lOjS  is  seen 
from  the  table  above  to  be  possible  when  T^j^  *  1500*R. 

There  seems  little  doubt  that  the  large-scale  performance  estim¬ 
ated  above  may  be  achieved  as  a  minimum.  It  seems  equally  likely  that 
considerably  better  performance  is  attainable,  the  reasons  being  as 
follows: 

(1)  Only  a  strange  coincidence  could  assure  that  the  Mach  Number 
dletrlbutions  involved  in  the  foregoing  calculations  are  the  optimum 
for  each  combination  of  Pq^,  Tq^,X2q,  and  f/Oj^.  What  is  more,  in 
large-scale  units  the  optimum  values  of  and  length  will  surely 

differ  from  those  in  smaller  units, 

(li)  The  friction  factor  should  be  less  for  large  units  than  for 
small  because  of  a  more  favorable  pipe  Reynolds  Number. 

(lii)  Little  attention  was  given  to  diffuser  design  in  the  small- 
scale  unit,  and  gains  may  very  likely  be  won  by  careful  attention  to  this 
element  of  the  apparatus* 

(iv)  Undesirable  effects  associated  with  water  near  the  walls 


be  minimized  in  larger  units. 


-  89 

Therefore,  the  extrapolatione  should  be  considered  conservative, 
and  the  extrapolated  performance  of  large-scale  units  shovra  in  the 
table  above  as  assured  levels,  with  a  considerable  margin  of  ir  cve- 
msnt  possible* 

As  explained  before,  the  extrapolated  performance  refers  to  a  large- 
scale  \mlt  having  an  area  variation  that  will  produce  the  same  Mach  Number 
distribution  as  in  the  small-scale  unit.  Considerable  practical  Inqjort- 
ance  is  therefore  attached  to  the  question  of  how  much  area  variation  la 
required  to  assure  the  attainment  of  extrapolated  performance.  This  may 
be  estimated  from  the  influence  coefficients  in  the  first  line  of  Table  III. 
To  obtain  the  same  Mach  Number  variation  in  a  variable-area  unit  without 
flriction  as  In  a  constant-area  unit  with  friction,  it  is  only  neoessary 
to  set  equal  the  terms  in  dh/k  and  in  4fds/D  appearing  In  the  equation 
for  dM^/M^,  remembering  that  the  terms  Involving  dw/w  will  be  the  same 
in  both  instances.  Thus  we  obtain 

clA  _  Ji!il  Af  ^ 

A  z  ^  D 

as  the  sought-for  criterion.  Lacking  precise  information  as  to  the  Mach 
Number  distribution  in  the  small-scale  unit,  we  may  only  compute  the 
over-all  change  in  area  In  terms  of  some  average  value  of  .  Thus  we 
may  write 

‘  #  S  -if/-. 

^  2  P 


Taking  lt»  f.35  ,  “f  ■  0.004-  ,  L/D^S^.?  ,  and  using  a  value  of 

M*  ^  0.7  f  we  obtain  =  —O.ZB  ,  This  means  that 

the  exit  area  must  be  approximately  three- fourths  of  the  inlet  area  in 
a  unit  so  large  that  friction  is  negligible*  Whether  this  is  the  best 
area  reduction  for  a  large-scale  unit  remains  problematical. 

The  actual  manner  of  area  variation  required  to  obtain  the  Kach  Number 
distribution  of  the  small-scale  unit  cannot  of  course  be  found  from  this 
simplified  calculation. 

6*  MBDUm-SCALB  AgtOTHERMOPRESSOR  EgERlMEMT 

At  tlie  time  of  this  writing,  there  is  under  preliminary  test  a 
medium-scale  Aerothermopressor,  11  Inches  in  diameter,  with  a  maximum 
available  air  flow  of  25  Ib/sec.  Pbr  these  preliminary  tests  the  evap¬ 
oration  section  is  of  constant  outside  diameter  and  is  thoroughly  instrum¬ 
ented.  Area  variations  will  subsequently  be  provided  by  means  of  internal 
tapered  plugs?  this  general  approach  is  thought  best  suited  to  determin¬ 
ing  optimum  area  variation  most  quickly  and  economically. 

In  Section  4  it  was  shown  that  the  greater  portion  of  the  evapora¬ 
tion  section  should  be  of  diminishing  area  in  order  to  establish  a  desir¬ 
able  Mach  Number  variation.  In  the  small-scale  constant-area  unit,  the 
comparatively  large  frictional  effects  accomplished  the  same  Mach  Num¬ 
ber  variation  as  that  which  would  result  from  ah  ai'ea  reduction  of  about  ons- 
fourth  in  the  absence  of  friction.  Therefore,  the  Mach  Number  distribu¬ 
tion  of  the  medium-scale  constant-area  unit  is  less  favorable  than  that 
of  the  small-scale  constant-area  unit.  Accordingly  it  may  be  expected  that 
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the  medluin-scale  constant-area  unit  will  not  have  as  good  performanca 
as  the  extrapolated  performance  shown  In  Table  VI. 

Despite  the  hemdicaps  accompanying  this  medium-scale  constant-area 
unit  mentioned  above,  the  preliminary  teats  have  already  yielded  the 
gratifying  result  of  a  net  rise  in  stagnation  pressure  of  about  356  ft*om 
the  upstream  stagnation  section  to  the  diffuser  exit.  What  Is  ln^jortant 
here  is  the  positive  proof  of  the  Aerothermopreesor  principle,  rather 
than  the  Bnallness  of  the  ne<crlee  achieved  under  unfavorable  circuie- 
stances.  there  is  evei7  reason  to  believa  that  substantial  improveioents 
lie  ahead,  ready  to  be  gained  by  intelligent  control  of  the  cross-sectional 
area, 

7.  ASSOCIATED  PROBLEMS 

The  research  and  development  discussed  in  this  report  has  brou^t 
to  the  fore  a  number  of  technical  problems  which  are  of  special  signifi¬ 
cance  to  the  Aerothermopressor.  These  are  named  below,  and  the  steps 
looking  toward  their  solution  outlined, 

7tl  ■  Determination  of  Drop  vSize 

The  weakest  link  between  theory  and  experiment  is  currently  a  know¬ 
ledge  of  the  drop  size  produced  during  atomization.  Two  methods  of  meas¬ 
uring  mean  droplet  size  are  under  development,  one  depending  on  optical 
scattering  of  a  light  beam  (ref,  10),  the  second  depending  on  the  pressure 
gradient  near  the  hose  of  an  impact  tube  (refs.  11,  12,  13), 

7.2  Atomization 

Nothing  would  better  a.-3sure  good  perforJiianco  of  the  Aerothermoprossoi- 
than  a  moans  for  producing  extromoly  small  droplets,  liicporimontc  wo 


bsing  planned  t,o  determine  the  effect  of  gas  speed  and  gas  properties 
on  droplet  size,  and  to  determine  whether  tricks  of  injector  design  can 
lead  to  smaller  droplets, 

1,3  Measurement  of  the  Properties  of  a  Droplet-Iaden  Stream 

An  understanding  of  what  happens  inside  the  Aerothermopressor  re¬ 
quires  the  noasurement  of  sufficient  properties  to  establish  the  point 
properties  of  the  gas  and  the  droplet  cloud.  Various  typos  of  devices 
have  been  developed  or  are  under  development  for  this  purpose: 

(i)  A  probe  for  drawing  a  sample  of  the  gas  phase  only  has  been 
successfully  developed  (ref.  9). 

fii)  Probe  designs  for  accurately  measuring  the  stagnation  pressure 
of  the  gas  phase  only  have  been  devised  and  calibrated  (ref.  13). 

(ill)  Techniques  for  sampling  the  liquid  flow  per  unit  area  are 
under  development  (refs.  13,14}. 

(Iv)  A  scheme  for  determining  the  droplet  velocity  has  received 
some  experimental  attention  and  appears  promising  (ref.  13). 

8.  CONCLUSIONS 

1,  The  wide  variety  of  gas  dynamics  phenomena  observable  within 
the  Aerothermopressor  are  explainable  in  terras  of  heat  transfer,  evap¬ 
oration,  droplet  drag,  wall  friction,  and  cross-sectional  area  variation. 

2,  Experiment  and  theory  both  show  that  the  Aerothermopressor  can  • 
produce  a  net  rise  in  stagnation  pressure  with  mass  flows  of  about  2  Ib/seo 
or  greater. 

3»  A  Net  rise  in  stagnation  pressure  has  already  bean  won  in  a 


25  lb/»©e  \mlt,  and  In  this  Bi*e  tha  theory  Buggeats  that  rlaoB  of 
perhapa  205<  a^e  achievable.  The  latter  figure  would  lead  to  improve- 
menta  in  both  fuel  econoay  and  power  capacity  of  a  aimple  gaa  turbine 
plant  of  about  20^^ 

4.  The  theoretical  analyala  of  the  Aerothermopresaor  gives  a  gen¬ 
erally  correct  picture  of  ita  performance. 

5.  With  high  subsonic  Mach  Nurabera  at  the  inlet,  the  required 
length  of  evaporation  aeotion  is  about  5  feet,  and  for  an  initial  stag¬ 
nation  temperature  of  1500®R,  the  optlawn  water-air  ratio  la  in  the  neigh¬ 
borhood  of  0.25. 
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10 «  APPENDIX  A  -  DERIVATION  OF  INFLUENCE  COEFFICIENTS 

The  analysis  is  most  conveniently  carried  out  by  casting  the  defin¬ 
itions  and  physical  laws  applicable  to  the  Infinitesimal  control  volume 
of  Fig.  4  into  the  form  of  logarithmic  dlfferentia],a.  The  gas  mixture 
is  assumed  to  obey  the  perfect  gas  law  and  the  Oibbs-Ealton  law,  but  with 
variable  specific  heats. 

Definition  of  Mach  Number,^  For  a  perfect  gas. 


M 


z 


(A-l) 


M  -  zH.  X.  dvj  ik  olT 

T?  ■  ^1/  ■  T  T 


(A-2) 


Perfect-gas  Law^ 

f  = 


dp  ^  _dW 

p  "  T  W 


(A-3) 

u  '  «• 


(A-4) 


1 


-  95  - 


ConBcrvation  of  Mass.  With  w 
phase,  an  A  the  net  cross-sectional 

■ur  = 

^  +  dA  4,^  (a-6) 

7tr  u,-  ^  ^  A  V 

Momentum  Theorem,  Taking  account  of  the  pressure  and  shear  forces 
acting  on  the  control  volume,  and  of  the  momentum  fluxes  entering  and 
leaving,  we  have(Pigrf  4b) 


the  mass  rate  of  flow  of  the  ga^ 
area  for  gas  flow, 

(A-5) 


fdA  "(KpA)  -  ZjrrDeb:  ^  (yr-\-durX^-^<^)  -f- 


Introducing  the  definition  of  local  skin-friction  coefficient. 


simplifying,  and  rearranging  in  dimensionless  form  with  the  help  of  the 

IcpM^ 


relation 


,  the  momontm  equation  bocomfi 


Considering  the  heat  transft*  dQ  ftroa  the  walls 
and  the  fluxes  of  enthalpY  and  kinetic  energy  (Fig.  4c )#  the  first  lav 
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of  thermodynamics  is  written  as 

W?  -  d{\fQ  ’hdivrj^  (L(wi^}  +  d(wVV2)  ^d(yrf,Vjz) 

For  a  mixture  of  perfect  gases  following  the  Glbbs-Dalton  law, 

>!areo7er,  from  continuity, 

T»r  ^ 

doj"  *  ctuJ^  =  -'d'UJ^ 


Combining  these,  simplifying,  and  rearranging,  the  energy  equation  becomes 


Formation  of  Influence  Coefficients.  Equations  A-2,  A-4,  A-6, 

A-7,  and  A-8  conatitute  five  simultaneous  linear  algebraic  equations 
among  twelve  differential  variables  i  dMVM*,  ^pjp  * 
iMpy  dW/W,  dur/nr,  dlc/lc^  dAA,  4fdz/D, 

Therefore,  seven  of  these  may  be  chosen  as  Independent,  and  the  remain¬ 
ing  five  may  be  computed  from  the  values  of  these  seven  Independent  vai’- 
iables.  Fco*  the  latter,  the  last  seven  are  selected,  inasmuch  as  these 
represent  most  directly  the  influence;  of  the  duct  geometry  and  of  the 
droplet  acceleration  and  evaporation  on  the  gas  properties.  Solving 
the  equations  simultaneously  in  the  manner  outlined,  the  influence 
coefficients  appearing  in  the  first  four  lines  of  Table  III  are  obtained. 

Gas  Stagnation  Temperature.  The  stagnation  temperature  of  the  gas  phase 
alone  is  defined  as  the  temperature  which  would  be  reached  if  the  gas 
phase  were  steadily  and  adlabatlcally  decelerated  to  zero  speed  without 
interacting  with  the  liquid.  Placing  its  exact  variation  within  the  scheme 
of  the  influence  coefficients  Is  complex  because  k  varies  during  the  docol- 
eration.  Since  the  effect  of’ this  variation  in  k  is  very  small,  and  for 
the  sake  of  obtaining  results  simple  and  instructive  in  form,  a  pseudo- 
stagnation  temperature  is  defined  as  above,  except  that  the  gas  is  assumed 
to  have  constant  k  during  the  deceleration.  Then  (ref.  2), 


T,/T 


1  + 


(A-9) 


±  <i!/f  j. 

T  ,  1-1  ,  M* 


1  +  Mm"  m’ 


¥  dt 
TTEv  T 


(A-IO) 
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from  which  the  fifth  line  of  Table  III  is  found. 

Gae  Stagnation  Preaenre.  The  stagnation  pressure  of  the  gas  phase 
alone  is  defined  as  the  pressure  which  would  be  reached  If  the  gas  phase 
were  steadily,  reversibly,  and  adiabatically  decelerated  to  sero  speed 
without  Interacting  with  the  liquid.  As  above,  we  de^Jne  a  pseudo- 
stagnation  pressure  based  on  the  assumption  of  unchanged  k  during  the 
deceleration.  Then  (ref.2), 

B/P  ‘  (*-u) 


from  which  the  seventh  line  of  Table  III  is  found. 

Migbure  Stagnation  Temperature,  ppr  certain  purposes  a  consideration 
of  the  stagnation  properties  of  the  mixture  of  gas  and  droplets  is  more 
revealing  than  that  of  the  gas  phase  alone.  Here  the  question  arises  of 
how  to  define  the  mixture  stagnation  properties;  that  is,  what  are  the 
interactions  between  the  droplets  and  gas  during  the  fictitioua  decelera¬ 
tions?  After  reflection  on  the  various  possibilities,  what  seem  most 
useful  are  definitions  which  differ  from  the  previous  definitions  only 
In  that  they  take  account  of  the  kinetic  energy  of  the  droplets.  Therefore, 
the  mixture  stagnation  temperature  is  defined  as  the  gas  temperature  which 


would  be  reached  if  the  mixture  were  eteadily  and  adiabatically  decel¬ 
erated  to  aero  speed  with  only  mechanical  (i«e.y  no  heat  transfer  or  evap¬ 
oration)  interactions  between  the  gas  and  droplets^  Then^  again  using  the 
notion  of  a  pseudo  stagnation  temperature  for  constant  k  during  the  decel¬ 
eration, 


which,  after  some  manipulation  becomes 


or,  in  differential  form. 


7.  oC  -w  ZoC  ^  ^ 


(A-U) 


from  which  the  sixth  lino  of  Table  III  is  found. 
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Mixture  Stafmatlon  Preannra.  Analogously  with  tho  mixture  stag¬ 
nation  temperature,  ve  define  the  mixture  stagnation  pressure  as  tho 


gas  pressure  which  would  be  reached  if  tho  mixture  were  steadily,  revors 
ibly,  and  adiabatically  decelerated  to  zero  speed  with  only  mochanloal 
interactions  between  the  gas  and  droplets.  For  the  deceleration  to  ho 
reversible,  the  gas  alone  must  first  change  speed  until  it  is  at  tho 


same  speed  as  the  droplets;  then  the  two  phases  must  decelerate  together 
so  slowly  that  the  difference  in  velocity  between  the  phases  is  at  each 
point  vanishingly  small.  Since  the  gas  undergoes  isentroplc  changes 
during  the  defined  process,  the  pseudo-stagnation  pressure  obtained 
by  ignoring  variations  in  k  during  the  deceleration  is  found  with  the 
help  of  Eq.  A-13  to  be 


3. 

T 


T. 


4+i. 

f  M 


(A-13) 


(A. 


from  which  the  last  line  of  Table  III  is  found. 

■11.  APPENDIX  B-AWALYSIS  FOR  APPROXIMATB  OPTIMZi^TTQU 
PROCEDURE  IN  REGIME  II 

With  the  aim  of  obtaining  an  explicit  and  simple  design  procedirro, 
the  following  assumptions  are  made  which  are  approximately  valid  in  Reglxie  II: 
(i)  y  =  Ij  (li)  Tg  is  constant;  and  (ill)  the  Influence  of  the 


t 


i 

I 

i‘ 

% 

I 

\ 

I 

\ 

I 

V’ 


s 


-  101  - 


I  V 


dkA  t«nii  on  stagnation  prassoro  la  negligible,  tbr  tlie  purpoee 
te  which  the  present  analysis  Is  later  put.  a  conelderatlon  of  the 

gas  stagnation  presenre  f,  leads  down  a  blind  alley  (because  the 

enalysls  would  Indicate  epeclouelj  that  could  locally  be  nado 

to  rlsa  s»et  rapidly  by  extreme  decelerations  which  would  profit 
through  the  momentum  of  the  droplets),  and  so  wo  shall  deal  with 

p 

the  mixture  stagnation  pressure  'o  •  . 

n .  1  AppT9yJ,"*^Q  f SS,^  /R 

Using  the  assumptions  listed  above »  Table  III  yields  for 
the  expression 


T  -2 


t+^ 

It«*J 

'^Bq.B-l) 

w 


•f 


kM^ 


Ji  .  dW 


+ 


In  what  follows, we  shall  express  dW  and  In  terms  of  dur  , 
and  thus  we  shall  obtain  an  expression  showing  the  amount  of  stagna- 
tion  pressure  rise  per  unit  of  evaporation. 

J1.2  ihtrodnQtion  of  Soeciflo  Humidity; 

Using  the  definitions  of  local  specific  humidity {  CtJ  ^  ) 

and  of  initial  water-air  ratio  (IZ^sttg/aC  ),  we  may  formthe 
following  equations  relating  the  masa-now  terms  of  Bq.  B-1  to  the 
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apacific  humidity j 


TiC+w;;- 


T^-f  K 


I  +  il» 
I  ^fO 


(£ci*B*2) 


duT  _  duliS-  _  (d  (^r/vfQ 
^  1+  (W/Wl) 


duL 

i+cj 


(Eq.  B-3) 


Change  of  Droplet  Pjamatar 

Aasumlng  the  existence  of  a  droplet  cloud  of  uniform  sl^e,  and  that 
neither  further  atomlaation  or  agglomeration  ocoura  In  Regime  II,  the 
number  of  droplets  per  unit  time  passing  any  aeotlon  will  be  unehanged. 
The  mass  rate  of  water  flow  will  therefore  be  proportional  to  be  eube 
of  the  droplet  diameter: 


(i) 


_  J2^-a) 


(Eq.  B-4) 


Differentiating,  we  get 


^  ^  ^  I  y  >^•5) 


11.4  Evaluation  of  dV/ 


The  molecular  weight  of  the  mixture  of  gases  isorelat^d  to  the 
molecular  weights  of  the  component  species  by 


i+CJ 


(Eq»B-6) 


i+w(VilWr) 


Honoe  Eq.  5  becomes,  after  rearrangement 


(Eq.B.7) 

w  (vc/K-)+w  1+'*^ 

U.5  Heat.Tranafqir  tO  ProgM 

Since,  by  assumption,  the  droplet  temperature  is  unaltered,  and  the 
droplet  and  gas  have  the  same  speed,  the  amount  of  heat  transferred  between 
gas  and  droplet  during  the  time  interval  dt  may  be  equated  to  the  en¬ 
thalpy  change  of  that  portion  of  the  droplet  which  evaporates.  Thus 


The  distance  traveled  by  the  droplet  during  the  time 


At 


is 


ds  «=  = 


Combining  these  relations,  employing  Eqa.  B-4  and  B-5,  and  expressing 


the  h0at-trans'''i3r  coefficient  by  the  Nussolt  Number  through  the  definition 

Ha  -  f  \ 


iZiVM 


}\.fy  WorklnE  Formria  .tSL 

Setting  Equ.  B-2,&-3,  B--?,  and  B-8  into  Eq.  B-1,  Bimpllfying 

and  reai'ranging,  we  finally  obtain 

idl  ^  -  Trr^r— 1 

2  -f/p  si 

y 

{Bq.B-9) 


11,7  Calculation  of  T 


Con8ld«ring  ths  .nw-gy  eq^tion  ralatli^  tta  propertlM  at  tha 
Aarotherimpreaaor  Inlet  to  those  at  any  section  nhere  y  -  1,  v»  may  wlte 


.1 


y  Iloy 


kK  ■*■ 


(Eq.B-10) 


The  enthalpy  changes  may  be  closely  approximated  by 


and  th(  term  y^lt  may  be  expressed-  as 

t  . 


Introducing  these  into  the  energy  equation,  the  local  stagnation 
temperature  may  be  solved  for  as 


T* 


Op 


(Bq.S-il) 


11,8  Estimate  of  Nuaselt  Humber 

The  value  of  Nu  depends  upon  the  relative  Reynolds  Number  between 
droplet  and  gas  stream  (ref .7).  If  y  were  exactly  unity,  this  Reynolds 
Number  would  bo  zero,  and  Nu  would  be  2.00,  Although  y  may  be  so  close 
to  unity  that  it  may  be  so  assumed  for  the  purpose  of  Eq.  B-1,  its  deviation 

from  unity  may  generate  a  sufficiently  high  Reynolds  Number  to  give  Kussol 
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Numbers  substantially  greater  than  2.00.  To  find  the  local  value' of  Ku, 
therefore,  would  require  integrating  the  entire  set  of  governing  equations 
of  the  Aerothermopressor  after  the  manner  of  Section  4*7*  But  this  is 
exactly  what  we  wish  to  avoid  in  the  present  simplified  design  method. 
Fortunately,  the  numerical  integrations  of  the  type  of  Section  4*7  show 
that  in  Regime  II,  the  value  of  Nu  lies  between  2  and  4,  thus  setting  upper 
and  lower  limits  on  Nu  which  are  not  excessively  apart. 


Suppose  that  in  a  given  case  the  conditions  at  the  Aerothermopressor 
inlet,  i.e.,  \  ,  and  ,  are  known.  In  Regime  II 

the  value  of  may  be  set  equal  to  the  average  wet-bulb  temperature, 

which,  for  the  rsmge  of  pressures  and  temperatures  encountered,  is  generally 
not  far  from  140®P.  A  particular  value  of  Nu  between  2  and  4  is  also 
selected.  In  addition,  we  choose  an  appropriate  average  value  of  f/D, 
nfiglecting  the  variations  of  diameter,  inasmuch  as  these  amount  only  to 
soma  10  or  20^  in  the  range  of  interest. 

Then,  for  each  combination  of  M  and  dO  ,  Eqa*  B-'9  and  B-10  allow 
the  value  of  (j1^  to  be  calculsted.  Thus  it  is  possible 

to  construct  curves  of  /l^du)  versus  M,  with  ^  as  a 


parameter,  as  exemplified  by  the  chart  of  Fig.  7a. 


^2.  APPENDIX  G  -  governing  PHYSICAL  EQUATIONS  tOR  DROPLETS 

Consider  a  single  droplet  of  diameter  d,  velocity  '/g  ,  and 

temperature  TJ  ,  Journeying  in  a  gas  stream  of  velocity  V'  ,  temper, 
ature  T  ,  pressure  f  ,  and  specific  humidity  W  ,  As  the  mean 
dlsUnce  between  droplets  Is  of  the  order  of  10  droplet  dismeters  or 
more,  mutual  interactions  between  the  droplets  are  Ignored,  and  the 
gaseous  medium  surrounding  each  droplet  is  imagined  as  infinite  in 

extent  * 

12.1  Drool  at  Diameter  yfl.  Stream  HumMjt; 

Since  the  droplets  are  uniform  and  incompressible,  and  of  un¬ 
changed  number  per  unit  time,  the  mass  rate  of  liquid  flow  at  any  sec¬ 
tion  is  proportional  to  the  cube  of  the  droplet  diameter.  Thus  the 
local  humidity,  which  is  a  measure  of  the  losg.  in  liquid  mass  flow, 

is  given  by 


3 


or,  after  differentiation,  the  change  in  droplet  diameter  is  related 
tc  the  change  in  humidity  by 


M  .-i  ^ 

a:  ■  3  (d/a/  (d/ir. 


(Eq. 
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In  evaluating  the  driving  potential,  is  taken  aa  the 

reciprocal  of  the  opecific  volume  of  saturated  water  vapor  at  the 
temperature  TJ  )  is  computed  from  the  perfect  gas  lav  of 

mixtures  as 

^  "P  CJ _  {Eq*  C-4) 

/V«  "  ^ 

Heat  Transfer 

Fi-om  a  consideration  of  the  first  law  of  thermodynamics,  the 
rate  of  heai  transfer  to  the  droplet  is  equal  to  the  sum  of  the  rate 
of  enthalpy  change  of  both  the  evaporated  and  unevaporated  water.  Thus 


which  simplifies  to 


§  -  ^  * 


3£  tU 


(Eq.  C-5a) 


where  is  the  latent  heat  at  the  tomporature ■  1^  •  Using 

Eqs.  G-1  and  C-2,  this  may  be  written  in  the  alternate  form 


dT, 

dx  ^  /’v<-/5~] 


(Eki.C-5b) 
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The  heat-tranafer  coeffioiant  Hj.  has  also  been  taken  Area 


the  correlation  of  Rans  and  Marshall  (ref*  7), 

«  2+0.6  CT?t)  (l3e/) 


(Eq.  C-6) 


where 


{Bq.C-7) 


Bnploylng  the  conventional  definition  of  drag<  coefficient  for  blunt 
bodies  like  spheres,  vs  may  equate  the  net  foroe  on  the  particle  to  the 
product  of  nass  and  acceleration.  Thus 


from  which 


(Eq.  (V8) 


In  the  numerical  calculations  the  well-established  correlation 
between  and  Rey  for  solid  spheres  has  been  approximated  tgr 


i 


r 
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Cjj  =>  24/Rey 

for 

0  < 

Rey 

< 

1 

Cjj  =  24/Rey^/^ 

for 

1  < 

Rey 

390  (Eq.C-9) 

Cjj  =  0.45 

for 

390  i 

Rey 

00 

12 »5  Displacement  of  Droplet 

The  distance  along  the  duct  traveled  by  the  droplet  cloud 

in  the  time  dt  is  given  by 


d-j/dt  «  \ 


(Eq.  C-1 


I 


13.  APPENDIX  D  ~  APPROXIMATE  DISCONTIIs^ITY  ANALYSIS 

13 tl  Governing  Integral  Equations 

Let  section  1  be  the  inlet  plane  and  section  2  the  exit  plane  of 
the  Aerothermopressor,  and  let  the  entire  water  input  be  evaporated  com¬ 
pletely  at  section  2  (it  is  known  from  other  considerations  that  excess 
water  has  little  effect  on  performance).  Then  we  may  write  the  follow¬ 
ing  conservation  equations  between  sections  1  and  2. 

Energy; 

(^<1^  D  +  A  (4^"  +  A, 


X 


*Og  ^ 

T 

'c>( 


L-l  ^ 


1+ 


kri  2 


=  1 


(Eq.  D-1 
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i; 


ilbMU 


Continuity; 


Definition  of  Mach  Romber; 


Eouatlon  of  State; 


(£q.  D-2) 


(Eq,  D-3) 


Momentum  Equation; 


\A4  .  T[ 
v^i  Tz'f, 


(£q*D^/ 


Omitting  wall  friction  for  the  present, 


-  -i^A^  +  [-pc/A  “ 

To  compute  the  definite  Integral  requires  a  detailed  knowledge  of  the 
drop  history,  which  we  are  here  avoiding.  As  shown  by  Wadleigh  (ref.9)» 
the  results  are  insensitive  to  the  way  in  which  A  varies  with  p,  and 
so  we  have  arbitrarily  chosen  a  linear  variation; 


dfl  s  conshni  - 
dn  Az~^t 


•t 


I 


I 


Evaluating  the  Integral  in  the  momentuni  equation,  and  sijDplifying, 
ve  get 


S) 


1  -  (i+A.) 


(Bq.D-S) 


Stagnation  Preaeurei  Ihe  well-known  iaentropio  relation  glvea 


k 

(1+ 


(Eq.D-6a) 


As  shown  in  ref.  2,  p.  95,  the  isentropic  stagnation  pressure 
may  also  bs  expressed  approximately,  at  least  in  tenns  up  to  order 
as 


f 


Retaining  only  the  terms  of  order 


,  we  may  form 


_  _L  {Eq.D-6b) 


where  the  subscript  infinity  here  refers  to  infinite  duct  size,  for 
which  the  wall  friction  is  indeed  zero. 


nrat>Qrder  Solution.  Ei^llclt  algebraio  solution  of  the  foregoing 

M.^ 


set  of  equations  for  seams  In^sslble 


For  low  Mach  Numbers,  bowerer,  such  a  solution  may  be  approximated.  Hav¬ 
ing  in  mind  the  a8au]iq)tlon  of  low  Mach  Numbers,  we  may  according  to 
Eq.  D-5  set  Vzlf]  ^  ^  ^  ^  further 

1.11,  =  i/t;  I"  a.-  D-l.  Algebraic  combination  of 
Eqs.  D-2,  D.3  and  D-4  then  yields 


=  0+X2.#)g 


(Eq.D.7) 


Where 


(£q.D-6} 


In  addition,  algelralo  combination  of  Eqs.  D-2,  D-4  and  D-5  produces 


1+4^  L 

A, 


i  -("i+ry-is 


(Bq.D-9) 


The  calculation  procedure  is  now  as  follows,  for  given  values  of 
Ml,  Tj,  12  o»  and  Ag/A^  : 

(i)  T2  Is  calculated  l^om  the  low  Mach  Number  approximation  for 
Bq.  D-l. 

(H) 


Is  calculated  from  Eq.  D-7. 


1 
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(iil)  (p.-R)ATiMr  is  calculated  fTon  Eq.  D-9« 

(iv)  /|  yZ.  is  calculated  from  Eq*  D-6b. 

W.0  /n'R'A 

E?cact  Discontinuity  Solution*  If  more  accurate  results  than  those 
obtainable  from  the  low  Mach  Number  analysis  are  desired,  trial-and-error 
solution  of  liqs.  D-1,  D-2,  D-3»  D-4i  D-5  and  D-6a  is  necessary,  the  approx¬ 
imate  solution  serving  as  a  useful  starting  point.  ■  First  we  derive  sev¬ 
eral  required  formulas. 

Eliminating  ^/V^l  fi'om  Eq.  D-5  with  the  help  of  Eqs.  D-2  and  D-4, 


wo  obtain 


%-f, 


i  - 


Ci+I2o)-g _ 


(Eq.D-10) 


Solxiticn  of  this  for 


yields 


U  _  I  + 

V  1-1-^- 

'■  A, 


1  - 


(\+<X)z 


Then,  enploying  Eqc.  D-10,  D-2  and  D-4,  tho  value  of  ^2. ^*^2  I 


be  found  irorQ  Eq.  D-3  as 


(Eq.It.ll) 


IK  ^ 

|i+ 

1  ‘+^ 

Ai  -- 


\  i'  '..j  <■  >.  • 


The  exact  expression  for  the  stagnation  prefiisiare  ratio  may  be 


found  ftrom  Eq,  D-6a  as 


(Bq.D-13) 


in  which  the  value  of  is  that  giveh  hy  Eq,  Ml. 

Assuming  that  values  of  >  T,  ^  and  ^2^1 

been  selected,  the  numerical  calculations  may  be  carried  out  in  the 
foUoidng  order: 

(i)  Guess  using  the  low  Mach  Number  analysis  as  a  guide. 

(ii)  Compute  the  corresponding  Tg  from  Bq.  D-1. 

(iil)  Evaluate  Zft-om  Eq.  D-8, 

(iv)  Check  whether  the  guessed  value  of  is  correct  by 
computing  1/^  from  Eq.  D-12. 

(v)  Using  the  value  of  fTom  Eq*  D<»12,  repeat  steps  (1) 

to  (iv)  until  satisfactory  convergence  is  obtained.  Only 
two  trials  are  usually  necessary. 

(vi)  Calculate  (p^-fj  )/fJ  from  Bq.  D-10. 

(vii)  Solve  Eq.  D-13  for  ('fJt /foi)^. 

v.f facts  of  Wall  Friction  and  Diffuser  Loss.  We  may  now  take  account 
of  wall  friction  in  a  manner  which  is  in  the  same  vein, as  the  approxima¬ 
tions  already  made.  Taking  note  of  the  influence  coefficients  of  Table  II 
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the  change  In  naeoclated  with  wall  ftlctlon  mj  b.  eapreeaed  aa 


(Eq.D^U) 


where 


Besides  wall 


is  the  aritlmetic  moan  of  lc,M,  and 
friction,  there  is  a  loss  in  the  diffuser  which  may 


bo  expressed  as  a  percent  loss  of  stagnation  pressure  in  terms  of  a 
loss-coefficient  6  , 


(Eq.D-15) 


Linearly  superposing  the  losses  due  to  wall  friction  and  the  dif. 
fuser  on  the  previous  results,  we  obtain 


(Eq.D-16)  . 
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CAPTIONS  TO  FIGURES 


Illustrates  one-dimensional  process  on  which  Eq.  1  is  based. 

Fig,  2  Application  of  the  Aerothermopressor  to  the  simple  gas  turbine 

(a)  At  exit  of  turbine 

(b)  Between  burner  and  ttarbine 

Fig,  3  Effect  of  Aerothermopressor  on  performance  of  gas  turbine  plant 

(a)  Qycle  arrangement 

(b)  Fuel  economy  and  power  capacity  for  favorable  cycje 

parameters 

(c)  Fuel  economy  and  power  capacity  for  unfavorable  cycle 

parameters 

Fig.  4  Illustrates  one-dlraensionnl  analysis  of  Appendix  A. 

(a)  Nomonclature  and  control  volume. 

(b)  Forces  and  momentum  fluxes  for  control  volume 

(c)  Fluxes  of  enthalpy  and  kinetic  energy  for  control  volume 
Schematic  representation  of  bokavior  of  sxibsonic,  constant-area 

Aerothermopressor 

Continuous  pannage  tlirough  the  speed  of  sound 

(a)  From  subsonic  to  £'i..^;ersonic 

(b)  From  supersonic  to  subsonic 

Fig.  7  .Approximate  optimi nation  method,  based  on  analysis  of  Appendix  D. 

(a)  T;/]3ical  curves  of  fractional  rate  of  chan  pa  of  mixU/ro 

stagnation  prossur.e  \7ith  speoifio  humidity  as  a  function 
of  Kach  Number,  with  local  spncific  humidity  es  a  para¬ 
meter,  and  for  constant  values  of  ,  Mu,  i/D, 

(b)  Optimum  curve  of  Mach  Humber  vs.  humidity,  corresponding 
to  "design  curve"  of  fig.  7a,  Dashed  curve,  sb'ovrs  eoiupar- 
ative  position  of  optimum  ‘curve  for  moro  favorablo  values 

of  parorr.olorvS, 


I 


(6)  Opt  latum  curve  of  iZkioi  VB.  humidity,  correapondlng  to 
"design  curve"  of  Fig.  7a.  Dashed  curve  shows  comparative  position 
of  optimum  curve  for  more  favorable  values  of  parameters. 

(d)  Illustrates  optimal  Mach  Number  variation  of  entire  Aerothermopreasor. 

Fie,  8  Typical  theoretical  calculations  for  subsonic,  constant-area  operation. 

Fig,  9  Typical  theoretical  calculations  for  supersonic  constant-area  operation.  ^ 

Wtii  o.»aiwo„..  ;  T„m500-«;  =  ^/D,-o.oo4ft, 

12,-0.15  ;  1J|  ” 6?o'R;  lilAr .  4)“  • 

a,.  10  Th.orrtioC  of  .ioe  on  conrt«.t-«r«.  »*.onlo 

operation.  Initial  conditions? 

M,  -  0.5  ;  T.,  -  I500"R  ;  -  W.7  pni.  ;  J2,-  o.z  ; 

T^i-sso’r;  i-ij-tV- 

Fie,  U  Theoretical  effect  of  initial  drop  sisa  on  small-scale,  constant- 
area,  subsonic  operation.  Initial  conditions?  ^ 

M,  -  0.5 ;  »  I50(f R ;  -  14.7  pViOf ,  f /D,  -  0.0226  ff; 

12, -0.2  ;  0  i 

Fig,  12  Theoretical  effect  of  initial  water-air  ratio  on  small-scale, 
constant-area  subsonic  operation.  Initial  conditions; 

M,  -0.5-  ;  =  l500“Ri 


J-3  Theoretical  curve  of  constant-area  operation  at  various  subsonic 

initial  Mach  Nuabers,  illustrating  continuous  passage  at  critical 

initial  Mach  N\unber  through  speed  of  sounds  and  normal  shocks  on 
on  supersonic  critical  branch.  Corresponding  curves  of  static  pres¬ 
sure  are  shown  in  Fig.  27. 

Fig.  14  Illustrates  theoretical  improvement  hy  control  of  cross-sectional 
area.  Initial  conditions: 

M,  -  0.727  ;  «  I500'r;  -f)!  -  I4.7  pwo/ ;  f/D,  =  0.004  ft*; 

TJ,  *  550*R,;  ;  ^0  *  doj*  13.4^  ;  12^ ’0.25. 

Solid  curves  represent  subsonic  and  supersonic  critical  operation. 
Dashed  curves  show  variable-area  operation  with  Mach  Number  distri¬ 
bution  shown  (beginning  at  x  =  0.30),  follow^  by  6^-lncluded  angle 
diffuser  (beginning  at  x  *  0.54)*  Dash-dot  curves  show  constant- 
area  operation  on  supersonic  critical  with  shock  at  x  »  0.30, 

^5  Typical  results  of  calculations  based  on  "low-velocity"  discontin¬ 
uity  analysis.  Dashed  curve  shows  results  of  "exact"  disccatinuity 
analysis* 

16  Arrangement  of  small-scale  experimental  unit. 

Fig.  17  Small-scale  experimental  results,  showing  compcurlson  between  water 
Injection  without  evaporation  (  ~  530®R.  and  with 

evaporation  {  ■=•  I500®R.  ), 

Small-scale  experimental  results,  showing  effect  of  water-air  ratio, 
riSi-lS  Small-scale  experimental  results,  showing  effect  of  initial  Ifech 
Number,  continuous  passage  through  speed  of  sound,  and  shocks  in 
'  '.per sonic  regime. 

^0  Illustrate.^  schematically  shock  positions  for  stable  and  unstable 


operation. 


Fig,  21  anall-Boale  experimental  results,  showing  effects  on  over-all 
performance  of  water-air  ratio,  initial  stagnation  temperature, 
and  Initial  stagnation  pressure.  Each  point  is  for  the  initial 
Mach  Number  or  shock  position  of  beet  performance. 

Fig,  22  Small-scale  experimental  results,  showing  effects  on  over-all 
performance  of  Initial  stagnation  temperature  and  pressure. 

Each  point  is  for  the  water-air  ratio  of  best  performance. 

The  right-hand  point  of  each  curve  represents  the  choking  Mach 
Number,  and  refers  to  the  shock  location  of  best  performance. 

Fig,  23  Small-scale  experimental  results,  shewing  effect  of  stagnation- 
pressure  level  on  static  pressure  distribution. 

Fig,  24  Snail-scale  experimenta],  results,  showing  humidity  profiles  at 

several  longitudinal  cross  sections, 

M,  »  0.  ^02  ]  *  0.  PyO  ;  '1^,  =  j  fj,,  =  14.6  psw . 

Fig,  25  Small-scale  experimental  results,  showing  fraction  evaporated 

vs,  longitudinal  distance.  Open  points  refer  to  direct  humidity 
measurements.  Closed  points  are  inferred  IVom  static  pressure  data. 

Fig,  26  Small-scale  experimental  results,  shoving  exit  humidity  (  2  «675  inches  ) 
as  function  of  water-air  ratio  and  initial  b5ach  Number. 

Fig,  27  Comparison  of  small-scale  experimental  results  (dashed  curves) 

with  theoretical  calculations  (solid  curves)  for  completely  subsonic 
operation. 

Fig,  28  Comparison  of  small-scale  experimental  results  (dashed  curves) 

with  theoretical  calculations  (solid  curves)  for  subsonic- supersonic 
operation. 
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